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Abstract 

Observations  of  current  velocity,  temperature,  salinity  and  pressure  from  a  2-year 
moored  array  deployment  and  four  hydrographic  cruises  conducted  by  the  United 
States  Southern  Ocean  GLOBEC  program  on  the  western  Antarctic  Peninsula  conti¬ 
nental  shelf  are  used  to  characterize  the  ocean  circulation  and  its  connection  to  fresh 
water  and  heat  fluxes  on  the  shelf.  Mean  velocities  on  the  shelf  are  of  the  order  of 
5  cm/s  or  less.  Tidal  motions  are  dominated  by  the  M2  and  S2  semi-diurnal  tides 
and  the  0\  and  K\  diurnal  tides,  although  the  tidal  velocities  are  typically  less  than 
2  cm/s.  Near-inertial  motions  are  relatively  large,  with  current  velocities  as  high  as 
26  cm/s.  It  is  shown  that  Marguerite  Trough,  a  large  bathymetric  feature  connecting 
the  shelf-break  to  Marguerite  Bay,  plays  a  critical  role  in  determining  the  circulation. 
The  mean  flow  is  strongly  steered  in  the  along-slope  direction,  and  the  tidal  currents 
also  show  increasing  current  polarization  at  depth  in  Marguerite  Trough.  At  time- 
scales  of  5  to  20  days,  the  observations  show  bottom-intensified  motion  in  Marguerite 
Trough  consistent  with  bottom-trapped  topographic  Rossby  waves.  The  subtidal  cir¬ 
culation  in  the  trough  has  a  significant  wind-driven  component  in  Marguerite  Trough, 
with  downwelling-favorable  winds  forcing  cross-shelf  flow  on  the  northern  side  of  the 
trough  and  along  the  shore  on  the  outer  shelf.  Upwelling-favorable  winds  force  roughly 
the  opposite  circulation.  The  cyclonic  circulation  on  the  trough  helps  advect  blobs 
of  salty,  warm  and  nutrient-rich  water  across  the  shelf.  These  intrusions  are  small 
(^4  km)  and  frequent  (4  events/month).  Also,  the  Antarctic  Peninsula  Coastal  Cur¬ 
rent  (APCC),  a  coastal  buoyant  current  which  is  described  for  the  first  time  here.  The 
APCC  is  a  seasonal  current  which  is  only  present  during  the  ice-free  season  and  is 
forced  by  freshwater  fluxes  associated  with  large  glacier  melt  and  precipitation  rates 
in  the  region. 
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Title:  Emeritus  Scientist,  Woods  Hole  Oceanographic  Institution 
Thesis  Supervisor:  W.  Brechner  Owens 

Title:  Senior  Scientist,  Woods  Hole  Oceanographic  Institution 
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Chapter  1 


Introduction 


Polar  continental  shelves  remain  one  of  the  least  studied  ocean  ecosystems  in  the 
world.  This  is  partly  due  to  the  fact  that  it  involves  collection  of  observations  in 
distant  places,  under  generally  harsh  weather  conditions  and  with  the  added  compli¬ 
cation  of  a  permanent  or  seasonal  ice  cover.  Many  polar  shelves  host,  however,  very 
diverse  and  highly  productive  biological  communities,  and  their  unique  physical  set¬ 
ting  and  characteristics  makes  them  distinctive  from  their  mid-latitude  counterparts. 
Moreover,  many  studies  have  shown  these  regions  to  be  especially  sensitive  to  recent 
changes  in  the  global  climate. 

The  continental  shelf  off  the  west  coast  of  the  Antarctic  Peninsula  (Figure  1-1), 
which  extends  for  some  1200  km  from  the  southern  tip  of  Alexander  Island  to  the  tip 
of  the  peninsula,  is  one  such  system.  It  hosts  a  rich  biological  community,  including 
large  colonies  of  penguins,  whales  and  other  mammals,  whose  main  food  source  is 
Antarctic  Krill  (Euphasia  super  bn).  Marguerite  Bay  and  the  surrounding  shelf  play  a 
central  role  in  the  life  cycle  of  krill  (Hofmann  et  al.,  2002.  2004),  and  it.  is  believed  that 
a  unique  combination  of  physical  and  biological  factors  in  th('  bay  and  its  surroundings 
explain  the  success  of  the  local  krill  population,  which  also  sustains  other  biological 
communities  in  the  southern  ocean,  especially  east  of  the  Drake  Passage  (Murphy 
et  al..  1998.  2007).  Studies  of  the  regional  climate  strongly  suggest  that  the  Antarctic 
Peninsula  has  one  of  the  fastest  changing  environments  in  the  world,  with  recent 
warming  of  3°C'  in  the  surface  air  temperature  and  1°C  in  the  ocean  since  1950 
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(Meredith  and  King.  2005),  fast  retreating  glaciers  (Cook  et  al.,  2005)  and  gradually 
shorter  winter  seasons  (Harangozo,  2006).  For  convenience,  we  will  call  the  region 
around  Marguerite  Bay.  from  the  southern  tip  of  Alexander  Island  on  the  south  to 
the  northern  tip  of  Adelaide  Island,  the  western  Antarctic  Peninsula  (wAP)  shelf. 
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Figure  1-1:  The  Antarctic  Peninsula  continental  shelf.  The  SO  GLOBEC  region  is 
outlined  by  the  dashed  black  line,  and  the  main  geographical  features  are  indicated. 


The  ocean  circulation  and  dynamics  of  the  wAP  shelf  are,  despite  its  relevance 
for  the  life  cycle  of  Krill,  poorly  understood.  The  observations  of  water  properties 
are  historically  mostly  limited  to  the  summer  months,  as  observations  during  the 
winter  are  impeded  by  weather  conditions  and  a  seasonal  ice  cover  which  forms  over 
most  of  the  wAP  shelf  from  July  to  October  (Stammerjohn  and  Smith,  1996).  As 
for  observations  of  current  velocity  on  the  shelf,  there  were  limited  to  the  region  of 
the  shelf  north  Adelaide  Island,  particularly  in  Bransfield  Strait  (Niiler  et  ah,  1990; 
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Huntley  and  Niiler.  1995:  Zhou  et  al.,  2002).  To  my  knowledge,  there  are  no  historical 
observations  of  current  velocity  on  the  shelf  of  the  study  region. 

The  need  to  understand  the  physical  processes  and  their  connection  to  biological 
dynamics  in  the  wAP  region  led  to  the  United  States  Southern  Ocean  Global  Ocean 
Ecosystems  Dynamics  (U.S.  SO  GLOBEC)  program,  a  multi-disciplinary  effort  whose 
field  work  component  included  several  broad-scale  cruises  during  the  fall  and  winter 
of  2001  and  2002.  the  deployment  of  the  first  long-term  mooring  array  on  the  shelf, 
and  a  host  of  biogeochemical  observations. 

Both  the  earlier  studies  and  initial  analysis  of  the  U.S.  SO  GLOBEC  hydrographic 
data  suggest  several  important  component  of  the  wAP  shelf  circulation  and  dynamics. 
To  motivate  this  thesis,  a  brief  summary  of  the  regional  circulation  and  previous 
studies  of  the  circulation  are  given  below. 


1.1  The  Circulation  on  the  wAP  Shelf  from  Prior 
Studies 

The  main  component  of  the  large  scale  circulation  in  the  Southern  Ocean  and  offshore 
of  the  wAP  shelf  is  the  Antarctic  Circumpolar  Current  (ACC  -  Figure  1-2),  which 
flows  towards  the  northeast,  roughly  parallel  to  the  direction  of  the  shelf-break.  The 
ACC  features  three  distinct  fronts:  the  Subantarctic  Front  (SAF).  the  Polar  Front 
(PF)  and  the  Southern  ACC  Front  (SACCF)  (Orsi  et  al..  1995).  Deflected  northwards 
by  the  presence  of  the  Antarctic  Peninsula,  the  ACC  moves  closer  to  the  continental 
margin,  and  the  Southern  ACC  Front  is  found  at  80  to  200  km  from  the  shelfbreak 
(Figure  1-2).  Notice,  however,  that  the  Southern  ACC  Front  is  not  considered  a 
shelf-break  front. 

Prior  to  U.S.  SO  GLOBEC,  only  geostrophic  velocities  from  hydrographic  sec¬ 
tions  had  been  obtained  on  the  shelf  surrounding  Marguerite  Bay,  the  study  region 
(Hofmann  and  Klinck,  1998;  Smith  et  al.,  1999).  Mid-depth  dynamic  height  fields 
referenced  at  400  db,  corresponding  to  the  mean  baroelinic,  geostrophic  component 
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of  the  circulation,  suggested  the  presence  of  a  large  cyclonic  circulation  off  Adelaide 
Island,  with  currents  of  the  order  of  1  cm/s  flowing  southward  near  the  shore  and 
northeastwards  along  the  shelf-break  (Figure  1-3).  The  flow  towards  the  north-east 
on  the  outer  branch  of  this  gyre  was  thought  to  be  “attached”  to  the  ACC  (Smith 
et  al..  1999).  This  shelf-scale  gyre  was  thought  to  provide  a  retention  mechanism  for 
plankton  on  the  shelf. 

Although  the  processes  determining  the  observed  hydrographic  structure  are  still 
being  investigated,  it  is  clear  that  transport  of  warm,  salty,  and  nutrient-rich  water 
from  the  Antarctic  Circumpolar  Current  (ACC)  onto  the  shelf  is  a  critical  factor 
(Klinck,  1998;  Smith  et  al.,  1999).  These  intrusions  are  also  believed  to  influence  the 
seasonal  ice  thickness  over  the  shelf  during  the  winter  (Smith  and  Klinck,  2002). 

The  above  studies  provided  much  of  the  motivation  for  the  physical  oceanogra¬ 
phy  component  conducted  as  part  of  SO  GLOBEC,  and  of  this  thesis.  The  present 
study  includes  a  basic  description  of  the  flow  and  the  main  scales  of  variability.  It 
also  explores  two  major  components  of  the  circulation  according  to  prior  knowledge: 
the  presence  of  a  southwards  circulation  which  forms  part  of  a  larger  gyre,  and  the 
intrusions  of  warm  water  onto  the  shelf. 

1.2  Thesis  Outline 

The  thesis  is  organized  as  follows:  Chapter  2  summarizes  the  SO  GLOBEC  obser¬ 
vations  relevant  to  this  thesis,  Chapter  3  deals  with  the  mean  flow,  tidal  motions, 
near-inertial  variability,  bottom-trapped  topographic  waves  and  the  wind-driven  re¬ 
sponse  on  the  shelf.  Chapter  4  explores  the  characteristics  of  warm  water  intrusions 
onto  the  shelf,  and  Chapter  5  provides  the  first  description  of  a  seasonal  buoyant  cur¬ 
rent  flowing  along  the  wAP  coast.  Finally,  Chapter  6  provides  the  overall  conclusions 
of  this  work. 
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Figure  1-2:  Location  of  the  Subantarctic  (SAF),  Polar  (PF),  Southern  ACC  (SACCF) 
fronts  and  the  Southern  Boundary  of  the  ACC  (SBdy).  Data  from  Orsi  et  al.  (1995). 
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Figure  1-3:  Dynamic  height  at  200  db  (referenced  at  400  db)  for  the  shelf  west  of  the 
Antarctic  Peninsula  (from  Smith  et  ah,  1999) 
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Chapter  2 


The  U.S.  Southern  Ocean 
GLOBEC  data  set  and  study  area 


The  physical  oceanographic  component  of  the  U.S.  Southern  Ocean  GLOBEC  (SO 
GLOBEC)  Program  included  the  deployment  of  three  mooring  lines,  four  broad-scale 
hydrographic  cruises  and  deployment  of  surface  drifters  and  subsurface  floats.  The 
following  is  a  brief  description  of  the  data  sets  used  in  this  thesis.  The  moored 
array,  automatic  weather  stations  and  the  bathymetric  survey  information  presented 
here  were  mostly  extracted  from  the  Southern  Ocean  GLOBEC  moored  array  and 
automated  weather  station  data  report  (Moffat  et  al.,  2005),  where  a  more  detailed 
description  is  given. 

2.1  The  mooring  array 

The  WHOI  SO  GLOBEC  mooring  locations  are  shown  in  Figure  2-1.  The  2001-2002 
array  consisted  of  one  transect  (the  “A-line”)  of  three  moorings  (A1-A3)  deployed 
across  the  shelf  west  of  Adelaide  Island  and  a  second  transect  (the  “B-line”)  of  three 
moorings  (B1-B3)  deployed  west  of  the  mouth  of  Marguerite  Bay.  The  2002-2003 
array  consisted  of  a  L-shaped  transect  (the  “C-line”)  of  three  moorings  (C1-C3)  de¬ 
ployed  slightly  inshore  of  the  B-line. 

The  A-line  was  designed  to  investigate  the  shelf-scale  gyre,  with  Al  in  the  coastal 
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southwestward  flow  and  A3  on  the  outer  shelf  in  the  northeastward  flow.  A2  was 
located  between  A1  and  A2  on  the  eastern  side  of  Marguerite  Trough,  where  deep 
intrusions  of  Upper  Circumpolar  Deep  Water  (UCDW)  were  thought  to  occur.  The 
B-line  was  deployed  across  the  mouth  of  Marguerite  Bay,  with  B1  and  B3  on  the  east 
and  west  flanks  and  B2  in  the  center.  This  line  was  designed  to  look  for  clockwise  flow 
around  the  basin,  the  southward  continuation  of  the  surface  coastal  current  towards 
the  mouth  of  Marguerite  Bay,  and  the  filling  and  flushing  of  the  deepest  part  of  the 
basin.  Based  in  part  on  the  preliminary  results  of  surface  drifters  deployed  during 
2001  austral  summer/fall  plus  a  desire  to  learn  more  about  flow  and  water  properties 
within  the  mouth  of  Marguerite  Bay,  the  C-line  array  was  set  across  Marguerite 
Trough,  with  Cl  and  C2  aligned  with  the  mouth  of  the  Bay  and  C2  and  C3  on  the 
east  and  west  sides  of  the  trough.  C3  was  deployed  at  the  B2  site  to  provide  a  two-year 
time-series  at  this  central  location. 


Table  2.1:  SO  GLOBEC  mooring  locations,  bottom  depths  and  dates  of  deployment 
and  recovery. 


Mooring 

Latitude 

Longitude 

Bottom  Depth 

Deployment 

Recovery 

A1 

67  1.134’ 

69  1.217’ 

509 

3/26/01 

2/13/02 

A2 

66  1.883’ 

70  0.683’ 

561 

3/30/01 

2/13/02 

A3 

66  5.002’ 

70  9.991’ 

480 

3/31/01 

2/13/02 

B1 

67  6.890’ 

69  4.398’ 

444 

3/30/01 

Not  Recovered 

B2 

68  6.091’ 

70  1.675’ 

811 

3/29/01 

2/14/02 

B3 

68  5.345’ 

70  9.853’ 

447 

3/29/01 

2/14/02 

Cl 

68  2.940’ 

69  1.790’ 

432 

2/18/02 

2/22/03 

C2 

68  3.331’ 

70  1.730’ 

859 

2/19/02 

2/26/03 

C3 

68  6.006’ 

70  1.799’ 

806 

2/21/02 

2/26/03 

The  moorings  were  designed  to  measure  the  following  physical  variables  at  com¬ 
mon  depths:  (a)  temperature  at  50,  100,  150,  200,  250,  400  m  and  close  to  the 
bottom  using  SeaBird  (SBE)  temperature  sensors;  (b)  conductivity  (salinity)  at  50, 
100,  250,  and  400  m  using  SBE  conductivity  sensors;  and  (c)  currents  at  250  and 
400  m  using  Vector  Averaging  Current  Meters  (VACMs).  To  obtain  vertical  profiles 
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of  currents  and  acoustic  back-scatter  intensity  in  the  upper  water  column  where  krill 
are  commonly  found  (Lascar a  et  al.,  1999),  300-kHz  broadband  RDI  Acoustic  Doppler 
Current  Profilers  (ADCPs)  would  be  deployed  at  100  m  on  all  moorings  except  the 
B2/C3  moorings  which  would  feature  a  150-kHz  broadband  RD  Instruments  (RDI) 
ADCP  moored  at  250  m.  These  design  depths  were  chosen  based  on  regional  mean 
profiles  of  T  and  S  (Figure  2-2),  the  desire  to  better  sample  the  upper  water  column, 
the  idea  from  hydrography  that  the  deeper  flow  would  be  more  barotropic,  and  the 
instrumentation  and  mooring  resources  available. 

The  A-  and  B-line  recovery  and  C-line  deployment  were  conducted  on  the  R/V 
L.M.  Gould  SO  GLOBEC  mooring  cruise  LMG02-1A  (Beardsley  et  ah,  2002).  A  CTD 
cast  was  made  before  each  mooring  recovery  for  later  use  in  checking  temperature  and 
conductivity  sensor  performance  and  drift.  The  three  A-line  moorings  were  recovered 
successfully  while  the  two  southernmost  moorings,  B2  and  B3,  were  recovered  in 
patchy  ice,  with  damage  to  the  upper  instruments  due  to  passing  icebergs.  The  Bl 
mooring  was  lost  with  no  trace. 

The  final  position,  bottom  depth  and  deployment  and  recovery  times  for  each 
mooring  are  listed  in  Table  2.1.  The  basic  instrumentation  used  in  the  moored  array 
is  listed  in  Table  2.2.  Details  on  the  quality  control  tests  and  corrections  of  the  data 
can  be  found  in  (Moffat  et  al.,  2005). 


2.2  Automatic  Weather  Stations 

To  help  interpret  the  moored  array  data  and  monitor  the  surface  forcing  during  the  SO 
GLOBEC  field  program,  Raytheon  Polar  Services  Corporation  (RPSC)  technicians 
and  WHOI  investigators  deployed  Automatic  Weather  Stations  (AWSs)  on  Kirkwood 
(AWS-K)  and  Dismal  (AWS-D)  Islands,  located  near  the  center  of  Marguerite  Bay 
(Figure  2-1).  The  two  self-contained  AWSs,  provided  by  the  University  of  Wisconsin 
Antarctic  Automatic  Weather  Station  Project  (AWSP),  transmitted  vector-averaged 
wind  speed  and  direction,  air  temperature,  relative  humidity  and  barometric  pressure 
data  back  to  AWSP  via  ARGOS,  where  the  raw  data  were  decoded  and  sent  to  WHOI 
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Table  2.2:  Instruments  used  in  the  moorings.  The  acronyms 
listed  in  the  second  column  are  used  throughout  the  rest  of  this 
thesis.  The  physical  variables  are  T  =  Temperature;  C  =  Con¬ 
ductivity;  P  =  Pressure;  UV  =  Horizontal  velocity;  W  =  Vertical 
velocity;  B  =  Backscatter  intensity;  and  I  =  Ice  draft. 


Type 

Variables 

Sample  Rate  (s) 

RDI  ADCP  (WorkHorse) 

UV,  W,  B 

1800 

RDI  ADCP  (BroadBand) 

UV,  W,  B 

3600 

ALS  Ice  Profiler 

I 

50 

SBE  MicroCat 

T,  C 

150  1 

SBE  SeaCat 

T,C,P 

900 

SBE  TempRecorder 

T 

225 

VACM 

UV,  T,  Pi  ,C* 

900 

SBE  SeaGauge 

P,T,  P 

300 

1  Except  Al  at  309m  and  A3  at  250.  Both  sampled  at  300  s 
1  Only  when  indicated 


for  final  processing.  The  two  AWSs  were  deployed  in  May  2001  and  produced  high- 
quality  data  through  March  2003,  with  only  a  few  short  periods  of  zero  wind  speed 
or  constant  wind  direction  thought  due  to  anemometer  icing. 

At  Kirkland  Island,  the  AWS  (#8930)  site  is  on  a  slab  rock  shoulder  on  the 
northwestern  side  of  the  island,  at  an  estimated  height  above  sea  level  (rough  estimate) 
of  25  m.  The  winds  from  the  south  might  be  distorted  by  the  main  snow  cap,  but  the 
site  has  open  exposure  from  west  through  northeast  (Wiebe  et  al.,  2001).  The  other 
AWS  (#8932)  was  installed  on  a  small  rocky  island  east  of  Dismal  Island.  The  AWS 
was  installed  on  the  north  end  of  the  island,  with  open  exposure  to  the  west  through 
southeast,  at  roughly  12  m  above  sea  level  (Wiebe  et  al.,  2001). 

The  AWS  data  are  self-consistent,  and  compare  favorably  with  that  from  the 
British  research  station  Rothera  (R-sta)  located  on  the  southeastern  side  of  Adelaide 
Island  (Figure  2-1),  except  for  wind  speed  and  direction  at  Rothera  which  differ 
strongly  from  the  AWS  winds  due  to  the  rugged  coastal  topography  around  Rothera. 
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2.3  Hydrographic  Cruises 


During  2001  and  2002,  four  hydrographic  cruises  were  conducted  off  the  west  coast  of 
the  Antarctic  Peninsula  aboard  the  RV/IB  N.  B.  Palmer  as  part  of  the  SO  GLOBEC 
field  program  (Table  2.3).  These  broad-scale  cruises  attempted  to  cover  the  entire 
study  region,  and  provided  CTD  (conductivity  for  salinity,  temperature  and  pressure 
for  depth)  profile  observations  on  a  grid  with  a  nominal  cross-shelf  resolution  of 
20  km  (Figure  2-3).  The  fall  surveys  provided  a  larger  spatial  coverage  than  the 
winter  surveys,  in  which  weather  conditions  and  sea-ice  coverage  impeded  access  to 
Marguerite  Bay  and  other  nearshore  areas.  Sea  surface  temperature  (SST)  and  sea 
surface  salinity  (SSS)  as  well  as  wind  speed  and  direction  data  were  also  collected 
throughout  the  cruises  using  the  RV/IB  N.  B.  Palmer  continuous  sampling  systems. 
Details  of  the  processing  of  the  data  set  can  be  found  in  Wiebe  et  al.  (2001)  and 
Klinck  et  al.  (2004). 


Table  2.3:  Summary  of  the  cruise  dates  and  number  of  CTD  stations  conducted. 


Cruise 

Objective 

Start 

End 

#  CTD 

LMG01-03 

Mooring  Deployment 

03/18/01 

04/13/01 

6 

NBP01-03 

Fall  2001  Survey 

04/24/01 

06/05/01 

84 

NBP01-04 

Winter  2001  Survey 

07/24/01 

08/31/01 

70 

LMG02-1A 

Mooring  Recovery/Deployment 

02/06/02 

03/03/02 

9 

NBP02-02 

Fall  2002  Survey 

04/09/01 

05/21/01 

92 

NBP02-04 

Winter  2002  Survey 

07/31/02 

09/18/02 

102 

LMG03-02 

Mooring  Recovery 

02/12/03 

03/07/03 

4 

As  part  of  these  broad-scale  cruises,  conductivity,  temperature  and  pressure  data 
from  depths  of  roughly  15  m  to  250  m  were  also  collected  using  the  Bio-Optical 
Multi- frequency  Acoustical  and  Physical  Environmental  Recorder  (BIOMAPER-II), 
a  towed  undulating  multi-sensor  system  (Wiebe  et  al.,  2002).  Unfortunately,  the 
conductivity  sensor  had  an  unknown  malfunction  during  the  winter  of  2001  and  both 
the  fall  and  winter  2002  surveys  and  these  data  are  not  of  sufficient  quality  to  calculate 
salinity  reliably. 

Velocity  data  were  obtained  during  the  cruises  with  a  shipboard-mounted  RD 
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Instruments  (RDI)  150-kHz  narrow-band  Acoustic  Doppler  Current  Profiler  (ADCP) 
sampling  from  31  to  431  m.  The  ADCP  was  configured  with  fifty  eight-meter  bins, 
an  eight-meter  pulse  length,  and  five-minute  ensemble  averages.  During  most  of 
the  cruise,  the  ADCP  operated  in  bottom-tracking  mode.  The  data  processing  was 
initially  done  during  the  cruise  using  the  Common  Oceanographic  Data  Access  System 
(CODAS).  Final  processing  and  quality  control  were  done  by  J.  Hummon  (University 
of  Hawaii)  and  S.  Howard  (Earth  Systems  Research)  (Klinck  et  al.,  2004).  ADCP 
data  return  during  both  winter  cruises  was  greatly  reduced  while  the  RV/IB  N.  B. 
Palmer  was  moving  through  sea-ice. 

2.4  Shelf  Bathymetry 

Since  much  of  the  study  area  was  uncharted  prior  to  2001,  the  SO  GLOBEC  program 
collected  as  much  multibeam  bathymetric  data  as  possible  during  the  field  program 
using  the  RV/IB  N.B.  Palmer  SeaBeam  swath  mapping  system.  These  data  were 
then  merged  with  other  available  digital  along-track  and  multibeam  data  to  create  a 
high-resolution  digital  bathymetric  data  set  for  the  study  area  (Bolmer  et  al.,  2004). 

The  bathymetric  survey  resulted  in  a  new  high-resolution  bathymetry  which  re¬ 
solved  many  previously  unknown  features  and  areas  on  the  shelf’s  sea  floor  (Bolmer 
et  al.,  2004).  As  can  be  seen  in  Figure  2-4,  the  wAP  shelf  is  characterized  by  a 
relatively  deep  shelf  with  very  rough  topography.  The  shelf  is  approximately  400  m 
deep  on  average,  with  the  maximum  depth  of  «1600  m.  The  cross-shelf  bathymetric 
features  tend  to  be  deeper  next  to  the  coast  and  shallow  towards  the  outer  shelf. 

The  large  channel  connecting  the  shelf  break  to  Marguerite  Bay,  Marguerite 
Trough,  is  the  most  prominent  of  these  features  (Figure  2-5).  The  structure  of  this 
fjord-like  canyon  is  rather  complex:  at  the  shelf  break,  it  is  characterized  by  a  shallow, 
gentle  depression,  which  becomes  steeper,  narrower  and  deeper  onshore,  extending  for 
roughly  250  km  where  the  canyon  reaches  its  maximum  depth  of  ~1600  m  (the  deep¬ 
est  location  on  the  shelf),  to  become  shallower  again  closer  to  George  VI  Sound  (Fig. 
2-4). 
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Figure  2-1:  Map  of  the  western  Antarctic  Peninsula  shelf  region  showing  the  loca¬ 
tions  of  the  WHOI  SO  GLOBEC  moored  arrays  and  the  AMRC  Automated  Weather 
Stations.  The  A-  and  B-line  moorings  were  deployed  during  2001/02  and  the  C-line 
moorings  was  deployed  during  2002/03.  Bl,  shown  in  the  figure,  was  not  recovered. 
The  two  AWSs  at  Dismal  Island  (D-aws)  and  Kirkwood  Island  (K-aws)  were  deployed 
in  2001.  Also  shown  in  the  figure  is  the  location  of  the  British  Antarctic  Survey  station 
Rothera  (R-sta). 
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Figure  2-2:  Vertical  Profiles  of  potential  temperature,  salinity  and  potential  density 
obtained  during  2001  and  early  2002.  The  CTD  profiles  were  taken  near  the  A2 
mooring  during  fall  (blue),  winter  (green)  and  summer  (red).  At  these  cold  temper¬ 
atures,  density  is  determined  primarily  by  salinity.  Note  the  cooling  and  increase  in 
salinity  and  density  of  the  surface  50  m  from  summer  to  fall /winter;  the  deep  winter 
surface  mixed  layer;  the  summer  and  fall  temperature  minimum  layer  at  80/100  in, 
the  remnant  of  last  winter's  mixed  layer;  the  permanent  pycnocline  (120  to  150  m); 
and  the  deeper  UCDW.  The  dotted  lines  indicate  the  design  measurement  depths  in 
the  moored  array. 
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Figure  2-4:  Bathymetry  of  the  SO  GLOBEC  study  region.  A  detailed  view  of  the 
area  indicated  by  the  dashed  white  line  is  plotted  in  Figure  2-5  (Data  from  Bolmer 
et  al..  2004). 
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Figure  2-5:  Bathymetry  of  Marguerite  Trough.  The  figure  shows  a  top  view  of  the  bathymetry  in  and  around  the  trough  (top 
panel).  The  depth  along  the  thin  white  line  along  the  axis  is  plotted  in  the  middle  panel,  and  four  cross-trough  sections  arc 
plotted  in  the  bottom  panels  at  locations  indicated  by  thin  white  lines  in  the  top  panel  and  red  dotted  line  in  the  middle  panel. 
The  cross- isobath  sections  are  oriented  so  that  the  observer  is  looking  onshore  from  the  shelf-break  (tlu*  north  side  of  the  trough 
is  on  the  left),  and  are  plotted  with  the  same  horizontal  scale.  Data  from  Bolmer  et  al.  (2004). 
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Chapter  3 


Description  of  current  and 
pressure  variability  in  the  SO 
GLOBEC  study  area 

3.1  Introduction 

This  chapter  explores  general  aspects  of  the  circulation  on  the  wAP  shelf  using  the 
SO  GLOBEC  moored  array  dataset.  Figure  3-1  shows  the  kinetic  energy  spectrum 
of  the  318-day  long  horizontal  current  at  247  m  on  the  A2  mooring.  The  spectrum 
shows  the  typical  flow  variability  on  the  shelf.  At  high  frequencies  (1  day  or  shorter), 
there  are  peaks  at  the  primary  four  tidal  components  (See  section  3.3  for  details) 
and  a  large,  broad  peak  centered  at  the  inertial  frequency.  At  subtidal  time  scales 
(the  spectrum  includes  variability  for  period  shorter  than  60  days),  the  spectrum 
has  a  broad  significant  peak  at  10  days  or  so.  In  this  chapter  and  elsewhere  in 
this  thesis,  these  spectral  peaks  in  the  current  variability  are  referred  to  as  ’tidal’, 
’inertial’  or  ’subtidal’.  To  obtain  the  latter,  all  the  velocity  records  were  low-pass 
filtered  with  a  parabolic-linear  filter  with  a  33-h  half-amplitude  period  (Flagg  et  al., 
1976;  Limeburner  et  ah,  1983). 

In  the  following  sections,  I  describe  the  basic  structure  of  the  mean  flow,  the  tidal 
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motions,  the  inertial  frequency  variability,  and  explore  some  of  the  subtidal  variability 
in  Marguerite  Trough. 


Figure  3-1:  Kinetic  energy  spectrum  of  the  VACM  velocity  at  247  m  on  the  A2 
mooring.  The  inertial  frequency  and  the  principal  tidal  components  (Oi,  K i,  M2 
and  S2 )  are  indicated  with  vertical  lines,  as  well  as  the  half- amplitude  period  of  the 
low-pass  filter  used  to  generate  the  subtidal  time  series. 


3.2  Mean  and  Variability  of  the  Subtidal  Flow 

In  this  section,  basic  statistics  of  the  subtidal  current  velocity,  including  mean,  stan¬ 
dard  deviations,  and  principal  axes  are  presented.  For  the  basic  description,  nominal 
depths  of  50,  100,  250,  400  and  800  m1  were  chosen.  The  statistics  were  calculated 

'The  A1  mooring  was  located  in  a  deep  depression  50  m  deeper  that  its  design  depth,  and 
therefore  the  deeper  records  are  at  309  and  459  m. 
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for  both  the  original  and  rotated  records.  The  coordinate  system  used  for  the  rotated 
records  is  explained  below. 

3.2.1  A  Coordinate  System  for  the  Moored  Array 

The  means  and  principal  axes  of  the  subtidal  currents  (Figure  3-2)  show  large  spatial 
variability.  Only  the  A1  mooring  off  Adelaide  Island  appears  to  have  mean  and 
principal  axis  aligned  with  the  coast,  while  across  the  mouth  of  Marguerite  Bay  the 
direction  of  the  coastline  is  ambiguous.  The  influence  of  the  main  bathymetric  feature, 
Marguerite  Trough,  on  the  orientation  of  the  axes  of  maximum  variance  is  evident  at 
several  mooring  sites,  suggesting  that  the  flow  variability  is  mostly  along  isobaths  in 
the  trough. 

In  order  to  rotate  the  velocity  records,  the  principal  axes  of  the  vertically-averaged 
velocity  were  calculated  for  each  mooring.  In  the  case  of  the  nearly  co- located  B2/C3 
mooring  pair,  the  resulting  orientations  differed  by  2°.  The  orientation  calculated 
from  the  C3  mooring  data  was  used  for  both  of  them,  as  this  mooring  provided 
better  depth  coverage  of  the  water  column.  For  Cl,  the  deepest  instrument  was  at 
232  m  (bottom  depth  at  859  m),  so  that  the  records  were  averaged  down  to  this  depth 
only2.  Except  at  A1  and  B3,  where  the  direction  of  the  local  isobaths  is  ambiguous, 
the  rotation  to  the  direction  of  the  vertically-averaged  principal  axes  results  in  the 
rotated  velocity  being  roughly  in  a  coordinate  system  aligned  with  the  local  isobaths. 
The  sign  convention  was  chosen  such  that  x  is  positive  towards  the  shelf-break  along 
the  axis  of  the  trough. 

The  resulting  rotated  coordinate  systems  for  each  mooring  are  shown  in  Figure  3- 
3,  together  with  high  resolution  maps  of  the  bathymetry  around  each  mooring.  At  Al 
(bottom  depth  509  m),  there  is  a  shallow  sill  300  m  deep  only  a  few  kilometers  to  the 
south-west  of  the  mooring.  The  A2  and  A3  axes  are  aligned  with  the  local  isobaths, 
but  notice  that  the  slope  in  the  vicinity  of  A3  is  much  gentler  than  at  A2.  At  Cl, 
C2  and  B2/C3,  the  rotated  axes  are  also  clearly  defined  by  the  local  bathymetry. 

2 Assuming  a  barotropic  velocity  from  232  m  to  the  bottom  results  in  a  10°  difference  in  the 
principal  axis  orientation. 
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The  Cl  mooring  was  located  between  a  small  island  ~14  km  to  the  north-east  and  a 
small  sea- mount  to  the  south-west.  The  same  is  true  for  C2,  located  in  a  section  of 
Marguerite  Trough  where  a  small  sea  mount  is  located  in  the  middle  of  the  trough.  At 
B2/C3,  the  axes  are  aligned  with  the  main  axis  of  the  trough,  and  the  moorings  were 
flanked  by  a  very  shallow  (<200  m)  mount.  Finally,  the  B3  mooring  was  located  in  a 
lateral  extension  of  the  trough,  where  the  direction  of  the  local  isobaths  is  ambiguous. 


3.2.2  The  Mean  Circulation 

Basic  statistics  of  the  low-pass  filtered  SO  GLOBEC  moored  array  velocity  records 
are  given  in  Table  3.1.  The  uncertainties  presented  in  Table  3.1  correspond  to  the 
95%  confidence  intervals  computed  as 


u  ±  1.96 


<7 

(A’<5t/e)1/2’ 


(3.1) 


where  u  is  the  record  mean  velocity,  a  is  the  standard  deviation,  and  NSt/e  is  the  num¬ 
ber  of  independent  observations  of  a  time  series  of  length  NSt  and  integral  time-scale 
e,  calculated  from  the  autocorrelation  function  of  each  set  of  the  hourly  observations 
before  any  filtering  (Emery  and  Thompson,  1998).  For  the  velocity  records,  e  is  typi¬ 
cally  of  the  order  of  1  to  2  days.  Overall,  the  mean  velocities  in  the  study  region  are 
of  the  order  of  5  cm/s  or  less,  and  the  variability  of  the  subtidal  currents  is  generally 
of  the  order  of  the  mean  or  larger. 

The  mean  circulation  (Figure  3-2)  is  characterized  by  two  main  components.  First, 
near  the  coast,  the  surface  flow  (above  «100  m)  is  characterized  by  southward  flow 
along  Adelaide  Island  (at  Al)  and  into  Marguerite  Bay,  which  is  explained  by  the 
presence  of  a  coastal  current,  the  Antarctic  Peninsula  Coastal  Current  (APCC,  see 
Chapter  5).  The  mean  records  above  100  m  at  Cl  also  show  flow  into  Marguerite  Bay, 
but  the  magnitude  is  smaller,  and  the  velocity  records  do  not  show  the  consistent  flow 
seen  at  Al.  This  might  be  because  the  APCC  was  weaker  during  2002  or  because  the 
current  shows  larger  lateral  displacements  at  Cl.  In  the  southern  half  of  the  mouth 
of  the  bay,  the  flow  above  100  m  is  out  of  the  bay.  At  B3  the  flow  has  a  strong  mean 
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with  considerable  polarization  of  the  axes  of  variance.  At  B2/C3  the  mean  also  tends 
to  be  out  of  the  Bay,  although  it  is  much  weaker  than  at  B3,  located  some  25  km  to 
the  south-west.  This  indicates  that  the  flow  scales  are  small  across  the  mouth  of  the 
bay,  which  is  consistent  with  a  internal  radius  of  deformation  of  only  4  to  7  km. 

The  second  fundamental  component  of  the  mean  circulation  is  determined  by 
the  local  bathymetry.  As  is  characteristic  of  polar  shelves,  the  wAP  shelf  is  deep 
and  with  variations  of  the  order  of  the  mean  depth.  At  Al,  the  mean  flow  below 
250  m  is  relatively  weak  and  towards  the  shelf  break,  likely  a  result  of  blocking  of 
the  flow  by  the  presence  of  a  300-m  sill  just  a  few  kilometers  south  of  the  mooring 
site  (Figure  3-3).  The  rest  of  the  SO  GLOBEC  moored  array  was  located  in  or  near 
Marguerite  Trough  (Figure  3-2),  one  of  the  main  features  of  the  wAP  bathymetry, 
and  its  influence  is  clear  in  the  orientation  of  the  major  axis  of  the  subtidal  flow  at 
several  deep  locations  around  the  trough. 

On  the  eastern  side  of  the  trough3,  at  A2,  Cl  and  C2,  the  principal  axes  have 
relatively  strong  polarization  (major/minor  axes  ratios  from  1.4  at  232  m  at  Cl  to 
3.7  at  C2  at  810  m)  and  the  major  axis  is  aligned  with  the  local  isobaths.  At  A2, 
the  mean  shows  a  small  intensification  with  depth  (and  anti-clockwise  veering),  but 
the  velocity  is  barotropic  in  the  along-isobath  direction,  and  the  mean  values  are 
largest  mean  values  measured  at  or  below  250  m  in  the  entire  array.  The  mean  is 
directed  onshore  at  these  three  moorings,  although  the  along-isobath  component  is 
indistinguishable  from  zero  at  232  m  at  Cl  and  at  810  m  at  C2.  At  these  two  moorings, 
the  mean  (0.2-3. 1  cm/s)  is  weak  compared  to  the  subtidal  variability  (0.9-5. 5  cm/s). 

On  the  western  side  of  the  trough,  the  mean  flow  is  toward  the  shelf-break,  except 
at  247  m  at  B3,  where  the  flow  is  weak  and  toward  Marguerite  Bay.  At  A3,  although 
there  is  a  weak  but  significant  offshore  mean,  the  principal  axes  are  weakly  polarized 
(unlike  at  A2)  and  the  major  axis  is  oriented  in  the  cross-isobath  direction,  that  is, 
roughly  in  the  direction  of  the  coastline,  which  would  be  expected  in  the  absence  of 

3The  main  axis  of  Marguerite  Trough  is  roughly  aligned  in  the  north-northwest  direction.  For 
simplicity,  the  right  (left)  side  of  the  trough,  looking  from  the  coast  to  the  shelf  break,  will  be 
referred  to  as  the  ’eastern’  (western)  edge 
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Marguerite  Trough4.  At  B2/C3,  the  flow  is  weakly  polarized  in  the  top  100  m  but 
becomes  strongly  polarized  at  the  deepest  records,  which  were  taken  below  200  m,  the 
depth  of  the  tip  of  a  sea  mount  located  ~12.5  km  to  the  southwest  of  the  mooring  site 
(Figure  3-3).  Notice  also  that  there  are  some  interannual  differences  at  this  site.  The 
flow  is  stronger  below  100  m  during  2001  (e.g.  3.1±0.5  cm/s  at  250  m  in  B2)  than 
during  2002  (2.0±0.2  cm/s  at  250  m  in  C3).  Finally,  the  B3  mooring,  which  shows 
a  strong  mean  with  weak  polarization  in  the  top  100  m  becomes  strongly  polarized 
at  depth  but  with  relatively  weak  (<1  cm/s)  means.  Notice  that  the  flow  becomes 
aligned  in  the  general  direction  of  the  small  side  extension  of  Marguerite  Trough,  and 
not  the  main  axis  of  the  trough  as  with  the  other  moorings  located  in  it. 

3.2.3  Consistency  with  Previous  Studies 

As  noted  in  the  introductory  chapter,  the  SO  GLOBEC  moored  array  provides  the 
first  set  of  direct  current  observations  spanning  several  seasons  on  the  wAP  shelf. 
However,  a  handful  of  previously  published  observations  are  available  for  comparison 
with  the  analysis  presented  here.  In  two  separate  studies  using  CTD  observations 
to  calculate  geostrophic  shear  (Smith  et  al.,  1999;  Klinck  et  al.,  2004),  it  was  sug¬ 
gested  that  the  nearshore  circulation  along  Adelaide  Island,  which  was  shown  to  be 
consistently  towards  the  southwest  during  the  ice- free  season,  might  be  part  of  a  shelf- 
scale  gyre  or  a  coastal  current.  This  was  shown  to  be  consistent  with  a  study  using 
drifter  observations  of  the  near-surface  lagrangian  currents  collected  during  the  SO 
GLOBEC  program  (Beardsley  et  al.,  2004).  The  coastal  current,  which  dominates 
the  flow  at  Al,  and  which  is  largely  baroclinic,  is  explored  further  in  Chapter  5. 

As  for  the  shelf-scale  cyclonic  gyre,  the  SO  GLOBEC  broad-scale  dynamic  height 
calculations  relative  to  a  zero-velocity  reference  velocity  at  400  m  showed  it  to  be 
centered  off  Adelaide  Island,  roughly  at  the  mid-shelf  point  (67°  S),  and  it  is  a  feature 
which  is  relatively  persistent  from  the  fall  to  the  winter  (Klinck  et  al.,  2004).  This  gyre 
is  consistent  with  the  cyclonic,  largely  barotropic  circulation  described  here  from  the 

4Note,  however,  that  the  weak  polarization  implies  calling  an  axis  ’major’  and  the  other  ’minor’ 
axis  somewhat  ambiguous. 
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moored  array  data  taken  around  the  trough,  although  the  predicted  velocity  from  the 
dynamic  height  calculations  underestimate  the  total  velocity.  Even  near  the  coast  of 
Adelaide  Island,  where  the  flow  is  largely  baroclinic,  the  dynamic  height  calculations 
gave  maximum  velocities  of  the  order  of  3  cm/s  during  the  fall  2001  season,  although 
the  measured  surface  velocities  at  Al  were  closer  to  40  cm/s.  This  appears  to  be  a 
direct  consequence  of  the  sampling  grid  being  too  coarse  in  the  cross-shelf  direction 
(20  km)  to  resolve  the  small  dynamical  scales  (the  internal  radius  of  deformation  is 
roughly  4  km). 

3.2.4  Blocking  Effects  in  the  Mean  Flow 

The  above  description  of  the  flow  highlights  the  obvious  influence  of  the  local  topog¬ 
raphy  on  the  alignment  of  both  the  mean  and  axes  of  variability  around  Marguerite 
Trough,  where  most  of  the  moorings  were  located.  However,  much  of  the  area  of  the 
shelf  is  occupied  by  small-scale  bathymetric  features  which  are  also  likely  to  have  a 
profound  effect  on  the  circulation. 

Consider,  for  example,  the  vertical  structure  of  the  mean  flow  at  Al.  The  flow  in 
the  top  200  m  is  dominated  by  the  APCC,  a  strongly  baroclinic  current.  Downstream 
(to  the  southwest)  from  the  mooring,  there  is  a  ridge  extending  from  the  coast  which 
the  flow  below  300  m  (the  approximate  sill  depth)  must  negotiate.  If  the  height  of 
an  obstacle  in  the  path  of  a  stratified,  steady  flow  is  relatively  small,  the  flow  will 
generate  internal  waves.  If  the  obstacle  is  large  enough,  the  fluid  will  instead  be 
blocked  upstream.  This  critical  scale  height  is  given  by  (Gill,  1982;  Baines,  1995) 

=  TV’  (3-2) 

which  can  be  compared  to  the  height  of  the  obstacle  hs.  At  Al,  hs  =200  m,  and  a 
reasonable  estimate  for  the  velocity  above  the  sill  and  the  stratification  are  U  =  5  cm/s 
and  TV  =  1.3T0-3s_1  respectively,  which  results  in  hm  ~40  m.  Clearly,  h8  ^  hm  and 
the  flow  is  expected  to  be  blocked  below  the  sill.  This  is  evident  in  Figure  3-2,  showing 
the  mean  velocities  for  the  moored  array.  While  the  flow  above  the  sill  appears  aligned 
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with  the  coastline,  the  flow  at  309  m  (roughly  at  the  sill  depth)  shows  some  offshore 
veering,  and  the  deepest  instrument  shows  a  mean  perpendicular  to  the  coastline, 
and  parallel  to  the  orientation  of  the  ridge  located  downstream  of  the  mooring  site. 
Notice  that  relatively  small  magnitude  of  the  flow  at  and  below  the  sill  depth  is  not 
surprising  given  the  baroclinic  nature  of  the  APCC.  The  orientation  of  the  mean, 
however,  appears  consistent  with  the  idea  of  blocking.  As  we  will  see  later  in  this 
chapter,  the  tidal  velocities  have  similar  velocity  amplitudes  and  the  ridge  also  blocks 
their  motion  at  depth 
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Table  3.1:  Mean  north  and  east  velocity  (u,v),  standard  deviation  (SD),  mean  rotated  velocity  (ur  and  vr)  and 
standard  deviation  (SD,),  and  major  axis,  minor  axis  and  inclination  of  the  principal  axis  of  variability  of  the 
SO  GLOBEC1  moored  currents  records,  all  in  cm/s.  All  records  were  low-pass  filtered  (33-h  half-amplitude)  before 
computing  the  statistics.  The  uncertainties  of  the  mean  are  the  95%  confidence  intervals.  The  angle  of  rotation 
(counterclockwise  from  east)  is  indicated  next  to  the  mooring  name. 
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Short  records  (<6  months) 


Figure  3-2:  Mean  and  principal  axes  of  the  low-pass  filtered  (33-h  half-amplitude) 
velocity  records  at  50,  100  and  250  m  (top  panels)  and  400  m  and  800  m  (bottom 
panels).  The  red  crosses  indicate  the  2001-2002  deployment  (the  A-  and  B-lines)  and 
the  green  the  2002-2003  deployment  (C-line).  At  Al,  the  instruments  labeled  250 
and  400  m  were  at  309  and  459  m,  respectively. 
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Figure  3-3:  Bathymetry  and  orientation  moored  array  current  records.  The  center 
panel  shows  a  map  of  the  study  area  around  Marguerite  Trough,  with  the  500-,  1000-, 
2000-  and  3000-m  isobaths  plotted.  The  panels  on  either  side  show  detailed  maps  of 
the  bathymetry  (with  contours  every  100  m)  in  boxes  around  each  mooring,  indicated 
by  the  dashed  lines  in  the  center  panel.  At  each  site,  the  positive  x-axis  is  aligned 
with  the  local  isobaths  and  pointing  towards  the  shelf-break,  with  the  exception  of 
Al,  where  the  y-axis  is  aligned  with  the  coastline,  and  B3.  where  the  direction  of  the 
local  isobaths  is  ambiguous  (see  text). 
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3.3  Tidal  Variability 


In  this  section,  an  overview  of  the  tidal  current  variability  on  the  wAP  shelf  is  pre¬ 
sented.  The  analysis  was  carried  out  using  a  subset  of  the  SO  GLOBEC  array  velocity 
and  pressure  records  using  the  T.TIDE  toolbox  for  MATLAB  (Pawlowicz  et  al.,  2002). 
The  toolbox  is  used  to  perform  harmonic  analysis  of  real  (i.e.  pressure)  or  complex 
(i.e.  velocity)  time  series.  The  product  of  the  analysis  of  a  velocity  time  series  is  a  pair 
of  complex  amplitudes  (a*.a_*)  for  each  tidal  constituent  k,  and  they  are  converted 
to  the  standard  tidal  parameters  reported  here: 


Majfc  — 

|a*|  +  Ia-fc| 

(3.3) 

Min*  - 

K-l  - 1«-*| 

(3.4) 

Inc*  - 

an3(at)  +  an9(a_,)modl80 

(3.5) 

Phase*  = 

Vk  -  ang(cik)  +  Inc*, 

(3.6) 

which  describe  an  ellipse  with  semi-major  and  semi-minor  axes  Maj  and  Min  and 
inclination  Inc ,  which  is  measured  counter-clockwise  from  east.  If  Min  is  positive 
(negative),  the  ellipse  is  traced  counter-clockwise  (clockwise).  The  phase  is  given 
as  "Greenwich  phase”,  i.e.  the  phase  of  the  response  referenced  to  the  equilibrium 
response  at  0°  Longitude  (e*  is  the  equilibrium  phase).  In  the  case  of  a  pressure 
record,  the  values  reported  are  the  amplitude  of  the  response  and  the  Greenwich 


In  order  to  separate  neighbouring  constituents  separated  in  frequency  by  A  f.  the 
time  series  must  be  of  length 

|A/|r  >  R  (3.7) 

where  T  is  the  record  length  and  R  (called  the  Raleigh  criterion)  is  typically  equal 
to  unity  (Emery  and  Thompson,  1998).  For  example,  separation  of  the  AT  and  Pi 
diurnal  components  with  R  =  1  requires  a  record  length  of 


T  > 


1 


1 


/(AT)  -  /(Pi)  0.0418  -  0.0416  cpd 


=  182  days. 


(3.8) 
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For  the  analysis  presented  here,  only  long  records  of  between  302  and  364  days  were 
selected,  which  allows  for  the  separation  of  the  most  energetic  components  observed. 
Notice  also  that  at  67°  S,  the  inertial  frequency  is  0.0226  cpd,  so  that  motions  at 
this  frequency  can  be  separated  from  motions  at  the  M2  tidal  frequency  with  a  time 
series  1 1-days  long  or  so,  with  the  caveat  that  motions  forced  for  example  by  the  wind 
usually  occupy  a  broad  band  near  /. 

The  instruments  recorded  velocity  and  pressure  at  nominal  depths  of  100,  250, 
400  and  800  m.  Significant  components  are  defined  as  having  a  signal-to-noise  (SNR) 
ratio  of  2  or  more.  The  SNR  is  calculated  as  the  squared  ratio  of  the  amplitude 
to  the  amplitude  error  (Pawlowicz  et  al.,  2002).  The  results  of  the  pressure  record 
analysis  are  reported  in  centimeters  of  water,  which  were  converted  from  pressure 
units  assuming  hydrostatic  balance  and  a  density  of  1027  kg/m3. 


Figure  3-4:  Spectrum  of  the  pressure  record  at  47  m  at  A2  showing  the  main  tidal 
constituents  for  the  SO  GLOBEC  study  region. 
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Overall,  the  analysis  reveals  that  tidal  currents  are  typically  less  than  2  cm/s. 
The  four  largest  tidal  components  are  the  semi-diurnal  M2  and  S2  and  the  diurnal 
0\  and  K\  (Figure  3-4).  These  four  components  explain  between  5  (at  100  m  at  A2) 
and  55%  (at  depth  at  B3)  of  the  variance  (18%  on  average)  of  the  velocity  records. 

In  the  following  sections,  the  semi-diurnal  and  diurnal  components  are  considered 
separately. 

3.3.1  Semi-diurnal  Components 

Velocity 

The  semi-diurnal  components  show  a  few  important  similarities.  Tidal  currents  are 
of  the  order  of  1  to  2  cm/s  everywhere.  On  the  mid-  and  outer-shelf  at  A2  and  A3, 
the  currents  (Table  3.3  and  Figures  3-5  and  3-6)  are  barotropic,  while  the  rest  of 
the  locations  show  significant  vertical  structure.  Overall,  the  S2  component  shows 
stronger  polarization  at  all  depths. 

Down  to  250  m,  the  M2  tidal  ellipses  are  roughly  aligned  with  the  along-shelf 
direction  north  of  Marguerite  Bay  (the  clockwise  rotation  of  the  ellipses  with  offshore 
distance  is  not  statistically  significant),  while  at  the  B-  and  C-mooring  lines  they  are 
roughly  oriented  east-west,  in  the  direction  of  the  mouth  of  the  bay.  For  the  same 
upper  layer,  the  S2  component,  which  in  the  nearshore  mooring  north  of  the  bay  (Al) 
is  also  strongly  aligned  with  the  coast,  turns  anti-clockwise  with  distance  from  the 
coast  and  is  oriented  perpendicular  to  the  coast  at  A3.  Across  the  mouth  of  the  bay, 
at  B3,  Cl  and  C3,  the  S2  ellipse  is  also  roughly  aligned  with  the  M2  components, 
while  at  B2/C3,  the  S2  ellipses  are  oriented  north-south.  Given  that  B2/C3  and  B3 
are  only  25  km  apart,  this  suggests  that  there  is  significant  spatial  variability  in  the 
orientation  of  the  semi-diurnal  currents  around  Marguerite  Bay. 

At  Al  and  at  the  B-  and  C-mooring  lines  located  across  the  mouth  of  Marguerite 
Bay,  the  semi-diurnal  currents  have  a  complicated  vertical  structure.  At  Al,  the  tidal 
currents  at  100  and  250  m  are  roughly  aligned  with  the  coast  of  Adelaide  Island, 
while  the  deepest  instrument,  located  below  the  sill  of  the  deep  depression  where  the 
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mooring  was  deployed  (see  Section  3.2.1  and  Chapter  5),  shows  strong  veering  of  the 
currents  so  that  they  are  aligned  perpendicular  to  the  coast  and  roughly  parallel  to 
the  orientation  of  the  nearby  sill  (Figure  3-2).  There  is  also  a  roughly  100°  phase  dif¬ 
ference  between  the  currents  above  and  below  the  sill  for  both  components,  although 
the  magnitude  of  the  currents  does  not  change  significantly. 

The  strong  influence  of  local  bathymetry  is  also  evident  elsewhere  on  the  shelf.  At 
C2  and  B2/C3,  the  currents  at  and  below  400  m  appear  increasingly  aligned  with  the 
local  bathymetry  with  increasing  depth,  and  also  show  phase  differences  of  60-150° 
above  and  below  400  m  (with  the  exception  of  the  M2  at  B3).  The  magnitude  of 
the  velocity  also  is  reduced  with  depth  at  C2,  while  at  B3  the  velocity  increases  with 
depth  for  both  components. 

The  vertical  structure  of  the  semi-diurnal  components  is  perhaps  most  interesting 
at  B2/C3.  The  estimates  are  consistent  from  the  first  (B2)  to  the  second  (C3)  year 
of  the  deployment,  and  show  the  velocity  decreasing  and  turning  anti-clockwise  from 
100  to  400  m  at  both  moorings.  However,  the  deepest  instrument  at  C3,  located 
at  800  m,  reveals  that  at  depth  the  velocity  increases  again,  so  that  the  velocity 
magnitude  shows  a  minimum  at  400  m. 

Sea-level  Amplitude 

The  tidal  sea-level  amplitude  on  the  shelf  (Table  3.5  and  Figure  3-11,  top  panels) 
is  maximum  along  the  shelf-break  for  M2  with  a  small  (1  cm)  drop  in  the  «110  km 
covered  by  the  A-mooring  line.  The  phase  estimations,  which  show  larger  errors,  are 
not  significantly  different  across  the  shelf.  Across  the  mouth  of  the  bay  going  from 
north  to  south,  from  Cl  to  B3  (Figure  3-12),  the  amplitude  also  drops  by  less  than 
1  cm,  with  a  relatively  constant  phase. 

The  S2  tide  show  an  opposite  structure  across  the  shelf,  with  a  maximum  near  the 
coast  and  a  1-cm  drop  toward  the  outer  mooring.  The  phase  estimates  also  show  an 
8°  to  10°  increase  in  the  phase  from  the  coast  to  A3.  Across  the  mouth  of  the  bay, 
the  structure  is  similar  to  the  M2  tide,  with  a  1-cm  or  so  drop  in  sea-level  from  Cl  to 
B3.  The  variability  in  the  phase  estimations,  even  at  the  B2/C3  pair,  is  too  large  to 
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detect  any  spatial  patterns  of  variability,  although  the  entire  range  is  only  5  °  (Table 
3.5). 

3.3.2  Diurnal  Components 

Velocity 

The  tidal  ellipses  and  phases  of  the  two  most  important  diurnal  components,  0\ 
and  Ki,  are  shown  in  Figures  3-7  and  3-8  and  Table  3.4.  The  structure  shares 
some  characteristics  with  the  semi-diurnal  components.  At  the  mid-  and  outer-shelf 
moorings  of  the  A-line  (A2  and  A3)  the  diurnal  components  are  barotropic,  while  at 
Al  the  motion  below  the  sill  shows  strong  veering  and  a  phase  difference  of  130°  or 
so  with  respect  to  the  overlying  water. 

The  0\  currents  have  a  similar  spatial  structure  to  M2  north  of  Marguerite  Bay. 
The  tidal  ellipses  are  strongly  polarized  near  the  coast  and  aligned  with  the  coastline. 
The  ellipses  show  clockwise  rotation  with  distance  from  the  coast  (the  uncertainties 
are  in  this  case  much  smaller),  with  the  ellipse  becoming  less  polarized  in  the  outer 
shelf.  Across  the  mouth  of  Marguerite  Bay,  the  tidal  ellipses  have  a  complex  horizontal 
structure,  with  the  currents  near  the  northern  edge  of  the  entrance  to  the  bay  oriented 
along  the  axis  of  the  mouth,  while  farther  to  the  south-west  at  B2/C3,  C2  and  B3  the 
ellipses  show  markedly  distinct  structures,  with  relatively  small  (1  cm/s)  velocities 
with  weak  polarization  at  B2/C3  and  strong  (3-4  cm/s)  velocities  at  B3,  oriented  in 
the  east-west  direction.  The  above  description  applies  also  for  the  K \  currents,  with 
a  few  exceptions:  the  tidal  ellipses  across  the  A-mooring  line  are  more  similar  to  S2  in 
that  the  ellipses  show  a  strong  90°  clockwise  rotation  from  the  nearshore  Al  mooring 
to  the  outer- shelf  A3  mooring.  Also,  the  B2/C3  moorings  show  stronger  polarization 
of  the  diurnal  ellipses. 

The  vertical  structure  of  the  diurnal  constituents  in  the  B-  and  C-mooring  lines 
also  shows  some  structure  with  depth,  but  is  generally  not  as  significant  as  that  shown 
by  the  semi-diurnal  components.  At  B3,  where  both  the  0\  and  K\  are  the  largest 
on  the  shelf,  diurnal  currents  increase  significantly  from  250  to  400  m,  although  there 
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is  no  significant  veering  and  only  a  weak  phase  difference  with  depth.  At  B2/C3, 
a  vertical  structure  resembling  the  semi-diurnal  components  is  also  evident,  with  a 
velocity  minimum  at  400  m.  Finally,  the  tidal  ellipses  at  C2,  which  are  roughly 
aligned  with  the  bathymetry  above  the  rim  of  Marguerite  Trough,  becoming  more 
strongly  polarized  with  depth  but  show  little  veering  and  phase  changes. 

Sea-level  Amplitude 

The  diurnal  sea-level  components  show  a  rather  simple  structure  along  the  A-line, 
with  both  components  showing  a  constant  amplitude  across  the  shelf.  The  phase 
estimations  have  too  much  spread  to  determine  trends,  but  in  both  cases  the  phase 
changes  are  less  than  5  °. 

Across  the  mouth  of  Marguerite  Bay,  the  diurnal  components  resemble  the  semi¬ 
diurnal  components,  with  the  highest  amplitude  on  the  northern  edge  of  the  mouth 
near  Adelaide  Island  and  the  0\  (K i)  tide  showing  a  1.5-cm  (1-cm)  drop  in  amplitude 
towards  the  southernmost  mooring.  The  phase  is  nearly  constant  along  this  line. 

3.3.3  Comparison  with  the  AntPen04.01  Tidal  Model 

A  regional  context  for  the  tidal  analysis  results  provided  in  the  previous  section  can 
be  obtained  from  comparisons  with  the  AntPen04.01  regional  tidal  model,  developed 
by  Laurie  Padman  and  colleagues  at  Earth  and  Space  Research  (http:/ /www. esr.org). 
The  following  sections  provide  a  brief  overview  of  the  model  and  the  first  comparison 
between  the  model  and  moored  currents  observations  for  this  region. 

A  brief  overview  of  the  model 

The  AntPen04.01  is  a  high-resolution  (1/30°  Longitude  x  1/60°  Latitude,  roughly 
2  km),  barotropic,  and  linear  forward  tidal  model  whose  domain  covers  the  entire 
Antarctic  Peninsula.  The  model  includes  the  best  bathymetry  available  to  date,  in¬ 
cluding  the  high-resolution  bathymetry  for  the  SO  GLOBEC  study  region  (Bolmer 
et  al.,  2004).  The  model  is  forced  at  its  lateral  boundaries  by  CATS02.01,  a  circum- 
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Antarctic,  coarser-resolution  forward  model  (Padman  et  al.,  2002)  and  by  astronom¬ 
ical  forcing,  and  is  based  on  the  linearized  shallow  water  equations. 


The  main  tidal  constituents  in  AntPen04.01 

The  model  sea-level  amplitude  and  phase  of  the  main  semi-diurnal  and  diurnal  com¬ 
ponents  are  shown  in  Figures  3-9  and  3-10.  For  all  of  the  components,  the  tidal 
amplitude  tends  to  be  smaller  and  with  larger  horizontal  scales  on  the  western  shelf 
of  the  peninsula  than  on  the  eastern  shelf.  The  diurnal  components  are  largest  all 
along  the  coast  except  near  the  tip  of  the  peninsula  where  M2  becomes  slightly  larger. 
In  and  around  Marguerite  Bay,  however,  the  tidal  phase  has  a  range  of  only  a  few 
degrees  for  each  component. 

Of  the  four  components  considered  here,  three  (S2,  0\  and  K\)  have  a  phase 
structure  consistent  with  anti-clockwise  propagation  around  the  tip  of  the  peninsula 
and  southward  along  the  coast.  The  diurnal  components  have  phase  fluctuations  with 
long  horizontal  scales  with  changes  of  the  order  of  10°  along  the  400  km-long  study 
region.  The  S2  tide  has  a  roughly  30°  phase  change  within  the  same  region. 

The  M2  tide  propagates  across  the  shelf,  turning  slightly  north-eastward  north  of 
Alexander  Island  (the  southern  boundary  of  Marguerite  Bay)  and  south-eastwards 
south  of  the  bay,  with  the  phase  increasing  by  roughly  10°  from  the  shelf-break  to 
the  coast  in  the  SO  GLOBEC  study  area. 


Comparison  with  the  observations 


To  compare  the  model  output  to  the  moored  array  observations,  the  current  and 
pressure  observations  were  vertically  averaged  prior  to  performing  the  tidal  analysis. 
The  vertically-averaged  velocity  and  pressure  records  are  defined  as: 


u  = 


(3.9) 
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and 


P'  =  JjfHP'dz »  (3-10) 

where  u  is  the  velocity,  p'  is  the  pressure  anomaly  (the  pressure  record  with  the 
time  mean  removed)  and  H  is  the  bottom  depth.  In  order  to  calculate  the  average 
quantities,  the  velocity  and  pressure  anomaly  are  assumed  constant  from  the  shallower 
(deepest)  instrument  in  the  water  column  to  the  surface  (bottom).  Several  of  the 
ADCP  records  in  the  top  50  m  or  so  have  long  periods  of  time  (in  particular  during 
the  winter)  with  no  data,  and  B2  and  B3  had  a  few  instruments  change  depth  due 
to  loss  of  buoyancy  of  their  top  instrument  (B2)  and  flotation  sphere  (B3)  (Moffat 
et  ah,  2005).  The  vertical  averages  were  therefore  computed  using  a  subset  of  the 
data  that  provided  measurements  throughout  the  entire  deployment  and,  in  the  case 
of  the  ADCP  records,  had  at  least  90%  of  data  return  during  that  period.  Table  3.2 
presents  a  list  of  the  instrument  records  used  to  calculate  the  vertical  averages. 

Figures  3-5  to  3-8  show  the  tidal  ellipses  from  the  AntPen04.01  model  together 
with  the  tidal  analysis,  and  Table  3.6  presents  the  tidal  ellipse  parameters  for  each 
mooring  site  from  both  the  model  output  and  the  vertically- averaged  records.  Let 
us  consider  first  the  tidal  currents  along  the  A-line.  For  the  four  components  consid¬ 
ered,  the  model  predicts  the  semi-diurnal  components  well  for  all  the  sites,  although 
it  overestimates  the  magnitude  of  M 2  at  Al.  For  0\,  the  model  tends  to  underesti¬ 
mate  the  magnitude  of  the  velocity  at  the  three  moorings  by  25%  to  50%.  The  0\ 
inclination  is  predicted  correctly  at  Al,  but  not  at  A2  and  A3,  and  the  phase  is  within 
error  at  Al  and  A2  and  is  too  small  at  A3  by  7°  (some  30  minutes).  The  model  also 
slightly  underestimates  the  K\  currents  at  Al,  and  the  inclination  and  phase  is  off  by 
70  and  100°  respectively.  At  A2  and  A3,  the  model  correctly  predicts  the  magnitude 
and  phase  of  the  K\  currents,  while  the  predicted  inclination  is  slightly  outside  the 
estimated  errors. 

Across  the  mouth  of  Marguerite  Bay,  the  model  shows  mixed  results  in  predicting 
the  tidal  currents.  The  observed  M2  tide  current  magnitude  is  overestimated  by  the 


55 


Table  3.2:  Depth  of  the  current  velocity  and  pressure  records  used  to  compute  the 
vertical  averages.  ADCP  records  are  given  as  a  range,  with  the  bin  separation  in 
parenthesis  (97-165  (2)  means  from  97  to  165  m  every  2  m),  and  the  instruments 
are  separated  by  commas.  Also  given  are  the  bottom  depth  at  each  mooring  and 
the  length  of  the  record  used. 


Mooring  Velocity 

Record  Depth 

M 

Pressure 

Record 

Depth  [m] 

Bottom  Depth 

M 

Length  [days] 

A1 

97-165  (2),  309, 

309,  502 

509 

323 

459 

A2 

107,  247,  397 

47,  247,  554 

561 

319 

A3 

103,  253,  403 

103,  253,  480 

480 

319 

B2 

91.7-289.7  (6), 

804 

811 

321 

400,  811 

B3 

247,  400 

247,  400 

322 

Cl 

26-88  (2),  232 

232 

859“ 

368 

C2 

260,  410,  810 

260,  852 

859 

367 

C3 

13-109  (2),  252, 

252 

806 

364 

395,  762 

°  The  integration  was  conducted  assuming  H  =  250  m 


model  at  B3  and  underestimated  at  Cl  (in  both  cases  by  «25%)  but  is  predicted 
correctly  everywhere  else.  The  M2  phase  is  within  the  error  at  B3  and  Cl  but  it  is 
not  predicted  correctly  at  B2,  C2  and  B2/C3.  The  inclination  from  the  model  for 
the  M2  tide  is  outside  the  confidence  intervals  at  B3  and  C2.  The  amplitude  of  the 
S2  tide  is  in  general  predicted  correctly,  with  small  (10%-20%)  overestimations  at  B2 
and  C2  and  a  bigger  underestimation  (25%)  at  C2.  The  inclination  and  phase  tend 
to  be  off  by  10  to  30°  (for  phase,  20  minutes  to  one  hour),  with  as  much  as  150°  of 
error  in  the  phase  at  B3. 

The  diurnal  components  are  also  predicted  reasonably  well  by  the  model.  Dif- 
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ferences  between  the  model  and  the  data  larger  than  the  confidence  intervals  are 
generally  within  15%  or  so  of  the  estimations  from  the  observations.  The  only  sig¬ 
nificant  exception  is  at  B3,  where  the  model  does  the  worst  of  all  the  mooring  sites 
for  the  diurnal  components,  underestimating  the  amplitude  of  the  velocity  by  around 
67%  and  predicting  the  phase  incorrectly. 

Comparisons  of  the  estimated  tidal  sea-level  amplitudes  and  phases  with  the  model 
predictions  are  shown  together  in  Figures  3-11  and  3-12.  At  the  A-line,  the  model 
underestimates  the  amplitude  of  the  M2,  Oi  and  K\  constituents  by  ~5%  (0.5  to 
1.5  cm),  while  overestimating  the  52  amplitude  by  1%  (0.2  cm)  or  so.  The  phase 
predictions  are  either  within  error  or  within  a  few  degrees  of  the  estimated  values. 
This  translates  into  an  error  of  2  to  4  minutes,  which  is  in  general  larger  than  the 
clock  drifts  found  upon  instrument  recovery,  although  a  few  instruments  were  found 
to  have  clock  errors  of  a  few  minutes  over  the  1-year  deployment  (Moffat  et  al.,  2005). 

Across  the  mouth  of  Marguerite  Bay,  the  model  amplitudes  also  predict  the  ob¬ 
servations  well.  The  model  consistently  underestimates  the  M2  and  K\  amplitude  by 
1-4%  (0.2  to  1  cm)  while  errors  in  the  other  two  components,  where  they  are  signifi¬ 
cant,  are  of  the  same  magnitude.  The  phase  estimates  are  significantly  different  than 
the  model  in  a  few  locations,  with  the  K\  phase  predictions  being  consistently  too 
small  across  the  bay  by  1  °  or  so  (4  minutes),  but  overall  the  errors  are  smaller  than 
a  few  degrees. 

3.3.4  Discussion  and  Summary  of  the  Tidal  Analysis 

The  analysis  presented  in  the  above  sections  shows  the  tidal  velocities  are  small  in  the 
study  region,  with  velocities  as  large  as  4  cm/s  but  in  general  less  than  2  cm/s.  The 
tidal  currents  are  fairly  barotropic  on  the  mid-  and  outer-shelf  off  Adelaide  Island,  but 
the  analyses  in  the  above  sections  show  that  the  tidal  currents  can  have  significant 
vertical  structure. 

At  C2,  for  example,  the  semi-diurnal  components  are  oriented  at  a  significant 
angle  with  respect  to  the  axis  of  Marguerite  Trough,  and  shows  strong  veering  and 
phase  shifts  below  the  rim  of  the  trough.  The  diurnal  components,  however,  are 
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roughly  aligned  with  the  axis  of  the  trough  and  one  of  the  effects  of  the  bathymetry 
is  to  increase  the  polarization  of  these  currents  with  depth.  At  Al,  where  all  the 
tidal  components  are  aligned  with  the  coastline  above  the  sill,  the  deepest  instrument 
always  show  strong  veering  and  phase  shifts  to  align  the  deep  currents  along  the 
direction  of  the  nearby  sill,  oriented  roughly  perpendicular  to  the  coastline. 

Other  properties  of  the  tidal  velocity  field  are  more  complex.  At  B2/C3,  for 
example,  the  tidal  velocities  show  a  minimum  at  400  m,  and  although  the  first- 
year  deployment  lacked  a  deep  instrument,  the  vertical  structure  down  to  400  m  is 
consistent  in  both  years,  suggesting  this  mid-depth  minimum  is  not  due  to  instrument 
malfunction. 

Some  of  the  ideas  about  flow  blocking  discussed  before  (Section  3.2.4)  also  apply 
to  oscillatory  flows  (Legg,  2004;  Garrett  and  Kunze,  2007)  and  therefore  blocking  of 
the  flow  is  also  expected  at  Al.  To  understand  the  vertical  structure  at  some  of  the 
mooring  sites,  simple  dynamical  scales  based  on  models  of  flow  past  isolated  obstacles 
are  likely  not  appropriate  and  more  complicated  models  of  tidal,  stratified  flows  will 
be  required. 

As  the  barotropic  tide  interacts  with  topography,  internal  gravity  waves  known  as 
internal  tides  might  be  forced.  The  internal  tides  extract  energy  from  the  barotropic 
tide  and  this  energy  could  be  propagated  away  or  be  dissipated  on  the  shelf.  Although 
a  detailed  study  of  internal  tides  on  the  wAP  shelf  is  beyond  the  scope  of  the  basic 
description  presented  here,  a  brief  exploration  of  the  relevant  scales  of  internal  tide 
dynamics  might  shed  some  light  on  their  relevance  in  the  study  region.  Radiating 
waves  can  be  expected  for  motions  with  frequencies  f  <  u  <  N.  For  the  wAP 
region,  both  the  M2  and  S2  semi-diurnal  components,  which  would  act  as  the  forcing, 
are  within  this  range  and  therefore  radiation  of  internal  waves  at  these  frequencies 
can  be  expected.  In  fact,  the  diurnal  components  tend  to  be  more  barotropic  than 
the  semi-diurnal  components,  which  suggest  that  the  transfer  of  energy  from  the 
barotropic  to  the  baroclinic  tide  might  be  more  relevant  for  the  semi-diurnal  tides.  If 
one  assumes  a  simple,  two-dimensional  stratified  ocean  with  sinusoidal  topography, 
where  there  is  a  tidal  flow  with  scale  uQ,  an  important  parameter  of  the  flow  is  u()/uj. 
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the  tidal  excursion.  For  the  relatively  small  tidal  velocities  observed  on  the  wAP 
shelf  (1-2  cm/s),  the  tidal  excursion  is  only  of  the  order  of  a  kilometer  or  so  for  the 
semi-diurnal  components.  Moreover,  the  scale  of  the  topography  k -1  can  be  much 
larger  than  the  tidal  excursion,  so  that  the  advective  terms  in  the  governing  equations 
need  to  be  retained  and  lee  waves  might  be  generated  (Garrett  and  Kunze,  2007). 
Another  parameter  of  importance  is  the  topographic  steepness,  defined  as 


7  = 


kho 

c 


(3.11) 


where 

/w2_  ra\l/2 

■  (312> 

where  ho  is  a  height  scale  of  the  topography,  and  c  is  the  slope  of  the  particle  motions 
and  of  the  internal  wave  rays  with  frequency  u>  propagating  in  a  fluid  of  stratification 
N  (Baines,  1974).  When  the  steepness  parameter  7  is  near  1,  the  topographic  slope  is 
said  to  be  critical,  and  regions  of  near-critical  slope  have  been  shown  to  have  enhanced 
baroclinic  tidal  motion  (Leaman,  1980;  Holloway,  1985;  Garrett  and  Kunze,  2007). 

On  the  wAP  shelf,  c  is  roughly  0.03  for  the  M2  tide  (with  N  =  1.3 10-3  s_1).  With 
the  complicated  bathymetry  typical  of  the  region  (see  Chapter  2),  critical  slopes  are 
ubiquitous  on  the  shelf.  As  an  example,  the  side  walls  of  Marguerite  Trough  (Figure 
2-5)  have  scales  h0  of  the  order  of  300  to  1000  m  over  lateral  scales  AT1  of  only  10 
to  20  km,  which  result  in  values  for  7  of  3  to  22,  clearly  supercritical.  This  simple 
analysis  suggests  that  internal  tides  forced  by  the  semi-diurnal  barotropic  tide  might 
be  important  on  the  wAP  region,  and  that  even  though  the  tidal  currents  are  small, 
the  strongly  sloping  bathymetry  might  provide  an  effective  mechanism  for  the  transfer 
of  energy  from  the  barotropic  tide  to  internal  waves,  which  might  in  turn  enhance 
mixing  on  the  shelf. 

As  for  the  diurnal  tidal  currents,  which  are  of  lower  frequency  than  the  local  Cori¬ 
olis  frequency  /,  the  evidence  indicates  they  tend  to  be  more  barotropic  everywhere, 
which  is  consistent  with  the  absence  of  internal  tide  generation  based  on  simple  mod- 
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els  (Garrett  and  Kunze,  2007).  However,  some  of  the  vertical  structure  which  is 
observed  at  some  sites  might  be  due  to  the  blocking  effects  by  bathymetry  discussed 
above. 

Finally,  a  comparison  of  the  tidal  analysis  results  with  the  AntPen04.01  model 
output  revealed  that  the  model  predicts  reasonably  well  the  amplitude  and  phase  of 
the  four  major  constituents  considered  here  to  within  a  few  percent.  However,  the 
observed  tidal  currents  have  significant  vertical  structure  at  several  of  the  mooring 
sites  which  appears  to  be  influenced  by  the  complicated  bathymetry  on  the  shelf.  The 
model  predicts  reasonably  well  the  tidal  velocity  at  the  A-mooring  line,  in  particular 
the  offshore  moorings  where  the  tides  are  barotropic,  but  it  has  mixed  results  across 
the  mouth  of  Marguerite  Bay,  where  the  tidal  currents  have  a  more  complex  vertical 
structure. 


60 


Table  3.3:  Major  axis,  minor  axis,  inclination,  phase  and  length  of  the  records  of  the  largest  semi-diurnal  com¬ 
ponents,  M2  and  S2 •  The  inclination  is  measured  counter-clockwise  from  east,  and  the  phase  is  given  in  degrees 
Greenwich. 
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Table  3.4:  Major  axis,  minor  axis,  inclination,  phase  and  length  of  the  records  of  the  largest  diurnal  components, 
Oi  and  K\. 
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Figure  3-5:  Tidal  ellipses  for  the  M>  component  as  a  function  of  depth  from  the  SO  GLOBEC  moored  array.  Ihe  left  panel 
shows  the  tidal  ellipses  for  each  mooring  as  a  function  of  depth  (see  also  the  maps  in  the  right  panel),  as  well  as  the  Ant  Pen04.01 
model  prediction  (top  row). 


Figure  3-6:  Tidal  ellipses  for  the  S2  component  as  a  function  of  depth  from  the  SO  GLOBEC  moored  array.  The  left  panel 
shows  the  tidal  ellipses  for  each  mooring  as  a  function  of  depth  (see  also  the  maps  in  the  right,  panel),  as  well  as  the  Ant.Pen04.01 
model  prediction  (top  row). 
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Figure  3-7:  Tidal  ellipses  for  the  01  component  as  a  function  of  depth  from  the  SO  GLOBEC  moored  array.  The  left  panel 
shows  the  tidal  ellipses  for  each  mooring  as  a  function  of  depth  (see  also  the  maps  in  the  right,  panel). 


Figure  3-8:  Tidal  ellipses  for  the  K1  component  as  a  function  of  depth  from  the  SO  GLOBEC  moored  array.  The  left  panel 
shows  the  tidal  ellipses  for  each  mooring  as  a  function  of  depth  (see  also  the  maps  in  the  right  panel),  as  well  as  the  Ant.Pen04.01 
model  prediction  (top  row). 
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Figure  3-9:  Tidal  sea-level  amplitude  (in  cm,  top  panels)  and  Greenwich  phase  (bot¬ 
tom  panels)  of  the  largest  semi-diurnal  components  from  the  AntPen04.01  tidal  model. 
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Figure  3-10:  Tidal  sea-level  amplitude  (in  cm,  top  panels)  and  Greenwich  phase 
(bottom  panels)  of  the  largest  diurnal  components  from  the  AntPen04.01  tidal  model. 
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Table  3.5:  Estimated  tidal  sea-level  amplitude  ( Aobs )  and 
Greenwich  phase  ((fiobs),  calculated  from  the  SO  GLOBEC 
array  pressure  records.  The  uncertainties  correspond  to 
the  95%  confidence  interval  for  each  estimate  (Pawlowicz 
et  al.,  2002). 


Mooring 

Depth 

M 

Aobs 

[cm] 

(fiobs 

[°G] 

Aobs 

[cm] 

(fiobs 

[°  G] 

Al 

309 

m2 

11.5±0.3 

262±1 

s2 

15.5±0.3 

55±1 

502 

11.5±0.4 

265±2 

15.5±0.4 

59±1 

A2 

47 

12.1±0.4 

263±2 

14.9±0.3 

59±1 

247 

12.1±0.2 

266±1 

14.9±0.2 

63±1 

554 

12.1±0.2 

266±1 

14.9±0.2 

63±1 

A3 

103 

12.4±0.4 

264±2 

14.4±0.4 

65±2 

253 

12.5±0.2 

264±1 

14.5±0.2 

66±1 

480 

12.4±0.2 

262±1 

14.3±0.2 

64±1 

B2 

140 

11.6±0.3 

263±1 

17.1±0.3 

70±1 

804 

11.6±0.3 

263±1 

16.8±0.3 

70±1 

B3 

247 

11.2±0.3 

264±1 

16.9±0.3 

72±1 

440 

11.4T0.2 

263±1 

16.7±0.3 

72±1 

Cl 

232 

12.0±0.1 

263±1 

17.8±0.1 

71±0 

C2 

260 

11.9±0.1 

264±1 

17.6T0.1 

74±0 

852 

11.5±0.2 

262±1 

17.2±0.2 

71±1 

C3 

252 

11.4±0.3 

264±2 

16.6±0.2 

74±1 

Al 

309 

24.6±0.2 

71±1 

Kx 

29.5±0.2 

86±0 

502 

24.6±0.3 

73±1 

29.6±0.3 

88±1 

A2 

47 

24.7±0.5 

70±1 

29.5±0.5 

87±1 

247 

24.7±0.2 

72±1 

29.3±0.2 

89±0 

554 

24.6±0.2 

72±0 

29.4±0.2 

89±0 

A3 

103 

25.0±0.2 

71±0 

29.6±0.2 

89±0 

253 

24.9±0.2 

71±0 

29.5±0.2 

89±0 

480 

24.7±0.2 

70±0 

29.4±0.1 

88±0 

B2 

140 

23.9±0.3 

76±1 

29.5±0.2 

92±1 

804 

23.8±0.2 

75±0 

29.5±0.2 

91±0 

B3 

247 

23.0±0.3 

76±1 

29.1±0.3 

92±1 

440 

23.3±0.2 

76±1 

29.1±0.2 

92±0 

Cl 

232 

24.5±0.1 

75±0 

30.1  ±0.2 

91±0 

C2 

260 

24.2±0.2 

77±0 

29.7±0.2 

92±0 

852 

23.7±0.3 

75±1 

29.2±0.3 

91±1 

C3 

252 

23.4±0.2 

76±0 

28.8±0.2 

92±0 
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Table  3.6:  Tidal  ellipses  parameters  for  each  SO  GLOBEC  mooring  sites.  For  each  mooring,  the 
tidal  ellipses  from  the  AntPen04.01  output  (top  row)  and  calculated  from  the  vertically-averaged 
velocity  records  (bottom  row)  are  presented.  The  semi-diurnal  components  are  presented  in  the 
top  half  of  the  table,  and  the  diurnal  components  on  the  bottom  half.  The  inclination  is  measured 
anticlockwise  from  east  and  the  phase  is  given  in  degrees  Greenwich. 


Mooring 

Major 

(cm/s) 

Minor 

(cm/s) 

Incl.  (°) 

Phase 

(°G) 

Major 

(cm/s) 

Minor 

(cm/s) 

Incl.  (°) 

Phase 

(°G) 

A1 

1.5 

0.1 

m2 

57 

240 

0.8 

-0.2 

s2 

65 

213 

1.2±0.2 

0.0±0.2 

52±9 

249±10 

0.6±0.2 

-0.3±0.2 

47±27 

228±28 

A2 

1.2 

0.2 

42 

233 

0.9 

-0.3 

117 

163 

1.0±0.3 

0.2±0.3 

43±14 

248±16 

0.9±0.3 

-0.3±0.3 

118±19 

166±19 

A3 

1.5 

0.6 

29 

231 

1.2 

0.0 

135 

146 

1.3±0.3 

0.5±0.4 

37±19 

233±20 

1.1T0.3 

-0.0±0.4 

147±18 

135±17 

B2 

0.6 

0.2 

173 

11 

1.1 

-0.4 

115 

189 

0.6±0.1 

0.3±0.1 

175±18 

44±18 

0.9±0.1 

-0.2±0.1 

108±7 

196±8 

B3 

1.4 

0.6 

3 

199 

1.5 

0.1 

156 

155 

1.1±0.2 

0.3±0.2 

43±15 

187±13 

1.4±0.2 

0.0±0.2 

5±8 

308±10 

Cl 

1.0 

0.3 

156 

357 

2.0 

0.2 

143 

173 

1.3±0.2 

0.9±0.1 

157±14 

352±14 

2.7±0.1 

1.1±0.2 

153±4 

173±5 

C2 

0.6 

0.0 

158 

8 

1.2 

-0.4 

118 

191 

0.5±0.1 

0.2±0.1 

15±13 

158±14 

1.0±0.1 

-O.lzbO.l 

145±6 

171±6 

C3 

0.6 

0.2 

174 

11 

1.1 

-0.4 

115 

190 

0.5±0.1 

0.3±0.1 

165±11 

36±12 

0.8±0.1 

-0.3T0.1 

108±5 

200±5 

Al 

0.7 

0.0 

Oi 

51 

340 

0.2 

-0.1 

Aj 

107 

195 

1.3±0.1 
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Figure  3-11:  Tidal  sea-level  amplitudes  and  Greenwich  phase  for  the  A-line  moorings. 
The  solid  (hollow)  dots  are  the  amplitude  (phase)  calculated  from  the  vertically- 
averaged  pressure  anomaly  records,  and  the  error  bars  are  the  95%  confidence  limits 
for  the  estimates.  The  solid  (dotted)  line  shows  the  amplitude  (phase)  from  the 
AntPen04.01  model  output. 
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Figure  3-12:  Tidal  sea-level  amplitudes  and  Greenwich  phase  for  the  B-  and  C-line 
moorings.  The  solid  (hollow)  dots  are  the  amplitude  (phase)  calculated  from  the 
vertically-averaged  pressure  anomaly  records,  and  the  error  bars  are  the  95%  confi¬ 
dence  limits  for  the  estimates.  The  solid  (dotted)  line  shows  the  amplitude  (phase) 
from  the  AntPen04.01  model  output. 
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3.4  Near-inertial  Motions 


Examination  of  the  spectra  of  the  moored  array  current  records  over  the  wAP  shelf 
reveal  a  large,  broad  peak  at  or  near  the  inertial  frequency  (e.g.  Figure  3-1).  Some 
of  the  specific  issues  mentioned  here,  like  the  dynamics  of  wind,  ice,  and  near-inertial 
currents,  are  treated  in  Hyatt  (2006).  The  following  section  provides  a  basic  descrip¬ 
tion  of  the  horizontal,  vertical,  and  temporal  variability  of  the  near-inertial  motions 
over  the  shelf. 

3.4.1  Rotary  Spectra 

Rotary  spectra  of  the  velocity  records  show  strong  counterclockwise  energy  at  near- 
inertial  frequencies  (Figure  3-13  shows  one  typical  spectra,  for  100  m  at  A2),  which 
is  expected  for  inertial-internal  waves  in  the  southern  hemisphere  (e.g.  Kundu  and 
Allen,  1976).  The  ratio  of  counterclockwise  to  clockwise  power  spectral  density  at  the 
near-inertial  peak  for  each  record  is  between  10  and  1000  for  all  the  records  except 
for  the  deepest  instrument  at  Al,  where  there  is  no  significant  near-inertial  energy 
peak  and  the  ratio  is  only  1.7. 

The  anticlockwise  components  of  the  rotary  spectra  of  the  detided  velocity  records 
at  nominal  depths5  of  100,  250,  400  and  800  m  reveal  the  basic  structure  of  the  near- 
inertial  motions  on  the  shelf  (Figure  3-14).  Across  the  shelf  at  the  A-line,  the  energy  at 
the  spectral  peak  increases  from  the  nearshore  Al  mooring  to  the  A3  mooring,  where 
the  energy  is  the  highest  measured  at  100  m  on  the  shelf.  At  depth,  the  structure  is 
similar  at  A2  and  A3,  with  the  energy  at  the  spectral  peak  diminishing  significantly 
between  the  top  (100  m)  and  the  bottom  two  instruments  (250  and  400  m),  which 
have  inertial  peaks  of  comparable  size.  At  Al,  the  vertical  energy  distribution  is  quite 
different.  The  energy  decays  markedly  between  100  m  and  309  m,  and  at  459  m,  there 
is  no  distinguishable  inertial  peak.  As  with  the  mean  flow  and  the  tides,  the  steep 
bathymetry  and  sill  around  Al  appears  to  be  blocking  the  near-inertial  motions. 

Across  the  mouth  of  Marguerite  Bay,  the  near-inertial  energy  also  decays  with 

5Note  the  depths  of  the  Al  mooring  instruments  are  as  100,  309  and  459  m 
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Figure  3-13:  Rotary  spectra  of  the  de-tided  current  record  at  100  m  from  A2.  The 
vertical  dotted  line  indicates  the  inertial  frequency  (1.84  cpd,  «13-h)  and  the  solid 
lines  the  half-amplitude  periods  for  the  bandpass  filter  (1.2  cpd  and  3  cpd)  used  to 
filter  the  current  records.  For  this  record,  the  ratio  of  anticlockwise  to  clockwise 
energy  at  the  inertial  frequency  is  391.  Notice  also  the  residual  of  tidal  energy  at  the 
S2  frequency. 

depth  but  all  of  the  records  show  a  significant  inertial  peak,  including  the  deepest 
instruments  at  C2  and  C3.  While  most  moorings  show  either  a  decrease  of  near- 
inertial  energy  with  depth,  or  a  lower  layer  (below  200  m)  with  nearly  constant 
energy,  C3  shows  a  significantly  higher  peak  at  800  m  than  at  400  m.  As  with  the 
tidal  analysis  presented  earlier  in  this  chapter,  the  structure  at  C3  is  consistent  with 
the  co-located  B2  down  to  400  m,  the  deepest  current  meter  on  B2. 
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3.4.2  Band-pass  Filtering 

In  order  to  compute  basic  statistics  of  the  near-inertial  motions,  the  hourly  current 
meter  and  ADCP  data  were  filtered  with  a  cosine-lanczos  filter  (Emery  and  Thomp¬ 
son,  1998)  with  100  weights  to  band-pass  variability  between  8  to  20-h.  The  filter  was 
applied  to  the  de-tided  time  series,  as  the  semi-diurnal  tidal  components  are  close  to 
the  inertial  frequency  and  cannot  be  easily  filtered  out.  The  choice  of  8-h  to  20-h  for 
the  band-pass  filter  retains  the  broad  near-inertial  variability  while  filtering  out  any 
remnants  of  diurnal  tidal  variability  (Figure  3-13). 

3.4.3  Near-inertial  Currents 

Basic  statistics  of  the  near-inertial  band  motions  are  listed  in  Table  3.7  for  a  subset 
of  the  current  records.  Overall,  RMS  velocity  magnitudes  are  ~3  to  ~7.5  cm/s  for 
the  top  100-m  records  and  around  1  cm/s  for  the  deepest.  The  largest  velocities  are 
found  on  the  outer  shelf,  at  A3,  where  the  largest  velocity  recorded  was  26.1  cm/s. 
Near-inertial  band  motions  explain  a  minimum  of  10  to  15%  percent  of  the  total 
current  variance6  at  the  deepest  instruments  across  Marguerite  Bay  (B-  and  C-lines) 
and  as  much  as  65%  of  the  variance  at  103  m  at  A3.  Notice  that  even  at  400  m  in 
the  mid-  and  outer-shelf,  the  near-inertial  frequency  variability  explains  25-30%  of 
the  variance. 

The  principal  axis  of  the  band-pass  velocity  shows  that  the  polarization  (the  ratio 
between  the  major  and  minor  axis)  is  very  small  almost  everywhere,  with  the  excep¬ 
tion  of  Al  at  459  m  (Maj/Min  =  1.8),  which  also  shows  the  weakest  inertial  peak  (see 
previous  section).  This  is  dynamically  consistent  with  the  Wentzel-Kramers-Brillouin 
(WKB)  solutions  for  horizontal  currents  for  near-inertial  waves  (Kundu  and  Allen, 
1976;  Pedlosky,  2003): 


u 


(3.13) 


6At  most  instruments,  the  variance  explained  is  similar  for  both  velocity  components,  so  only  an 
average  is  presented. 
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°  =  ±A(i 

where  A  is  the  amplitude,  k  is  the  x-wavenumber,  m  is  the  vertical  wavenumber,  u  is 
the  frequency,  and  the  x-coordinate  is  aligned  with  the  direction  of  horizontal  prop¬ 
agation.  At  u>  —  /,  there  is  no  polarization  in  the  velocity,  and  as  the  frequency  of 
the  waves  increases,  the  principal  axis  becomes  increasingly  ellipsoidal.  This  relation¬ 
ship  is  not  very  consistent  in  our  velocity  records  (Figure  3-15),  but  the  near-circular 
motions  axe  generally  consistent  with  the  idea  of  near-inertial  waves. 

3.4.4  Vertical  Structure 

In  order  to  study  the  vertical  structure  of  the  near-inertial  motions,  complex  em¬ 
pirical  orthogonal  functions  (cEOF,  Emery  and  Thompson,  1998;  Shearman,  2005) 
were  calculated  for  each  mooring  using  all  the  available  data  that  covered  the  entire 
deployment  period  (with  the  exception  of  B3,  where  only  the  first  78  days  while  the 
ADCP  at  150  m  was  sampling  were  used). 

The  first  cEOF  explains  between  44  and  82%  of  the  variability  of  the  band-pass 
filtered  detided  data.  On  the  mid-  and  outer-shelf  (Figure  3-16,  top  panel),  the  A2 
and  A3  moorings  have  a  similar  vertical  structure  resembling  a  first  baroclinic  mode, 
with  the  velocity  at  100  m  180°  out  of  phase  with  the  velocities  at  depth.  Consistent 
with  results  mentioned  earlier,  the  first  cEOF  has  a  larger  amplitude  at  A3  than  at 
A2.  Near  the  shore  at  Al,  the  deepest  instruments  show  near-zero  amplitude,  and 
the  structure  in  the  top  160  m  is  complicated,  with  a  maximum  amplitude  of  4.4  cm/s 
at  around  125  m  and  counterclockwise  veering  with  depth  of  0.6-1.5°/m  over  the  top 
165  m. 

Across  the  mouth  of  Marguerite  Bay  (Figure  3-16,  bottom  panel),  the  structure 
of  the  first  cEOF  is  complicated.  The  moorings  that  include  data  in  the  top  100  m 
(B3,  Cl  and  C3)  show  counterclockwise  veering  with  depth,  and  velocity  amplitudes 
of  the  order  of  3-5  cm/s,  while  velocities  below  250  m  in  all  the  moorings  are  less 
than  1  cm/s  .  Unlike  the  clear  first-mode  baroclinic  structure  observed  at  A2  and 
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A3,  the  overall  vertical  structure  at  B3  and  C3  is  complicated,  with  the  small  bottom 
velocities  showing  a  ~90  °  phase  difference  with  the  surface  velocities. 

3.4.5  Horizontal  Scales 

Complex  correlation  at  zero  lag  (Kundu  and  Allen,  1976)  between  pairs  of  eigenvectors 
of  the  first  cEOF  modes  were  computed  in  order  to  get  an  order  of  magnitude  estimate 
of  the  horizontal  scales  of  the  near-inertial  motions7.  The  largest  correlation  (0.38)  is 
between  the  A2  and  A3  records,  but  neither  of  them  shows  significant  correlation  with 
Al,  nor  with  the  B-line  records.  The  B2  and  B3  records  are  significantly  correlated 
(0.35),  although  no  significant  correlation  was  found  between  the  C-line  moorings. 
This  is  surprising,  as  B2  and  B3  are  separated  by  roughly  the  same  distance  as  the 
C-line  moorings  and  located  roughly  in  the  same  area.  Notice,  however,  that  the 
bathymetry  separating  the  C-line  mooring  sites  is  quite  complicated,  with  at  least 
one  large  sea  mount  between  C1-C2  and  C2-C3.  Comparatively,  the  flow  has  less 
restrictions  between  B2/C3  and  B3  (Figure  3-3). 

3.4.6  Seasonal  Evolution 

A  rough  idea  of  the  time  variability  of  near-inertial  motions  on  the  shelf  can  be 
obtained  by  calculating  monthly  averages  of  the  kinetic  energy  per  unit  mass  (e.g. 
Shearman,  2005): 

KE=l-(cl  +  ol)  (3.14) 

where  a 2  is  the  variance  for  a  given  near-inertial  velocity  component  (calculated  from 

the  band- pass  filtered  detided  data).  Figure  3-17  (top  panel)  shows  KE  for  all  the 

available  current  records  at  100  m.  The  variability  from  month  to  month  is  largest 

on  the  mid-  and  outer-shelf  and  smallest  across  Marguerite  Bay,  especially  in  2002. 

The  largest  values  of  KE  were  registered  during  the  spring  of  2001  and  well  into  the 

summer  of  2002,  with  values  as  large  as  three  times  the  ones  for  the  fall  and  winter, 

7One  could  use  different  methods  here,  like  calculating  correlation  indices  for  individual  velocity 
components,  or  calculating  the  coherence  between  velocity  pairs.  Those  methods  result  in  the  same 
qualitative  results  presented  here. 
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but  spatially  limited  to  the  A-line  moorings,  with  no  significant  increase  at  B2,  some 
150  km  to  the  south.  In  March  of  both  years,  KE  is  relatively  high  everywhere. 

The  middle  panel  of  Figure  3-17  shows  the  SSM/I  satellite-derived  ice-concentration 
(see  Chapter  5  for  more  information)  for  A2  and  B2/C3.  Large  spatial  and  interan¬ 
nual  changes  are  evident.  In  2001,  the  ice  remained  in  Marguerite  Bay  long  after  it 
had  retreated  from  the  open  shelf,  and  there  were  only  brief  periods  of  no  ice  before 
the  next  winter  arrived.  The  interannual  changes  are  also  large,  with  the  ice  arriving 
roughly  two  months  earlier  in  2002  than  2001,  and  covering  a  larger  extent  of  the 
shelf. 

The  evidence  of  suppression  of  near-inertial  motions  by  the  ice  cover  is  mixed. 
In  2002  the  largest  KE  values  occurred  when  the  shelf  was  either  ice-free  (April) 
or  a  clear  opening  occurred  in  the  ice  (October).  In  2001,  the  highest  values  also 
occurred  during  relatively  ice-free  conditions,  but  there  are  also  relatively  high  values 
in  August  (although  notice  the  simultaneous  drop  in  ice  concentration  in  both  sites). 

The  standard  deviation  of  monthly  averages  of  the  eastward  and  northward  wind- 
stress  at  Kirkland  Island  (Figure  3-17,  bottom  panel)  shows  that  the  wind  is  highly 
variable  on  the  shelf,  with  no  clear  seasonal  signal  in  these  two  years,  and  with  strong 
month-to-month  variability.  The  correlation  between  the  wind  variability  and  the 
KE  in  the  upper  ocean  is  low,  which  is  particularly  evident  in  the  spring  of  2001, 
when  the  wind  variability  was  low  but  the  energy  in  the  A-line  moorings  was  highest. 


3.4.7  Discussion  and  Summary  of  Near-inertial  Variability 

The  brief  description  of  the  near-inertial  motions  presented  above  suggest  they  are  an 
important  component  of  the  wAP  circulation.  On  average,  the  near-inertial  motions 
tend  to  be  larger  than  the  most  energetic  tides,  and  single  events  can  have  velocities 
an  order  of  magnitude  larger. 

Several  of  the  moorings  also  show  a  shift  of  the  highest  peak  towards  frequencies 
larger  than  /  (Figure  3-14),  with  the  shift  increasing  with  depth  at  some  sites.  If  one 
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considers  the  vertical  component  of  the  group  velocity 


(iV2  —  o>2)(/2  —  a;2) 
mu(N2  —  f2) 


(3.15) 


In  an  idealized  ocean  with  no  lateral  boundaries  where  where  the  winds  force  inertial 
motions,  one  would  expect  no  vertical  propagation  ( u  =  /).  The  presence  of  a  coast, 
or  wind  forcing  with  a  finite  horizontal  scale,  for  example,  can  generate  vertically 
propagating  near-inertial  waves.  The  shift  in  the  near-inertial  peak  with  depth  is 
thus  consistent  with  the  dominant  waves  having  a  larger  vertical  component.  No¬ 
tice,  however,  that  these  shifts,  both  towards  frequencies  larger  and  smaller  than  /, 
have  been  observed  when  the  waves  propagate  through  flows  with  significant  relative 
vorticity  (Weller,  1982;  Shearman,  2005). 

These  results  also  highlight  the  small  dynamical  scales  and  strong  spatial  vari¬ 
ability  characteristic  of  this  area.  In  the  relatively  open  region  where  A2  and  A3 
are  located,  the  near-inertial  motions  are  well  correlated,  suggesting  scales  of  45  km 
(the  mooring  separation)  or  more  for  the  near-inertial  motions.  However,  across 
Marguerite  Bay,  the  C-line  moorings  show  little  correlation  among  them,  suggesting 
horizontal  scales  of  less  than  25  km.  The  correlation  between  the  first  cEOF  at  B2 
and  B3  is  high,  which  seems  strange  considering  the  results  for  the  C-line  moorings, 
which  show  small  correlation  among  them. 

The  evolution  of  the  100-m  monthly-mean  neax-inertial  kinetic  energy  shows  very 
high  values  for  the  spring  of  2001  and  early  summer  of  2002.  This  might  be  connected 
with  anomalous  atmospheric  conditions  with  strong  northwesterly  winds  (as  seen  from 
pressure  maps)  which  helped  create  the  high  ice  concentration  conditions  shown  here 
(Massom  et  al.,  2006),  and  highlights  the  need  for  wind  measurements  in  the  open 
shelf  region. 

The  forcing  mechanism  of  the  near-inertial  motions  remains  unclear.  On  monthly 
to  yearly  time-scales,  the  wind  forcing  does  not  seem  well  correlated,  although  further 
exploration  of  the  dataset  is  needed.  Moreover,  the  discussion  of  the  tidal  current 
variability  (Section  3.3.4)  suggest  that  internal  tides  might  be  common  on  the  wAP 
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shelf.  These  internal  waves  can  have  a  broad  spectral  response  to  the  barotropic  tidal 
forcing,  and  therefore  the  band-pass  filtering  used  to  calculate  statistics  (Section  3.4.2) 
would  not  remove  this  motions.  Therefore,  the  variability  discussed  here  is  likely  to 
be  a  mixture  of  wind-  and  tidally-driven  near-inertial  motions. 

Overall,  near-inertial  motions  are  an  important  component  of  the  flow  on  the 
wAP  shelf.  Because  these  motions  have  a  significant  decay  of  energy  in  the  vertical 
(which  results  in  vertical  shear)  and  are  a  relatively  important  component  of  the 
circulation,  near-inertial  motions  might  be  an  important  candidate  to  provide  the 
mixing  necessary  to  close  the  property  budgets. 
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Table  3.7:  Statistics  of  near-inertial  band  motions  on  the  wAP  shelf.  The 
columns  show  the  root  mean  square  (RMS),  maximum  band-pass  velocity 
measured  (Max),  the  variance  explained  by  the  band-pass  filtered  time 
series  (Var  -  see  text),  the  major  (Maj),  minor  (Min)  and  inclination  (Inc) 
of  the  principal  axis  of  the  band-pass  time  series,  and  the  record  length. 


Mooring 

Depth 

RMS 

Max 

Var 

Maj 

Min 

Inc 

Length 

M 

[cm/s] 

[cm/s] 

[%] 

[cm/s] 

[cm/s] 

[°] 

[days] 

A1 

103 

5.3 

17.4 

38 

3.8 

3.7 

18 

318 

309 

1.8 

7.8 

29 

1.3 

1.2 

3 

318 

459 

1.3 

8.6 

19 

1.1 

0.6 

-21 

318 

A2 

107 

5.3 

17.3 

48 

3.9 

3.7 

83 

318 

247 

2.7 

8.3 

27 

1.9 

1.8 

-77 

318 

397 

2.6 

8.8 

24 

1.8 

1.8 

-53 

318 

A3 

103 

7.4 

26.1 

65 

5.3 

5.1 

2 

318 

253 

3.1 

11.1 

31 

2.2 

2.2 

9 

318 

403 

3.2 

10.4 

31 

2.3 

2.2 

15 

318 

B2 

122 

3.8 

15.9 

37 

2.7 

2.7 

-6 

318 

254 

2.3 

6.8 

33 

1.6 

1.6 

75 

318 

400 

0.9 

3.5 

15 

0.7 

0.6 

85 

318 

B3 

51 

4.7 

11.9 

13 

3.5 

3.2 

-26 

79 

101 

5.7 

14.6 

40 

4.0 

4.0 

-79 

79 

247 

2.3 

8.1 

22 

1.7 

1.5 

23 

318 

400 

2.2 

8.1 

15 

1.7 

1.4 

-20 

318 

Cl 

50 

3.4 

12.7 

14 

2.5 

2.4 

-47 

365 

88 

3.5 

13.9 

21 

2.5 

2.4 

33 

365 
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Figure  3-14:  Counterclockwise  component  of  the  rotary  spectra  for  the  SO  GLOBEC 
moored  array  records,  calculated  from  the  de-tided  (not  band-passed)  data.  In  each 
panel,  the  vertical  dotted  line  is  the  local  inertial  frequency,  and  two  horizontal  lines 
are  plotted  as  reference.  The  horizontal  axis  spans  the  1-4  cpd  frequency  range. 
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Figure  3-15:  Ratio  of  major  to  minor  axis  of  variance  of  the  band-pass  velocity  records 
as  a  function  of  the  ratio  of  the  frequency  of  the  maximum  counterclockwise  spectral 
peak  to  the  inertial  frequency.  The  point  corresponding  to  the  deepest  instrument  at 
A1  (u/  f  =  1.01,  Maj/Min  =  1.88)  is  omitted.  The  dotted  line  is  the  1:1  relationship. 
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Figure  3-16:  First  cEOF  for  the  A-line  moorings  (top  panels)  and  the  B-  and  C-line 
moorings  (bottom  panels).  The  horizontal  line  indicates  the  bottom  depth  at  each 
mooring.  For  each  mooring,  the  percentage  of  the  variance  of  the  band-passed  records 
is  indicated. 
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Figure  3-17:  Monthly  averaged  kinetic  energy  at  100  m  (top  panel),  sea-ice  concen¬ 
tration  (middle  panel)  and  standard  deviation  of  the  monthly  averaged  wind-stress 
at  Kirkland  Island  (bottom  panel). 
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3.5  Bottom-trapped  Topographic  Waves  in  Mar¬ 


guerite  Trough 

Examination  of  the  vertical  structure  of  the  current  records  revealed  in  some  of  the 
moorings  bottom-intensified  subtidal  motions  in  scales  of  a  few  days  to  two  weeks  or 
so.  Figure  3-18  shows  the  power  spectra  of  the  along-  and  cross-slope  currents  for 
two  of  these  moorings,  A2  and  C2  (the  third,  B3,  shows  a  similar  structure).  At  A2, 
there  is  consistently  more  energy  at  the  deepest  instrument  (397  m)  than  at  247  m  in 
subtidal  scales  of  a  few  days  and  longer  than  9  days  or  so,  with  a  gap  in  this  pattern 
around  seven  days.  The  surface-most  instrument,  at  107  m  is  less  consistent,  and 
shows  significantly  more  energy,  in  particular  in  the  cross-slope  direction.  At  C2,  all 
the  records  show  energy  increasing  with  depth  (notice  this  mooring  had  a  failure  in  its 
surface- looking  ADCP  so  no  records  of  the  surface  currents  were  recovered).  Because 
the  presence  of  the  relatively  strong  bottom  slope  in  Marguerite  Trough  together 
with  the  deep  stratification  can  sustain  baroclinic  vorticity  waves  (e.g.  Rhines,  1970; 
Thompson  and  Luyten,  1976;  Auad  et  al.,  1998),  we  speculated  that  these  motions 
could  be  due  to  bottom-trapped  topographic  Rossby  waves  (TRWs). 

In  the  following  sections,  I  use  the  A2  mooring  records  to  explore  this  hypothesis. 
The  results  at  C2  are  qualitatively  similar,  but  the  applicability  of  the  simple  model 
outlined  below  is  at  best  questionable:  simple  theories  which  can  be  tested  with  a 
single  mooring  assume,  a  simplified  bathymetry,  and  the  C2  mooring  site,  surrounded 
by  the  walls  of  Marguerite  Trough  and  a  large  sea  mount,  means  one  must  make 
arbitrary  decisions  about  the  coordinate  system  and  direction  of  propagation.  The 
A2  mooring  site,  instead,  provides  a  more  straight-forward  test  of  the  hypothesis. 


3.5.1  Dynamics 

In  the  context  of  linear,  unforced  quasi-geostrophic  dynamics  in  the  presence  of  a 
small  bottom  slope  and  constant  stratification,  the  leading  order  horizontal  velocity 
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field  can  be  rewritten  using  a  streamfunction  (e.g.  Rhines,  1977): 

•  ei(kx+ly~wt)  cosh (nNz/f)  (3.16) 

where  N  is  the  buoyancy  frequency,  k  is  the  magnitude  of  the  wavevector,  /  is  the 
Coriolis  parameter  and  z  is  measured  positive  upward  from  the  bottom.  This  func¬ 
tion  is  a  solution  to  the  quasi-geostrophic  potential  vorticity  equation  with  boundary 
conditions: 


w  =  0  at  z  =  H 


dh 

w  =  u  •  — 
dy 


at  2  =  0 


where  h  is  the  bottom  bathymetry  and  y  is  the  upslope  direction.  The  dispersion 
relation  for  these  bottom  trapped  modes  is  given  by  (for  a  complete  derivation,  see 
Rhines,  1970): 


u  =  N  a  sm  coth.(KN  H  / /), 


sin  (j)  = 


k 

(/c2  +  Z2)1/2 


(3.17) 


where  a  is  the  bottom  slope,  H  is  the  water  depth,  (p  is  the  orientation  of  the  wavevec¬ 
tor  with  respect  to  the  upslope  direction  (and,  as  k  ■  u  =  0,  also  the  orientation  of  the 
velocity  vector  to  the  isobaths).  In  the  short-wave  limit,  where  the  internal  radius 
of  deformation  is  larger  than  the  scale  of  the  motion,  coth (kNH/ /)  «  1  and  the 
dispersion  relation  becomes: 


u  =  Na  sin  <p 


(3.18) 


that  is,  independent  of  the  wavevector  magnitude.  In  this  case,  the  highest  possible 
frequency  is  given  by  u  =  Na.  Notice  that  both  the  phase  speed  and  group  velocity 
propagate  to  the  right  facing  shallow  water  in  the  southern  hemisphere: 

Na  (k\  Na 

1  (&  +  p)v 2  ’  y  \i)(k2  +  py /2 
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NaP  _  Nakl 

9X  ~  (k2  + 12) 3/2  ’  9y  ~  (k2  +  P) 3/2 

In  the  case  of  Marguerite  Trough,  both  phase  and  energy  would  therefore  propagated 
clockwise,  towards  the  coast  on  its  northern  side  and  towards  the  ocean  on  its  south¬ 
ern  side.  Figure  3-19  depicts  the  vertical  structure  of  a  single  wave  packet  and  the 
orientation  of  the  wave  structure. 

In  summary,  the  assumptions  involved  in  the  derivation  of  the  above  dispersion 
relation  and  derived  quantities,  are: 

•  small  slope  ( L  =  1/k):  5  =  aL/H  <C  1 

•  linear:  Ra/5  =  UH/(fL2a)  1 

•  topographic  beta  larger  than  planetary  beta:  a/5  =  /30/(fa/H)  -C  1 

•  short  waves8:  \/L  =  NH/(fL )  <C  1 

•  low  frequency  motions:  u/ f  =  0(5) 

•  Constant  stratification 


A  few  properties  of  the  dispersion  relation  and  velocity  field  can  be  used  to  evaluate 
the  dynamics  described  above  with  single  mooring  data  (Thompson  and  Luyten, 
1976).  From  Equation  (3.16),  the  ratio  of  kinetic  energy  at  depths  Z\  and  z 2  is  given 
by: 


cosh(N /  f  kzi) 
cosh(N/ /kz2) 


(3.19) 


where  R,  if  significant,  can  be  evaluated  from  observations  to  obtain  the  horizontal 
scale  of  the  waves  k.  As  the  dispersion  relation  for  the  short-wave  limit  (Equation 
3.18)  is  independent  of  the  magnitude  of  the  wavevector,  it  provides  a  prediction  for 
its  orientation  (and  thus,  for  the  orientation  of  the  velocity  vector)  as  a  function  of 
frequency  for  a  given  stratification  and  bottom  slope.  This  theoretical  estimate  can 

8Not  strictly  required  unless  one  wants  to  use  (3.17)  as  the  dispersion  relation. 
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be  compared  to  the  orientation  of  the  principal  axis  of  the  velocity  records  (again,  as 
function  of  frequency).  Finally,  the  density  field  and  the  up-slope  velocity  component 
are  expected  to  be  in  quadrature,  with  the  velocity  leading. 

3.5.2  Topographic  Waves  at  A2 

As  described  in  section  3.5.1,  linear  quasi-geostrophic  theory  predicts  that  on  time 
scales  longer  than  ( aN )-1,  bottom  intensified  modes  can  be  excited  in  the  presence 
of  stratification  and  a  bottom  slope.  The  power  spectra  shown  in  Figure  3-18  show 
several  frequency  bands  with  such  vertical  decay  in  kinetic  energy.  Although  the 
vertical  decay  is  not  significant  in  the  spectra  plotted  (based  on  the  95%  confidence 
F-distribution  test,  see  Emery  and  Thompson,  1998),  band- averaging  of  the  spectral 
estimates  over  synoptic  time-scales  (7  to  20  days)  to  increase  the  number  of  available 
degrees  of  freedom  showed  that  the  decay  is  significant  on  average  over  that  frequency 
band. 

The  ratio  of  the  variance  for  the  along-slope  velocity  component  of  the  bottom  two 
instruments  is  0.71,  which  can  be  used  to  estimate  the  horizontal  scale  of  the  waves 
using  Equation  (3.19).  A  reasonable  estimate  of  the  lower-layer  stratification  at  A2 
can  be  obtained  from  the  density  records  collected  on  the  mooring  (which  axe  not 
significantly  different  from  CTD  casts  collected  at  or  near  the  mooring  site)  and  is, 
on  average,  N  =1.3T0~3  s-1.  The  local  bottom  slope  a  is  uncertain,  but  a  reasonable 
value  is  0.02  over  a  scale  of  50  km  or  so.  With  /  =  1.35  •  10~4  s-1,  the  horizontal 
scale  of  the  waves  is  k  =2. 1  -10— 4  m. 

Consider  now  the  orientation  of  the  velocity  vector.  At  low  frequencies,  motion 
is  predicted  to  be  mainly  along  isobaths,  with  the  velocity  slowly  turning  at  higher 
frequencies  until  we  reach  the  highest  allowable  frequency,  aN,  when  motions  are 
purely  up  and  down  the  slope.  In  order  to  compare  the  observations  to  the  theoretical 
prediction,  the  velocity  records  at  the  A2  mooring  were  band- pass  filtered  around  a 
central  period  that  was  varied  from  around  60  days  to  (aN)-1.  Figure  3-20  shows 
the  orientation  of  the  band-passed  velocity  records  at  the  three  measurement  depths. 
Also  plotted  is  the  theoretical  estimate  of  the  orientation  from  Equation  3.18  using 
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local  values  for  N  and  a. 

The  results  are  encouraging  for  the  bottom  two  instruments,  which  are  the  pairs 
that  shows  consistent  variance  decay  from  the  bottom.  The  misfit  of  the  theoretical 
estimate  and  the  100-m  velocity  record  suggests  that  the  waves  have  a  vertical  scale 
of  less  than  400  m,  or  more  likely,  that  the  large  increase  in  N  associated  with  the 
permanent  pycnocline  found  centered  near  150  m  reduces  the  wave  amplitude  in  the 
surface  layer  above  the  maximum  stratification.  The  flow  is  highly  coherent  between 
250  m  and  400  m,  with  zero  temporal  phase  lag  between  the  band-pass  currents 
currents  at  250  and  400  m  (not  shown),  as  is  predicted  by  linear  theory. 

An  independent  test  of  the  presence  of  this  waves  comes  from  the  relation  between 
the  density  and  the  upslope  velocity.  At  A2,  the  density  field  and  the  upslope  velocity 
are  roughly  in  quadrature  for  periods  around  10  days  (Figure  3-21),  as  expected  from 
the  theory,  with  significant  but  low  (0.4)  coherence. 

3.5.3  Discussion  and  Summary  of  Bottom-trapped  Topographic 
Waves 

Several  assumptions  are  made  in  the  derivation  of  the  theory  used  here.  A  number  of 
them  can  be  revisited  a  posteriori  and  are  shown  in  Table  3.8.  The  horizontal  scale 
of  the  waves  is  not  much  larger  than  the  radius  of  deformation,  so  in  this  analysis  the 
short-wave  assumption  is  the  weaker  one. 


Table  3.8:  Values  of  the  non-dimensional  parameters  rele¬ 
vant  to  the  derivation  of  the  linear  quasi- geostrophic  vorticity 
waves.  All  values  should  be  much  less  than  1. 


Site  (Period/Scale)a 

8  Ro/8 

a/S  L/\ 

“>// 

A2  (10.9d/30km) 

0.17  0.43 

1.8  •  10-3  0.88 

0.05 

0  Period:  Center  of  the  band  where  significant  decay  of  power  was 
found  (see  text).  Scale:  For  the  given  period,  L  =  l//t,  as  calculated 
from  Eq.  (3.19). 


Two  further  assumptions  are  made,  namely  constant  stratification  and  bottom 
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slope.  The  latter  is  somewhat  broken  in  the  study  area  as  the  slope  changes  signifi¬ 
cantly  in  the  WKB  sense.  A  derivation  of  TRW  theory  relaxing  this  assumption  or 
pursuing  a  numerical  solution  is  required  to  determine  the  effects  of  a  variable  slope 
on  these  waves.  As  for  the  constant  stratification  assumption,  this  can  be  explored 
by  solving  the  vertical  structure  of  the  quasi-geostrophic  potential  vorticity  equation 
numerically  for  a  given  stratification  function.  The  result  (not  shown)  indicates  that 
in  the  study  region  one  should  expect  the  TRW  modes  to  decay  extremely  fast  and 
not  to  penetrate  the  pycnocline,  found  at  150  m  or  so. 

The  above  results  assume  that  the  spectral  properties  do  not  change  in  time. 
An  idea  of  the  evolution  of  the  velocity  field  can  be  obtain  by  performing  a  wavelet 
analysis  of  the  A2  records  (Torrence  and  Compo,  1998).  The  wavelet  transform 
allows  for  the  energy  (or  variance)  to  be  decomposed  both  in  the  time  and  frequency 
domains.  Figure  3-22  shows  the  wavelet  variance  for  the  along-isobath  components 
of  the  velocity  at  A2,  averaged  over  periods  from  7  to  20  days,  which  corresponds 
closely  to  the  band  used  to  calculate  R.  The  result  is  interesting  because  Figure  3-22 
shows  the  bottom-intensified  events  to  be  highly  episodic.  There  is  one  large  event 
in  mid-August  to  mid-September  and  another  one  the  first  few  weeks  of  November. 
Moreover,  there  are  clearly  other  dynamics  acting  at  these  frequency  bands.  There 
is  a  high-energy  period  during  June,  and  another  period  during  October  that  shows 
some  interesting  vertical  structure.  This  complexity  is  not  too  surprising,  as  the  10  d 
time-scale  is  some  of  the  most  energetic  in  the  atmospheric  forcing. 

In  summary,  the  results  shown  here  are  consistent  with  the  presence  of  bottom- 
trapped  topographic  Rossby  waves  in  Marguerite  Trough.  At  A2,  the  waves  have  a 
wavelength  of  30  km  or  so  and  a  vertical  scale  spanning  several  hundred  meters  from 
the  bottom,  although  the  pycnocline  prevents  their  penetration  into  the  surface-most 
layer. 
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Figure  3-18:  Power  spectra  of  the  along-slope  velocity  (left  panel)  and  the  cross¬ 
slope  velocity  (right  panel)  for  the  A2  (top  panels)  and  C2  (bottom  panels)  mooring 
records. 
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Figure  3-19:  Schematic  showing  the  structure  of  a  bottom-trapped  topographic 
Rossby  wave.  The  top  panel  shows  the  velocity  structure  of  a  wave  over  a  bot¬ 
tom  slope  a,  and  with  a  vertical  velocity  structure  decaying  with  distance  from  the 
bottom.  The  bottom  panel  shows  the  velocity  vector  u  with  respect  to  the  wavevector 
k  and  the  group  velocity  Cg.  <t>  is  the  angle  between  the  upslope  direction  and  the 
wavevector.  Modified  from  Rhines  (1970) 
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Figure  3-20:  Orientation  of  the  band- pass  velocity  as  a  function  of  filter  frequency 
for  the  107  (dotted),  247  (dashed)  and  397  (solid  black)  A2  velocity  records.  The 
angle  is  positive  (negative)  in  the  upslope  (downslope)  direction.  The  prediction  of 
the  orientation  of  the  currents  records  from  wave  theory  (solid  gray)  is  also  plotted. 
The  vertical  line  corresponds  to  the  maximum  frequency  possible  for  the  short-wave 
limit  considered  (See  text). 
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Figure  3-21:  Coherence  (top  panel)  and  phase  (bottom  panel)  between  the  downslope 
velocity  v  and  the  density  p  at  247  m  at  the  A2  mooring.  A  phase  value  of  -90  ° 
indicates  that  the  density  lags  the  velocity  by  90  °,  as  predicted  by  the  theory  discussed 
in  the  text.  The  horizontal  dashed  line  is  the  95%  confidence  level  for  the  coherence. 
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Figure  3-22:  Wavelet  variance  averaged  over  periods  of  7  to  20  days  for  the  along- 
isobath  current  records  at  A2. 
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3.6  Alongshore  Wind  Response  in  Marguerite  Trough 


The  relationship  between  ocean  circulation  and  wind  forcing  in  this  region  has  not 
been  investigated  until  recently,  when  the  relationship  between  ice  motions,  wind, 
and  the  surface  ocean  at  near-inertial  frequencies  were  studied  (Hyatt  et  al.,  2006). 
The  subtidal  circulation,  on  the  other  hand,  has  not  been  studied.  This  section 
aims  at  providing  a  lowest-order  description  of  the  subtidal  current  wind  response  in 
Marguerite  Trough. 

3.6.1  Wind  Data 

As  part  of  SO  GLOBEC,  two  automatic  weather  stations  (AWSs)  were  deployed  on 
Kirkland  and  Dismal  Islands,  in  Marguerite  Bay  (see  Chapter  2  for  details).  The  AWS 
were  not  deployed  until  late  May  of  2001,  almost  two  months  after  the  beginning  of 
the  2001  mooring  deployment.  The  wind  speed  and  direction  were  used  to  calculate 
the  wind  stress  following  Large  and  Pond  (1981). 

Figure  3-23  shows  the  mean  and  principal  axes  of  the  wind  stress  at  the  two 
AWSs,  which  are  separated  by  roughly  30  km.  The  wind  stress  has  a  mean  vector 
roughly  parallel  to  the  direction  of  the  main  landmass  of  the  Peninsula,  and  both 
these  measurements  and  atmospheric  model  results  suggest  this  is  the  general  direc¬ 
tion  of  the  mean  throughout  this  region  (van  Lipzig  et  al.,  2004),  making  the  winds 
downwelling-favorable  for  most  of  the  region9.  The  wind  stress  at  the  two  sites  is 
highly  correlated.  Complex  correlation  at  zero  lag  (Kundu  and  Allen,  1976)  between 
the  wind-stress  time  series  at  Dismal  at  Kirkland  Islands  is  0.86  with  a  phase  angle  of 
7.3  °.  The  wind  stress  at  Dismal  has  principal  axes  which  are  markedly  more  polarized 
than  at  Kirkland,  and  a  mean  wind  almost  twice  as  large  (Table  3.9).  This  variability 
suggests  the  influence  of  the  surrounding  topography  and  of  the  complicated  coastline 
of  Marguerite  Bay,  which  appear  to  steer  the  wind  and  thus  implies  that  the  AWS 
wind  measurements  might  not  represent  the  wind-field  in  the  open  shelf.  Also,  the 
two  instruments  were  located  at  different  heights  on  small  islands  (see  Chapter  2). 

9Naturally,  along  the  southwest  coast  of  Marguerite  Bay,  the  winds  are  upwelling-favorable 
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These  two  factors  might  also  contribute  to  the  observed  differences  between  the  two 
AWS  records. 


Table  3.9:  East  and  north  means  and  principal  axes 
of  the  wind  stress  from  the  two  SO  GLOBEC  au¬ 
tomatic  weather  stations  on  Kirkland  and  Dismal 
Islands  and  from  the  Polar  MM5  model  output  at 
the  A2  mooring  site  (statistics  from  05/25/02  to 
02/13/02). 


Location 

Tx 

Nm“2 

Ty 

Nnr2 

Maj 

Nm-2 

Min 

Nm”2 

Inc. 

o 

Kirkland  I. 

-0.01 

-0.05 

0.12 

0.09 

20 

Dismal  I. 

0.02 

-0.09 

0.17 

0.09 

-62 

MM5.A2 

-0.02 

-0.06 

0.12 

0.09 

39 

The  apparent  local  effects  of  coastal  topography  on  the  AWS  winds  and  the  lack 
of  wind  during  the  first  few  months  of  the  mooring  deployment  led  us  to  find  other 
sources  of  wind  information  for  the  area.  Output  from  the  lowest  level  of  the  Polar 
MM5  atmospheric  model,  run  by  the  Polar  Meterology  Group  at  The  Ohio  State 
University  (Guo  et  al.,  2003),  were  obtained  by  M.  Caruso  at  WHOI.  The  model 
produces  a  3-day  forecast  of  the  wind  field  every  12  hours,  with  a  temporal  resolution 
of  three  hours,  on  a  grid  with  a  horizontal  resolution  of  60  km.  Comparisons  of  the 
model  output  with  the  available  shipboard  wind  data,  AWS  data,  and  the  scarce 
satellite  data  were  good  (M.  Caruso,  personal  communication).  Figure  3-23  shows 
the  mean  and  principal  axes  calculated  from  the  model  output  at  the  A2  mooring  site 
(MM5.A2),  for  the  same  period,  and  calculating  the  stress  using  the  same  formulation 
as  with  the  AWS  data.  The  model  wind  stress  show  a  magnitude  and  variability 
similar  to  Kirkland  Island  (with  complex  correlation  0.75),  but  it  is  more  closely 
aligned  with  the  coastline  of  Adelaide  Island  (Table  3.9)  and  are  also  downwelling- 
favorable  on  average,  although  the  variability  is  high  compared  to  the  mean  at  all  the 
sites.  Based  on  the  reasonable  agreement  found  between  model  and  the  observations 
by  M.  Caruso,  we  therefore  assume  that  the  MM5.A2  time  series  is  an  appropriate 
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representation  of  the  wind  stress  on  the  shelf  off  Adelaide  Island. 

In  the  rest  of  this  section,  the  response  of  the  currents  at  the  B-  and  C-line 
moorings  is  briefly  described  using  the  model  wind-stress  in  the  along-shelf  direction 
at  A2.  The  orientation  of  the  principal  axis  differs  from  the  local  coastline  orientation 
(f«620  counterclockwise  from  east),  which  is  the  direction  where  one  would  expect 
the  maximum  response  in  the  along-shelf  flow  in  the  case  of  a  much  simpler  two- 
dimensional  coastal  ocean  domain.  The  MM5.A2  wind-stress  was  therefore  rotated 
into  the  coastline  orientation,  which  is  then  called  ’along-shelf  wind-stress’.  The 
mooring  data  will  be  kept  in  the  same  along-  and  cross-isobath  directions  defined  in 
Section  3.2.1,  and  the  current  response  will  be  described  in  that  coordinate  system. 

3.6.2  Correlations  with  the  Wind-stress 

Overall,  cross-correlation  between  the  along-shelf  wind-stress  and  the  mooring  veloci¬ 
ties  is  low,  but  statistically  significant  in  several  locations.  The  95%  significance  level 
for  most  of  the  records  is  0.13-0.15,  and  the  cross-correlation  is  highest  at  A2  and  C2 
and  B3  in  the  along-isobath  direction  (0.40,  0.23  and  0.20  respectively)  with  phase 
lags  (wind  leading)  of  30-  to  50-h  (1  to  2  days).  The  correlation  at  all  these  moorings 
is  roughly  equal  with  depth.  The  other  moorings  show  no  significant  correlation  with 
the  along-shelf  wind-stress  in  the  along-isobath  direction.  For  the  cross-isobath  ve¬ 
locity,  the  highest  correlation  is  found  at  B2  (0.30),  Cl  (0.23),  A2  and  A3  («-0.20), 
with  phase  lags  of  the  same  order  as  for  the  along-shelf  flow. 

3.6.3  Magnitude  and  Structure  of  the  Response 

To  illustrate  the  response  of  the  circulation  in  Marguerite  Trough,  the  along-  and 
cross-isobath  subtidal  currents  records  were  bin-averaged  as  a  function  of  the  along- 
shelf  wind-stress.  For  all  the  data  falling  into  a  given  wind  range,  the  average  along- 
isobath  and  cross-isobath  velocity  are  calculated,  together  with  a  95%  confidence  level 
based  on  the  total  number  of  data  points  in  the  bin  and  the  integral  time-scale  of  the 
wind  for  the  entire  period  (f»45  h).  A  minimum  of  20  events  must  be  contained  in 
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each  bin  for  statistics  to  be  calculated  (the  usual  trade-off  is  to  lose  binning  data  at 
more  extreme,  but  rarer,  wind- stress  magnitudes). 

Figures  3-24  and  3-25  show  the  response  of  the  circulation  in  the  along-  and 
cross-isobath  direction  in  Marguerite  Trough.  The  range  of  along-shelf  wind-stress 
magnitudes  for  which  enough  events  were  found  is  -0.2  to  0.2  Nm'2.  Somewhat 
consistent  with  the  cross-correlation  results,  the  response  is  evident  in  a  subset  of 
the  moorings  only.  For  downwelling-favorable  wind-stress  (positive  in  the  coordinate 
system  used),  the  flow  is  forced  towards  the  shelf-break  (Figure  3-24)  and  downslope 
(Figure  3-25)  at  A2  and  only  along  the  isobath  at  C2  (the  northern  side  of  the  trough), 
with  maximum  velocities  of  2-4  cm/s  at  A2  and  1-2  cm/s  at  the  deeper  C2.  Notice 
that  there  is  a  hint  of  a  baroclinic  structure  at  A2,  with  the  vertical  shear  increasing 
with  the  wind-stress  magnitude  (the  uncertainties,  however,  are  too  large  for  this 
tendency  to  be  significant).  All  of  the  velocity  records  at  C2  are  from  well  below  the 
pycnocline  and  the  response  to  the  wind  is  barotropic  to  within  our  error  estimates. 
For  upwelling- favorable  wind-stress,  the  opposite  is  true,  with  both  of  the  moorings 
showing  shoreward  (negative)  flow  and  a  tendency  for  upslope  flow  at  A2. 

The  results  at  A3  have  large  uncertainties,  and  there  is  no  clear  response  in  the 
along-isobath  direction,  but  downwelling-favorable  wind-stress  tends  to  generate  pos¬ 
itive  flow,  which  is  roughly  in  the  along-shelf  direction  and  surface  intensified.  The 
steering  of  the  flow  in  the  along-isobath  direction  seen  at  A2  does  not  affect  the 
flow  at  A3,  most  likely  because  the  slope  is  much  gentler  on  the  southern  side  of 
the  trough  where  A3  was  deployed  (see  Figures  3-3  and  2-5).  For  upwelling-favorable 
wind-stress,  the  tendency  is  for  the  flow  to  be  downslope  (to  the  north-east),  although 
the  uncertainties  are  too  large  to  make  a  definitive  statement  about  this. 

3.6.4  Summary  of  the  Subtidal  Current  Wind  Response  in 
Marguerite  Trough 

When  the  winds  are  downwelling-favorable,  the  results  indicate  that  the  circulation 
is  anticyclonic  along  the  northern  side  of  Marguerite  Trough  and  southwestwards 
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at  A3  (roughly  in  the  along-shelf  direction).  Upwelling-favorable  wind-stress  causes 
shoreward,  upslope  How  at  A2  and  into  Marguerite  Bay  at  C2.  The  response  is  of 
the  scale  of  1  to  4  cm/s,  which  is  comparable  to  the  mean  flow  at  several  moorings 
(See  section  3.2).  Most  of  the  long  records  considered  here  were  located  deep  in  the 
water  column,  although  the  surface-most  instruments  both  at  A2  and  A3  show  some 
tendency  for  baroclinicity.  At.  C2.  however,  instruments  spanning  550  m  (all  below 
the  pycnodine)  show  a  clear  barotropic  response. 

Interestingly,  the  mean  flow  in  the  canyon  was  shown  (See  section  3.2)  to  be 
strongest  at  A2  and  cyclonic  around  the  canyon,  that  is.  in  the  opposite  sense  of 
what  would  be  expected  from  a  mean  downwelling-favorable  wind-stress  as  shown 
here. 
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Figure  3-23:  Mean  and  principal  axes  of  the  wind  stress  from  the  AWSs  at  Kirkland 
and  Dismal  Island  and  from  the  output  of  the  Polar  MM5  model  at  the  A2  mooring 
site.  The  statistics  are  calculated  for  the  period  starting  05/25/02  (the  beginning  of 
the  AWS  time-series)  until  02/13/02  (when  the  A-  and  B-line  moorings  were  recov¬ 
ered). 
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Figure  3-24:  Along-isobath  current  response  to  along-shelf  wind-stress  in  Marguerite 
Trough.  The  vertical  axis  is  the  along-isobath  velocity  binned  as  a  function  of  the 
along-shelf  wind  stress.  The  error  bars  are  the  95%  confidence  intervals  for  each  bin. 
Downwelling-favorable  wind-stress  is  positive. 
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Figure  3-25:  Cross-isobath  current  response  to  along-shelf  wind-stress  in  Marguerite 
Trough.  The  vertical  axis  is  the  cross-isobath  velocity  binned  as  a  function  of  the 
along-shelf  wind  stress.  The  error  bars  are  the  95%  confidence  intervals  for  each  bin. 
Downwel ling-favorable  wind-stress  is  positive. 
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Chapter  4 


Circumpolar  Deep  Water 
Intrusions 

4.1  Background 

Early  studies  of  the  western  Antarctic  Peninsula  (wAP)  shelf  revealed  a  complex 
hydrographic  structure  with  a.  strong  seasonal  cycle.  The  surface  layer  is  occupied 
by  Antarctic  Surface  Water  (AASW).  a  relatively  fresh  and  cold  water  mass  that 
is  usually  found  all  around  the  Antarctic  continent  (Figure  4-1).  During  the  late 
fall  and  winter,  the  surface  mixed  layer  is  deepened  by  surface  cooling  and  by  brine 
rejection  as  a  result  of  ice  formation.  The  resulting  mixed  layer  is  roughly  at  freezing 
temperatures  and  reaches  a  depth  of  100  m  or  so.  As  the  air  temperature  increases  at 
the  end  of  winter,  the  ice  breaks  up,  melt  water  from  the  coast  increases  and  a  posit  ive 
heat  flux  into  the  ocean  all  lead  to  rest  ratification  of  the  deep  mixed  layer,  leaving  a 
layer  of  cold  water  at  100  in  or  so  called  Winter  Water  (WW).  Below  the  WW.  the 
shelf  is  occupied  by  a  warm  and  salty  water  mass  called  modified  Circumpolar  Deep 
Water  (inCDW),  a  slightly  colder  and  fresher  version  of  the  oceanic  Circumpolar 
Deep  Water,  found  along  the  shelfbreak  of  the  wAP. 

Although  precise  estimates  of  the  heat  budget  are  lacking,  hydrographic  surveys 
and  modeling  studies  suggest  that  the  integrated  heat  budget  for  the  surface  layer 
results  in  a  net  heat  loss  to  the  atmosphere.  As  heat  is  readily  available  in  the  lower 
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Figure  4-1:  Evolution  of  the  hydrographic  structure  on  the  mid-shelf  of  the  wAP  (A2 
mooring  site).  The  panel  shows  a  6- S  diagram  of  profiles  taken  throughout  2001  and 
2002.  The  black  boxes  illustrate  the  typical  range  of  the  Antarctic  Surface  Water 
(AASW),  Winter  Water  (WAV)  and  Circumpolar  Deep  Water  (CDW).  The  dashed 
blue  line  indicates  the  freezing  point  of  seawater.  Notice  that  CDW  is  composed  of 
a  warm  branch  (UCDW,  with  temperatures  >1.5  °C)  which  is  evident  in  the  July 
profile,  and  a  lower,  saltier  branch  (salinities  larger  than  34.72  or  so). 


layer  of  the  shelf  in  the  form  of  mCDW,  and  because  this  water  mass  has  clearly  lost 
heat  after  moving  onto  the  shelf,  it  has  been  hypothesized  that  a  vertical  flux  from 
the  lower  layer  is  responsible  for  balancing  the  net  heat  loss  to  the  atmosphere.  The 
magnitude  of  this  upward  flux  has  been  estimated  from  observations  using  a.  simple 
two-dimensional  model  equating  lateral  to  vertical  fluxes,  resulting  in  an  estimate  of 
1.5  W/m2  (Klinck.  1998),  and  from  a  1-D  mixed-layer /ice  model  which  resulted  in 
an  estimate  closer  to  5  W/m2  (Smith  and  Klinck,  2002). 

In  order  to  close  the  heat  budget,  an  onshore  flux  of  heat  must  compensate  the 
cooling  of  the  lower  layer  of  the  shelf.  Hydrographic  data  from  a  summer  cruise  con¬ 
ducted  in  1993  suggested  that  CDW7  intruded  in  four  specific  locations  along  the  shelf, 
two  of  which  are  within  the  SO  GLOBEC  study  region  (Prezelin  et  ah,  2000:  Prezelin 
et  al.,  2004).  Remnants  of  CDW'  were  found  along  Marguerite  Trough  during  the  SO 
GLOBEC  cruises,  which  also  suggested  that  this  bathymetric  feature,  connecting  the 
shelfbreak  to  Marguerite  Bay,  provides  a  path  for  CDW  intrusions  to  move  across 
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the  shelf  (Klinck  et  al.,  2004).  On  the  basis  of  these  hydrographic  surveys,  it  has 
been  suggested  that  CDW  moves  into  Marguerite  Trough  and  across  the  shelf  into 
Marguerite  Bay  four  to  six  times  a  year  (Klinck  et  ah,  2004). 

In  tliis  chapter,  the  questions  of  the  spatial  scales,  frequency  and  circulation  of 
CDW  intrusions  onto  the  shelf  are  addressed  using  the  SO  GLOBEC  broad-scale 
hydrographic  cruises  and  moored  array  (See  Chapter  2  for  details  on  the  dataset).  In 
the  next  sections,  I  provide  a  characterization  of  the  source  waters  for  the  intrusions,  a 
definition  for  the  intrusions,  a  characterization  of  their  spatial  scales,  frequency,  and 
intensity,  and  a  rough  estimation  of  their  contribution  to  the  horizontal  heat  flux. 
Finally,  the  implications  of  this  work  for  the  wAP  system  and  biological  communities 
are  discussed. 


4.2  Characteristics  of  Oceanic  Waters 

A  hydrographic  section  conducted  in  the  austral  summer  of  1992  along  the  S04P 
World  Ocean  Circulation  Experiment  (WOCE)  line  illustrates  the  key  hydrographic 
structure  of  the  region.  The  section  was  conducted  starting  in  the  mid-shelf  re¬ 
gion  and  moving  offshore,  with  four  relatively  close  CTD  casts  across  the  shelfbreak 
(Figure  4-2).  Figure  4-3  shows  the  potential  temperature  (0),  salinity  and  potential 
density  (as)  sections  from  the  cruise.  Away  from  the  shelf,  the  vertical  structure  is 
characterized  by  a  deep  temperature  maximum  at  500  m  or  so.  corresponding  to  the 
core  of  the  Upper  Circumpolar  Deep  Water  (UCDW).  The  colder  and  saltier  water 
mass  below,  with  densities  >27.8  and  salinity  maximum  around  34.72  corresponds  to 
Lower  Circumpolar  Deep  Water  (LCDW),  which  is  also  characterized  by  higher  dis¬ 
solved  oxygen  values.  Both  branches  of  Circumpolar  Deep  Water  are  carried  by  the 
Antarctic  Circumpolar  Current  (ACC)  along  the  shelfbreak,  and  as  is  characteristic 
of  the  southern  boundary  of  the  ACC  around  Antarctica,  both  the  isotherms  and 
isohalines  shoal  towards  the  continent  (Sievers  and  Nowlin,  1984:  Orsi  et  ah,  1995). 

Because  the  Circumpolar  Deep  Water  is  carried  by  ACC  fronts  (see  Chapter  1 
for  a  brief  description  of  the  regional  circulation)  that  lack  a  distinctive  upper  layer 
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Figure  4-2:  Location  of  the  first  ten  CTD  casts  conducted  during  the  S04P  cruise 
(blue  squares),  the  SO  GLOBEC  moorings  (red  diamonds),  CTD  stations  conducted 
off  the  shelf  during  the  SO  GLOBEC  broad-scale  cruises  (open  circles)  and  location 
of  CTD  casts  from  the  cross-isobath  section  conducted  by  the  L.  M.  Gould  during 
the  fall  of  2001  (magenta  dots). 


signature  that  can  be  used  to  trace  its  presence  on  the  shelf,  surface  data  -  including 
satellite  images  -  are  of  little  use  in  determining  the  boundary  of  water  masses  carried 
by  the  ACC  near  the  wAP  shelf,  and  therefore  of  the  UCDW  intrusions  of  interest 
here.  Below,  we  use  the  hydrographic  survey  data  from  the  SO  GLOBEC  cruises  to 
provide  a  suitable  definition  of  an  UCDW  intrusion. 


4.2.1  Definition  of  an  Oceanic  Intrusion 

Although  the  SO  GLOBEC  sampling  grid  included  only  a  handful  of  CTD  casts  in 
the  open  ocean  during  each  of  the  broad-scale  surveys,  together  they  provide  a  good 
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Figure  4-3:  Color  contours  of  potential  temperature  (top  panel)  and  salinity  (bottom 
panel)  for  the  S04P  section,  conducted  between  February  22  and  24,  1992.  Contours 
of  oo  (black)  are  superimposed  in  both  panels.  The  white  vertical  lines  represent 
where  CTD  stations  were  taken. 
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picture  of  the  variability  of  the  source  water  for  the  deep  water  intrusions.  Figure 
4-4  shows  a  0- S  diagram  of  all  the  C’TD  casts  collected  off  the  shelf  and  shelf-break 
at  locations  deeper  than  2000  m.  The  figure  reveals  the  basic  features  of  the  seasonal 
changes  in  the  region.  The  relatively  warm  and  fresh  water  found  in  the  surface 
layer  during  the  fall  is  replaced  by  cold  (near  freezing)  and  salty  water  during  the 
winter.  At  depth,  the  variability  is  not  seasonal  but  mostly  spatial.  The  mid-depth 
potential  temperature  maximum  corresponding  to  the  core  of  the  UCDW  (Figure  4-4, 
right  panel)  is  characterized  by  a  relatively  wide  range  of  temperatures  (from  1.55 
to  2.1  °C)  and  salinities  of  34.62  to  34.68.  Below  UCDW,  the  6- S  properties  tend  to 
collapse  as  the  salinity  increases  and  temperature  decreases.  The  maximum  salinity 
of  around  34.72  and  temperatures  of  1.2  °C  or  so  corresponds  to  the  core  of  the  other 
branch  of  CDW,  Lower  Circumpolar  Deep  Water  (LCDW). 
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Figure  4-4:  0-S  diagram  of  all  the  CTD  casts  from  open  ocean  locations  collected 
during  the  fall  2001  (green  dots),  winter  2001  (red),  fall  2002  (blue)  and  winter  2002 
(magenta)  SO  GLOBEC  broad-scale  cruises.  Open  ocean  is  defined  as  locations 
deeper  than  2000  in.  The  left  panel  shows  the  full  depth  (the  dashed  line  is  the 
freezing  temperature)  and  the  right  panel  is  a  detail  of  the  0-S  space  associated  with 
Circumpolar  Deep  Water  (the  dashed  line  marks  1.5  °C). 


The  vertical  hydrographic  structure  on  the  shelf  is  qualitatively  similar  to  the 
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open  ocean,  except  for  a  weaker  subsurface  temperature  maximum  as  a  result  of  the 
modification  of  UCDW  into  inCDW  and  also  the  absence  of  LCDW  from  most  of 
the  shelf,  although  it  has  been  reported  as  a  thin  deep  layer  in  Marguerite  Trough 
(Klinck  et  al.,  2004).  Therefore,  a  reasonable  detection  criteria  for  UCDW  intrusions 
on  the  shelf  is  water  with  potential  temperature  of  1.5  °C  or  more  below  the  mixed 
layer.  This  is  consistent  with  what  is  known  about  the  source  waters  carried  by 
the  ACC  close  to  the  shelfbreak  from  observations  (Figure  4-4).  and  with  historical 
characterizations  of  the  ACC  water  mass  structure  and  boundaries  as  seen  in  Figure 
4-3  and  in  the  literature  (Orsi  et  al.,  1995;  Smith  et  al.,  1999). 

To  characterize  the  intrusions,  we  first  look  for  evidence  of  warm  water  in  the 
hydrographic  cruises  and  the  moored  array,  with  the  objective  of  determining  their 
amplitude,  spatial  scale,  frequency  and  circulation. 


4.3  Spatial  Distribution  of  UCDW  Events 

4.3.1  UCDW  in  the  Hydrographic  Surveys 

Evidence  of  UCDW  on  the  shelf  can  be  found  in  the  four  SO  GLOBEC  cruises.  The 
maximum  0  below  the  mixed  layer  (Figure  4-5)  shows  water  of  1.75  °C  or  more  along 
the  shelfbreak,  except  for  the  southern  half  of  the  study  region  during  the  fall  2001 
survey,  when  no  stations  were  conducted  there.  The  1.5  °C  isotherm,  representing 
relatively  unmixed  UCDW,  was  found  all  along  the  shelfbreak  and  on  the  outer  shelf 
for  all  cruises.  UCDW  was  found  penetrating  onto  the  shelf  at  two  locations:  during 
both  of  the  2001  survey  and  the  fall  2002  survey  around  G7.5°  S.  and  in  all  of  the 
cruises  around  67°  S.  along  the  eastern  side  of  Marguerite  Trough.  Warm  water  at 
these  two  sites  was  also  reported  in  previous  work  examining  data  collected  during 
hydrographic  surveys  conducted  in  1993  and  1994,  which  also  identified  two  other 
intrusion  sites  to  the  north  of  the  SO  GLOBEC  study  region  (Prezelin  et  al..  2000: 
Prezelin  et  al.,  2004). 

Broad-scale  hydrographic  surveys  sampling  the  shelf  with  a  typical  horizontal 
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Figure  4-5:  Maximum  potential  temperature  6  below  the  mixed  layer  for  the  four 
broad-scale  SO  GLOBEC  cruises.  The  nominal  CTD  horizontal  resolution  is  20  km. 
The  isotherms  separation  is  .25  °C.  Notice  the  figure  is  meant  to  show  presence  or 
absence  of  UCDW  (defined  here  as  8  >  1.5),  and  should  not  be  interpreted  as  repre¬ 
senting  contiguous  warm  intrusions  on  the  shelf  (see  text  for  details). 


resolution  of  20  km  and  showing  the  presence  of  UCDW  on  the  shelf  as  shown  in 
Figure  4-5  must  be  considered  carefully.  Figure  4-6  shows  a  detail  of  0-S  space  for 
the  CTD  casts  on  the  shelf,  showing  the  stations  where  UCDW  was  found.  In  most 
of  them,  the  signal  is  relatively  weak,  with  temperatures  between  1.5  and  1.6°C. 
In  the  fall  of  2001.  for  example,  the  only  station  where  water  warmer  than  1.6  °C 
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was  recorded  was  located  some  20  kin  from  the  shelf-break  on  the  eastern  side  of 
Marguerite  Trough.  During  the  winter  2001  cruise,  the  highest  temperature  recorded 
(>1.8°C)  on  the  shelf  was  also  located  on  the  eastern  side  of  the  trough,  between  the 
Al  and  A2  mooring  sites.  Another  isolated  patch  of  warm  water  can  also  be  seen  in 
the  winter  2002  cruise  (Figure  4-5),  as  well  as  in  the  1993  survey  reported  by  Prezelin 
et  al.  (2004). 


1  i _ _ _ jel  a _ cn _ i  i _ »  ac.  ■  _ j 
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Figure  4-6:  8- S  diagram  of  the  hydrographic  data  for  all  four  broad-scale  hydrographic 
cruises.  The  colors  represent  CTD  casts  collected  off-shelf  (gray),  on  the  shelf  showing 
UCDW  (red),  and  all  other  casts  on  the  shelf  (green). 

Although  UCDW  was  not  found  in  Marguerite  Bay,  there  is  more  evidence  Mar¬ 
guerite  Trough  is  a  path  for  the  intrusions.  Four  CTD  stations  collected  across  the 
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trough  by  the  L.  M.  Gould  during  tier  fall  2001  cruise  show  slightly  warmer  (though 
colder  than  1.5  °C)  and  salty  water  on  the  eastern  side  of  the  trough  (Figure  4-7), 
which  is  consistent  with  the  advection  of  warm  UCDW  and  LCDW  by  the  cyclonic 
flow  in  the  trough  (see  Chapter  3).  This  is  consistent  with  the  finding  of  LCDW  in 
the  deepest  part  of  Marguerite  Trough  -  and  within  Marguerite  Bay  -  reported  by 
Klinek  et  al.  (2004)  using  the  SO  GLOBEC  broad-scale  hydrographic  data. 


Figure  4-7:  Color  contours  of  potential  temperature  (left  panel)  and  salinity  (right 
panel)  for  a  section  across  Marguerite  Bay,  conducted  in  the  fall  of  2001.  Contours 
of  a,?  (black)  are  superimposed  on  both  panels.  The  figure  is  oriented  so  that  the 
eastern  edge  of  the  trough  is  on  the  right. 


4.3.2  UCDW  in  the  Moored  Array 

The  SO  GLOBEC  moored  array  (diamonds  in  Figure  4-5)  is  well  located  to  test,  the 
hypothesis  that  Marguerite  Trough  is  one  of  the  locations  for  UCDW  to  move  onto 
the  shelf.  The  A3  mooring  is  located  38  km  onshore  of  the  shelfbreak,  on  the  western 
side  of  Marguerite  Trough  (Figure  4-2).  A2  is  located  45  km  farther  onshore  (83  km 
from  the  shelfbreak),  on  the  eastern  side  of  the  Trough,  a  location  which  Figure  4-5 
suggests  is  a  path  for  the  intrusions.  The  analysis  shown  below  reveals  this  is  indeed 
the  case. 
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Inspection  of  the  detided  temperature  records  revealed  intense,  high-frequency 
events  of  warm  water  at  the  A2  and  A3  mooring  sites.  As  many  of  the  events  have  a 
duration  around  the  33-hr  half-amplitude  period  used  to  low-pass  filter  the  mooring 
records  elsewhere  in  this  thesis,  the  time  series  are  treated  differently  here.  All  of 
the  mooring  data  were  first  de-tided  using  T_TIDE  (Pawlowicz  et  al.  2002)  and 
subsequently  low-pass  filtered  with  a  100- weight  cosine- lanczos  filter  with  a  20- hr  half- 
amplitude  period  to  filter  out  the  inertial  frequency  variability  (Emery  and  Thompson, 
1998).  20-hr  is  a  good  compromise  between  filtering  the  broad  inertial  peak  centered 
at  12.9-hr  while  keeping  most  of  the  variability  associated  with  the  warm-water  events. 
Although  there  is  no  significant  tidal  signal  at  250  in  in  the  temperature  records  at 
either  mooring,  the  shallower  instruments  did  have  a  tidal  signal  removed  prior  to  the 
filtering.  All  of  the  records  here  are  for  the  entire  deployment  except  for  the  200  m 
temperature  record  at  A3,  which  is  246  days  long  as  the  temperature  recorder  failed 
on  12/03/01. 

At  A 2  and  A3,  water  warmer  than  1.5 °C  was  detected  in  instruments  located 
from  200  m  to  400  m.  Histograms  of  potential  temperature  for  the  instruments  where 
intrusions  where  detected  (Figure  4-9)  clearly  show  a  higher  prevalence  of  warmer 
than  1.5  °C  observations  at  the  A2  mooring  at  all  depths.  The  colder  anomalies 
at  A2.  particularly  at  200  m.  are  typical  of  the  moorings  closer  to  the  coast.  At 
250  m.  where  the  mid-depth  temperature  maximum  is  usually  found  on  the  shelf, 
temperature  higher  than  1.5  °C  was  recorded  2%  of  the  time  at  A3,  compared  to  15% 
at  A2.  At  Al,  two  brief  events  were  recorded  at  the  end  of  July,  when  the  APCC  is 
absent  (see  Chapter  2),  but  the  rest  of  the  time  series  shows  no  evidence  of  UCDW. 
None  of  the  remaining  moorings  (the  B-  and  C-lines)  recorded  waters  warmer  than 
1.5  °C  at  depth. 

In  order  to  calculate  the  statistics  of  the  events,  individual  events  were  defined  as 
follows,  for  each  instrument  recording  temperature  at  each  mooring,  times  when  the 
potential  temperature  was  1.5 °C  or  more  were  found.  The  beginning  (end)  of  the 
event  at  each  depth  was  defined  as  the  time  when  the  temperature  started  (finished) 
increasing  (decreasing).  An  example  of  one  such  event  is  shown  in  Figure  4-8.  If 
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Figure  4-8:  Example  of  a  UCDW  intrusions  as  recorded  at  the  A2  mooring  at  250  m. 
The  horizontal  dotted  line  indicates  the  mean  potential  temperature  for  the  entire 
record,  and  the  basic  scales  (AO.  Duration)  are  also  indicated. 

those  times  overlapped  with  events  recorded  at  another  instrument,  it  was  assumed 
they  corresponded  to  one  event  spanning  the  entire  depth  range  covered  by  those 
instruments. 


4.4  Frequency  and  Duration  of  Warm  Events 

Statistics  of  the  events  reveal  the  basic  features  of  the  intrusions  of  UCDW  onto  the 
shelf.  At  A2,  40  events  were  recorded  in  the  10.5-month  time  series,  compared  to  only 
11  events  at  A3.  At  both  moorings  sites,  the  events  were  relatively  short,  typically 
lasting  1  to  3  days,  with  events  as  long  as  7  days  at  A2  (Figure  4-10).  There  does  not 
seem  to  be  a  significant  seasonal  signal  in  the  occurrence  of  events.  At  A2,  there  are 
as  few  as  zero  events  per  month  (in  October)  to  seven  (in  May),  with  an  average  of 
four  events  per  month  (Figure  4-11).  At  A3,  all  the  events  occur  during  the  first  few 
months  of  the  deployment,  from  April  to  June,  resulting  in  a  similar  average  of  events 
for  the  fall  season  as  at  A2.  From  July  on,  no  events  of  warm  UCDW  were  recorded 
at  A3.  However,  it  must  be  noted  that  during  the  winter  2002  cruise,  UCDW  was 
found  at  the  location  of  the  A3  mooring  (Figure  4-5),  which  suggest  that  intrusions 
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Figure  4-9:  Histogram  of  potential  temperature  for  the  A2  (red)  and  A3  (blue)  moor¬ 
ings  at  approximately  200  m  (top  panel),  250  m  (middle  panel)  and  400  m  (bottom 
panel).  The  200  m  temperature  record  at  A3  is  24G-d  long  (the  instrument  failed  on 
12/03/01)  and  the  rest  are  ~318-d  long. 


at  this  location,  although  less  frequent  than  at  A2,  can  also  occur  during  different 
seasons. 


4.5  Circulation  Associated  with  the  Warm  Events 


The  A2  mooring  recorded  currents  at  100.  250  and  400  m.  Figure  4-12  shows  time 
series  of  potential  temperature  at  250  m  (where  the  signal  of  the  warm  events  is 
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Figure  4-10:  Histogram  of  the  UCDW-event  duration  for  the  A2  (red)  and  A3  (blue) 
moorings.  Also  indicated  are  the  total  number  of  events  (n)  for  each  mooring. 


Figure  4-11:  Number  of  intrusions  per  month  at  A2  (red)  and  A3  (blue)  throughout 
the  2001/2002  SO  GLOBEC  mooring  deployment  (March  31,  2001  to  February  12, 
2002). 

strongest)  and  along-isobath  and  cross-isobath  velocity  records.  Times  when  intru¬ 
sions  were  detected  at  any  depth  are  marked  in  red.  There  is  not  a  clear  relationship 
between  the  velocity  field  and  the  occurrence  of  warm  events.  At  both  moorings, 
correlations  between  u  or  v  and  0  are  small  (maximum  <0.20  with  lags  of  one  to 
two  days)  and  statistically  not  significant  at  the  95%  confidence  level.  The  subtidal 
current  variability  has  its  strongest  component  in  the  5-10  days  band  (see  Chapter 
3)  and  the  coherence  with  temperature  oscillations  in  the  1-3  day  band  where  the 
intrusions  occur  is  not  significant.  However,  it  is  important  to  remember  that  both 
cross-correlation  and  coherence  analysis  average  the  variability  over  the  entire  time 
series,  and  might  not  be  appropriate  when  studying  relatively  isolated  events. 

Average  velocities  during  the  warm  events  are  shown  in  Figure  4-13.  All  but  three 
(93%)  of  the  events  have  a  mean  average  along-isobath  velocity  toward  the  coast,  with 
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Figure  4-12:  Time  series  of  det.ided,  low-pass  filtered  time  series  of  potential  temper¬ 
ature  (top  panel),  along-isobath  velocity  (middle  panel)  and  cross-isobath  velocity 
(bottom  panel).  The  warm  UCDW  events  are  indicated  in  red. 


an  overall  mean  of  around  5.8T2.7  cm/s  (±1  standard  deviation).  The  magnitude 
and  vertical  structure  of  the  velocity  during  the  events  are  very  similar  to  the  average 
structure  during  the  remaining  part  of  the  records.  The  mean  velocity  during  the 
events  has  very  small  vertical  shear  with  a  total  of  10°  counterclockwise  veering  with 
depth.  For  all  the  intrusions,  the  flow  is  dominated  by  the  low-frequency  component 
and  the  velocities  in  the  frequency  band  of  the  UCDW  events  are  typically  less  than 
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10%  of  the  low-frequency  flow. 

The  velocity  and  temperature  records  from  the  moored  array  are  illuminating 
when  interpreting  the  broad-scale  survey  data  (Figure  4-5).  The  time  series  indicate 
that  warm  events  occur  frequently  and  are  short-lived.  This  is  consistent  with  the 
survey  data  if  one  considers  some  important  aspects  of  the  survey  resolution  and 
the  resulting  hydrographic  fields.  First,  the  warm  water  found  during  the  survey  is 
associated  with  different  density  surfaces  at.  different  sites,  and  therefore  the  results 
of  Figure  4-5  should  be  interpreted  as  presence  of  UCDW.  not  a  contiguous  warm 
intrusion.  Second,  the  UCDW  signal  is  generally  weak  when  found  on  the  shelf,  with 
temperatures  slightly  over  1.5 °C  in  thin  layers,  with  a  few  exceptions  where  the  CTD 
cast  shows  a  very  strong  0-S  signature  (Figure  4-6).  Third,  in  the  five  days  or  so 
that  it  takes  to  conduct  the  CTD  survey  over  the  northern  part  of  the  study  region 
where  the  warm  water  was  found,  one  or  two  distinct  events  lasting  1  day  or  so  each 
were  registered  at  the  A2  mooring.  This  is  inconsistent  with  the  idea  of  a  large-scale 
intrusion  moving  onto  the  shelf,  because  typical  velocities  on  the  shelf  are  of  the 
order  of  5  to  10  cm/s,  which  results  in  an  advective  time-scale  of  8  to  20  days  for  an 
intrusion  moving  across  the  shelf  as  far  as  is  typically  observed  in  Figure  4-5. 

These  results  suggest  that  the  UCDW  is  intruding  as  frequent,  short  events,  and 
not  as  large-scale  intrusions.  These  events  might  be  instead  eddies  moving  onshore, 
and  leaving,  as  they  mix  away  on  the  shelf,  the  broad  regions  of  weak  UCDW  signature 
which  are  typically  observed  in  the  mid-shelf.  In  the  next,  section,  I  explore  this 
hypothesis  using  a  simple  model  applied  to  the  A2  mooring  data. 


4.5.1  Eddies  Advected  Pass  a  Mooring 

Let.  us  assume  the  UCDW  events  are  eddies  carried  by  the  shelf  circulation  across  the 
shelf  via  Marguerite  Trough.  The  flow  as  observed  by  the  mooring  would  be  given  by 
(i.e.  Lilly  and  Rhines,  2002): 

u  =  U  +  tic  (4.1) 
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Figure  4-13:  Average  along-isobath  (left  panel)  and  cross-isobath  (right,  panel)  veloc¬ 
ity  for  each  of  the  40  UCDW  intrusions  recorded  at  the  A2  mooring  site.  The  red 
lines  are  the  ensemble  averages  for  the  along-  and  cross-slope  velocity  components, 
and  the  green  lines  the  average  velocity  for  times  when  no  intrusions  were  recorded. 


where  u  is  the  observed  flow,  U  is  the  background  flow,  which  will  be  assumed  to  be 
spatially  uniform  but  slowly  evolving  in  time,  and  ue  is  the  flow  of  the  eddv  being 
adverted  by  the  background  flow  pass  the  mooring.  A  sensible  option  for  the  eddy 
flow  is  given  by  a  Rankine  vortex,  where  the  core  is  in  solid  body  rotation,  and  outside 
of  which  the  radial  velocity  decays  with  distance  from  the  center  of  the  eddy.  In  polar 
coordinates  the  azimuthal  velocity  is  given  by 


v(r) 


VrR~l  ,  r  <  R 
Vr~xR  ,  r  >  R 


(4.2) 
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where  r  is  the  distance  from  the  center  of  the  eddy,  V  is  the  maximum  azimuthal 
velocity  (at  r  =  R),  and  R  is  the  radius  of  the  eddy.  V  is  taken  here  to  be  positive 
for  an  anticyclonic  eddy.  The  velocity  field  associated  with  an  anticyclonic  eddy  is 


X-distance  [km] 


Time  [days]  Time  [days] 


Figure  4-14:  A  Rankine  vortex.  The  left  panel  shows  the  horizontal  flow  associated 
with  an  eddy  of  R  —  8  km  and  V  =  2  cm/s.  The  contours  are  the  pressure  field.  The 
four  panels  on  the  right  show  the  zonal  (solid)  and  meridional  (dashed)  velocity  as 
would  be  recorded  by  a  fixed  mooring  as  the  eddy  is  advected  along  the  horizontal 
lines  shown  in  the  left  panel  by  a  constant  flow  of  10  cm/s  from  left  to  right. 

shown  in  Figure  4-14.  In  t  his  example,  the  eddy  has  a  radius  of  8  km  and  a  maximum 
azimuthal  velocity  V  of  2  cm/s.  The  four  panels  on  the  right  show  the  velocity 
components  as  they  would  be  seen  by  a  mooring  at  different  y-locations  represented 
by  the  horizontal  lines  on  the  left  panel,  assuming  the  eddy  is  advected  by  a  constant 
10  cm/s  zonal  (from  left  to  right)  flow  .  The  zonal  eddy  velocity  is  symmetric  and  it 
goes  from  negative  values  south  of  the  center  of  the  eddy  to  positive  north  of  it.  The 
meridional  velocity  is  antisymmetric  and  always  positive  first  when  measured  at  our 
hypothetical  moorings  and  it  decreases  (increases)  linearly  in  time  (space)  inside  the 
core  of  the  eddy.  The  opposite  is  true  in  the  case  of  cyclonic  eddies. 

If  the  UCDW  events  are  eddies  being  advected  past  the  moorings,  one  would 
expect  the  "cross-stream’'  velocity,  the  component  perpendicular  to  the  direction  of 
advection,  to  have  the  characteristic  signature  shown  in  Figure  4-14:  early  positive 
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(negative)  velocities  for  an  anticyclone  (cyclone).  However,  the  example  above  has 
the  eddy  adver  ted  by  a  constant  zonal  velocity.  This  is  bound  to  be  poor  assumption 
in  most  realistic  flows,  and  therefore  one  must  estimate  the  velocity  and  direction  of 
the  adverting  flow  U  as  it  evolves.  The  eddy  flow  then  becomes 

iTe  =  (u-Utf)e-l<t,{t)  (4.3) 

where  Uif  is  the  low-pass  velocity  (5-day  half-power  lanczos  filter)  with  direction  q. 
Ue  =  lie  +  ive  is  therefore  rotated  so  that  ve  is  the  cross-stream  component  equivalent 
to  the  meridional  velocity  in  Figure  4-14.  Any  small  linear  trends  remaining  in  t>  are 
removed. 

There  is  an  important  caveat  to  keep  in  mind  as  the  results  are  presented  in  the 
next  sections:  a  single  mooring  recording  passing  eddies  would  necessarily  tend  to 
underestimate  the  intensity  and  size  of  the  eddies  as  chances  are  that,  the  "sampling” 
is  off-center.  Therefore,  the  average  properties  over  many  events  are  biased  and  do  not 
represent  the  true  properties  of  the  eddy  “population”,  even  under  the  assumption 
that  A2  is  recording  all  the  eddies  moving  through  Marguerite  Trough. 

4.5.2  Eddy-like  Events  in  the  A2  Records 

In  order  to  test  the  simple  model  outlined  above,  the  events  must  be  relatively  isolated. 
As  was  described  above,  the  A2  records  contain  40  warm  events;  inspection  of  the 
records  reveal  that  most  of  the  longer  (>3  days)  events  appear  to  be  closely  spaced 
short  events.  Several  temperature  peaks  larger  than  1.5  °C  are  evident  at  each  depth 
in  those  cases,  and  individual  events  are  difficult  to  separate. 

For  the  analysis  below.  24  warm  events  which  appear  to  be  reasonably  isolated 
in  time  were  selected.  Of  those,  19  (79%)  had  a  velocity  signature  consistent  with 
an  anticyclonic  eddy,  four  (17%)  with  a  cyclonic  eddy  and  one  was  not  consistent 
with  either.  Examples  of  a  cyclonic  and  anticyclonic  event  are  shown  in  Figure  4-15. 
Both  of  them  have  a  strong  temperature  anomaly  at  250  m  and  a  velocity  structure 
spanning  the  300-m  depth  range  of  the  velocity  records.  The  events  have  scales  of  a 
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Figure  4-15:  Potential  temperature  (color  contours),  and  cross-stream  (t?e)  velocity 
records  (black  contours)  for  a  cyclonic  (top  panel)  and  an  anticyclonic  (bottom  panel) 
event  recorded  at  A2.  Positive  (negative)  cross-stream  velocities  are  shown  as  solid 
(dotted)  contours. 


day  or  so  and  average  cross-stream  velocities  of  1  to  5  cm/s. 

Figure  4-1G  shows  the  potential  temperature  and  velocity  records  for  the  23  events 
with  a  consistent  "eddy'  signature.  In  order  to  compare  events  of  different  duration, 
the  time  coordinate  is  normalized  using  the  250-m  records,  where  the  temperature 
signal  is  strongest,  as 


=  tm)/Tvmax  (4.4) 

where  t  is  the  time,  tm  is  the  center  time  between  the  two  absolute  maximum  ve 
during  an  event  (which  corresponds  to  the  limits  of  the  eddy  chord),  and  Tvmax  is  the 
time  period  between  those  velocity  maximums. 

On  average,  the  anticyclonic  events  have  a  maximum  temperature  signal  of 
1.6±0.09°C  (±1  standard  deviation),  with  a  weak,  but  consistent,  circulation  with 
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maximum  velocities  as  large  as  4  cm/s  at  250  m  or  so  but  only  slightly  over  0.9T0.7  cm/s 
on  average.  The  eddies  have,  on  average,  a  significant  velocity  structure  at  250  and 
400  m,  with  stronger  velocities  at  mid-depth.  The  velocity  signal  at  100  in  is  not  very 
consistent.  The  cyclonic  events  have  a  circulation  and  temperature  signal  of  com¬ 
parable  magnitude,  although  because  they  are  rarer  their  statistics  are  less  certain. 
Also,  notice  that  the  velocity  does  not  tend  to  diminish  at  the  beginning  and  end  of 
the  period  sampled  (Figure  4-16),  which  is  inconsistent  with  the  simple  ranking  eddy 
model  (Figure  4-14). 


4.6  Events  Intensity  and  Depth  Scale 

A  measure  of  the  intensity  of  the  eddies  can  be  obtained  from: 

(A0(z)>  =  (Mz)  -  MO)  (4.5) 

where  0,n  is  the  maximum  potential  temperature  during  an  event  at  depth  z  and  0,/  is 
a  profile  of  reference  potential  temperature  calculated  from  averages  of  the  time  series 
in  times  when  no  intrusions  were  detected.  The  ()  indicate  an  ensemble  average  of 
all  the  events.  The  same  procedure  is  used  to  calculate  equivalent  profiles  for  salinity 
((AS))  and  potential  density  ((Arr#)). 

Figure  4-17  shows  the  averaged  vertical  profiles  of  (A#),  (AS)  and  (A^)  for 
the  same  events  plotted  in  Figure  4-16.  On  average,  the  anticyclonic  events  have  a 
temperature  signature  at  mid-depth  with  intensities  of  0.26±0.03°C  and  salinity  of 
0.03±0.003  at  250  m.  The  slightly  lower  0  associated  with  the  events  in  the  bot¬ 
tom  instrument,  though  small,  is  consistent  with  the  hydrographic  structure  offshore, 
which  has  the  colder  LCDVV  located  below  UCDW. 

Despite  having  a  significant  temperature  and  salinity  signature,  these  warm  and 
salty  events  have,  on  average,  a  significant,  density  signal  only  at  100  m,  where  the 
difference  is  0.02±0.01.  This  signature  is  consistent  with  the  anticyclonic  velocity 
signature  shown  in  Figure  4-16.  The  weak  density  signal  is  clearly  shown  in  the  0-S 


125 


2 


1.8 


E 

o 


I 


0Lt 

I  / 


\  / 

■\ . *■■■ 

'% ,  y 


-2 


Figure  4-16:  Potential  temperature  (top  panel)  and  cross-stream  (ve)  velocity  records 
at  100.  250  and  400  m  (bottom  panel)  for  warm  eddy-like  events  at  A2.  The  left 
panels  show  the  anticyclonie  events  and  the  right  panel  the  cyclonic  events.  The  gray 
dots  indicate  individual  records,  and  the  red  dashed  lines  are  ensemble  averaged  for 
all  the  events.  The  time  is  normalized  using  equation  4.4. 

diagram  for  the  events  shown  in  Figure  4-18.  The  average  density  at  250  m  is  roughly 
27.75  regardless  of  the  presence  of  the  intrusions. 

For  the  cyclonic  events  (bottom  panels  in  Figure  4-17),  the  hydrographic  structure 
is  more  ambiguous.  They  also  show  a  stronger  temperature  signal  at  mid-depth,  with 
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the  salinity  and  temperature  anomalies  nearly  compensating  there.  And  although 
the  density  signature  is.  on  average,  negative  at  100  in.  the  errors  are  large  (there  are 
only  4  events  being  averaged  here),  and  therefore  no  firm  conclusions  can  be  reached 
about  their  vertical  structure. 

The  depth  scale  of  the  events  is  harder  to  study  using  the  moored  array  data.  In 
both  A2  and  A3,  there  are  only  a  few  instruments  in  the  300-m  depth  range  where 
UCDW  was  detected.  However,  the  vast  majority  of  the  events  were  recorded  at 
250  m  in  both  moorings  (Figure  4-19),  with  only  a  handful  being  detected  only  at 


200  m  (2  at  A2  and  none  at  A3)  or  at  400  m  (2  at  A2  and  1  at  A3).  Also,  12%  of  the 


events  (5)  at  A2  were  detected  at  both  200  m  and  400  m.  For  the  events  detected  by 
at  least  two  instruments,  the  results  indicate  that  20%  of  the  events  have  a  minimum 
thickness  of  150  m-200  m  and  a  maximum  of  «400  in,  as  the  temperature  records  at 
147  in  and  554  in  at  A2  did  not  register  UCDW. 


4.7  Horizontal  Scale  of  the  Eddies 

An  estimate  of  the  horizontal  scale  of  the  eddies  can  be  obtained  from 


(4.6) 


Horizontal  scales  for  the  intrusions  are  roughly  4.2±2.5  km  and  4.3±3.5  km  for  the 
anticyelonic  and  cyclonic  events,  respectively.  A  rough  estimate  of  the  internal  radius 
of  deformation  for  the  deep  layer  using  N  =  1.3  •  10-  5  s-1,  H  =  400  m  and  /  = 
1.35  •  10-J  s_l  results  in  A  =  3.9  km.  Notice  that  because  the  “sampling''  of  the 
eddies  by  the  moorings  is  biased  towards  underestimating  their  size  and  intensity, 
this  estimate  appears  to  be  rather  reasonable  and  dynamically  consistent  with  the 
idea  that  the  UCDW  events  are  indeed  eddies  recorded  by  the  mooring. 


(a  e>  [°C]  <a  s>  <A  °6>  [k9  m_3l 


Figure  4-17:  Vertical  profiles  of  (A 0),  (AS)  and  (Aae)  for  the  A2  mooring.  The  top 
panels  correspond  to  19  events  identified  as  anticyclones,  and  the  bottom  panel  the  4 
events  identified  as  cyclones.  The  gray  shading  indicates  the  95%  confidence  interval 
for  the  ensemble  averages. 
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Figure  4-18:  0-S  diagram  of  average  A2  mooring  records  during  times  in  the  absence 
of  UCDW  intrusions  (blue),  when  UCDW  was  detected  and  the  velocity  signal  was 
anticyclonie  (red)  and  when  the  velocity  signal  was  cyclonic  (green).  Also  indicated 
are  the  instrument  depths  (in  meters). 


4.8  Heat  Content 


As  mentioned  previously,  estimates  from  hydrographic  surveys  and  modeling  studies 
for  the  region  suggest  that  a  vertical  heat  flux  of  1.5-5  W/m2  is  necessary  to  dost' 
the  heat  budget  on  the  shelf  in  the  SO  GLOBEC  study  region  (Klinck,  1998:  Smith 
and  Klinck,  2002).  Making  an  estimation  of  the  heat  flux  due  to  the  warm  UCDW 
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Figure  4-19:  Depth  span  of  intrusions  at  A2  (red)  and  A3  (blue).  Events  are  grouped 
by  the  instruments  at  which  single  events  were  recorded. 

intrusions  studied  here  provides  at  best  an  order  of  magnitude  estimate  for  only  one 
of  the  locations  for  cross-shelf  exchange  of  heat,  but  it  is  important  even  to  have  a 
rough  idea  of  the  contribution  of  these  processes  to  the  overall  heat  budget. 

To  make  this  calculation,  I  wrill  assume  that  each  warm  eddy  is  adverted  onto  the 
shelf  and  is  mixed  away,  and  that  when  it  moves  onto  the  shelf,  an  equivalent  volume 
of  water  with  some  reference  temperature  characteristic  is  advected  away  from  the 
region  of  interest  (Lilly  and  Rhines,  2002).  Consider  first  the  heat  content  of  a  warm 
event: 


(4.7) 


where  pa  is  a  reference  density,  cp  is  the  specific  heat  capacity,  h  is  the  bottom  depth 
and  6  is  the  potential  temperature.  On  average,  the  heat  gained  with  respect  to 
a  reference  temperature  profile  is  (A  H)  =  (Hev  —  Hrf).  Using  the  vertical  profiles 
of  potential  temperature  for  the  times  with  and  without  events  (the  profile  of  the 
difference  is  shown  in  Figure  4-17),  (A H)  —  1.72±0.35-108  J/in-2.  For  this  calcula¬ 
tion,  both  profiles  are  extrapolated  to  the  surface  and  the  bottom  using  the  c  losest 
available  value. 

To  estimate  the  contribution  over  the  entire  shelf  from  this  single  average  event,  a 
ratio  of  the  area  of  a  single  “eddy”  to  the  entire  shelf  is  calculated.  For  the  event,  we 
consider  a  circle  of  radius  r  =  4.2  ±  2.5  km,  which  is  the  lateral  scale  for  the  eddies 
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we  calculated  in  the  previous  section.  For  the  area  of  the  entire  shelf,  we  consider  the 
study  area  of  the  SO  GLOBEC  program  (the  same  used  for  the  published  estimates 
of  the  vertical  heat  flux),  that  is,  from  the  northern  tip  of  Adelaide  Island  to  roughly 
the  southern  tip  of  Charcot  Island,  an  area  of  around  1-105  km2.  The  ratio  of  both 
areas  is  R  =  0.001  ±  0.0005,  which  results  in  an  average  for  the  heat  gain  on  the  shelf 
of  (NH)  =  1.72  ±  0.93  •  105J/m2. 

The  equivalent  vertical  heat  flux  associated  with  these  events  is 

Qev  =  (AH)E  (4.8) 

where  E  is  the  number  of  events  per  year.  Assuming  there  are  46  events  a  year  (40  in 
10.5  months),  Qev  =0.25±0.15  W/m2.  With  all  uncertainties  considered,  the  recorded 
events  could  contribute  between  5%  to  15%  of  the  required  heat  flux  to  balance  the 
heat  budget. 

The  above  number  is  probably  a  lower  bound  for  a  number  of  reasons.  First,  it 
is  likely  that  there  is  more  than  one  location  between  Adelaide  Island  and  Charcot 
Island  where  UCDW  intrusions  occur.  Second,  Marguerite  Trough  is  some  80-km 
wide  where  A2  is  located,  so  eddies  with  an  horizontal  scale  of  5-10  km  can  be  easily 
missed  if  the  advective  flow  scales  with  the  bathymetric  scales.  Third,  because  of 
the  dependence  of  R  on  the  square  of  r,  the  sensitivity  to  the  horizontal  scale  of  the 
intrusions  is  rather  critical.  For  example,  an  “average  eddy”  with  a  horizontal  scale 
of  2A  would  contribute  Qev  =  0.50  W/m2.  Finally,  events  counted  as  a  single  eddy 
here  tend  to  group  together  several  eddies,  which  has  a  direct  effect  on  E  above. 

4.9  Discussion  and  Conclusions 

In  this  chapter,  I  have  characterized  intrusions  of  warm,  salty  water  onto  the  wAP 
shelf  which  are  thought  to  strongly  influence  the  heat  budget.  A  number  of  conclusions 
can  be  drawn  from  this  work. 

First,  the  analysis  presented  here  suggests  that  Marguerite  Trough,  and  in  partic- 
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ular  its  eastern  side,  is  a  path  for  warm  water  to  move  across  the  shelf.  Almost  four 
times  as  many  warm-water  events  were  recorded  on  the  eastern  side  of  Marguerite 
Trough  (at  A2)  compared  to  the  west  side  (at  A3). 

Second,  the  events  are  much  smaller  in  scale  and  more  frequent  than  had  been 
thought  previously  on  the  basis  of  large-scale  hydrographic  surveys.  The  sampling 
bias  introduced  by  the  20-km  or  so  CTD  cast  separation  typical  of  the  SO  GLOBEC 
and  other  seasonal  surveys  in  this  area  suggested  that  intrusions  had  large  horizon¬ 
tal  scales  (significantly  larger  than  the  internal  radius  of  deformation)  and  with  a 
frequency  of  six  or  so  per  year.  The  results  presented  here  suggest  otherwise.  The 
UCDW  intrusions  have  a  thickness  of  a  few  hundred  meters,  and  they  occur  roughly 
four  times  a  month  on  average.  A  simple  model  of  an  isolated  eddy  being  advected 
pass  a  mooring  was  shown  to  be  somewhat  consistent  with  the  observations  at  A2, 
and  results  in  eddies  with  an  horizontal  scale  of  4  to  5  km  and  a  associated  currents  of 
1  to  2  cm/s.  A  caveat  relating  to  the  band-pass  filtering  of  the  velocity  data  needs  to 
be  taken  into  account  when  considering  the  results  presented  here.  Figure  4-20  shows 
a  synthetic  velocity  time-series  of  a  Rankine  eddy  (see  Section  4.5.1)  being  advected 
pass  a  mooring.  The  eddy  has  a  velocity  scale  of  2  cm/s  and  the  event  lasts  for  a 
few  days,  which  is  of  the  same  order  as  the  events  recorded  at  A2.  The  solid  lines  in 
Figure  4-20  are  the  original  modelled  time-series,  and  the  dashed  lines  are  the  mod¬ 
elled  velocity  records  band-pass  filtered  using  the  same  method  used  for  the  moored 
observations.  The  filtering  appears  to  affect  the  original  signal,  although  the  basic 
scales  of  the  eddy  as  seen  from  the  cross-stream  velocity  remain  roughly  the  same  in 
the  example.  However,  it  must  be  noted  that  in  the  case  of  the  synthetic  time-series 
we  know  what  the  along-stream  and  cross-stream  directions  are.  In  the  case  of  the 
moored  observations,  this  direction  was  determined  by  rotating  the  velocity  in  the 
direction  of  the  low-pass  filtered  velocity,  and  thus  the  effect  of  the  filtering  might 
be  more  pronounced.  Other  methods  for  isolating  the  eddies,  and  observations  which 
resolve  the  spatial  scales  of  the  intrusions  might  be  required  to  resolve  the  question 
of  their  circulation. 

Although  this  work  does  not  attempt  to  answer  the  question  of  the  formation 
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process  of  the  intrusions,  it  does  help  to  provide  some  constraints  for  the  spatial 
and  temporal  scales  involved.  This  has  straightforward  implications  for  the  wAP 
ecosystem.  Prezelin  et  al.  (2000)  noted  that  diatom  populations  on  the  mid-  and 
outer-shelf  depend  on  the  renewal  of  the  silica-rich  UCDW.  The  work  presented  here 
suggests  that  this  renewal  is  frequent  and  more  or  less  continuous  throughout  the 
year,  which  has  important  implications  for  the  continuous  success  of  the  biological 
community. 

Finally,  the  contribution  of  the  warm-water  events  considered  here  to  the  heat 
budget  of  the  shelf  has  been  estimated.  The  results  suggest  that  Marguerite  Trough 
could  contribute  at  least  5%  to  15%  of  the  heat  flux  required  to  close  the  heat  budget. 
However,  it  must  be  noted  that  the  surface  net  ocean  heat  loss  to  the  atmosphere  on 
the  shelf  is  a  poorly  constrained  quantity,  that  the  rough  estimation  of  the  contribu¬ 
tion  of  Marguerite  Trough  is  based  on  observations  from  a  single  mooring,  and  that 
the  choice  of  a  reference  temperature  profile,  though  reasonable,  is  not  necessarily 
the  correct  one.  All  of  these  caveats  must  be  addressed  in  future  studies  of  the  heat 
and  nutrients  budgets  on  the  wAP  shelf. 

Overall,  these  results  suggest  that  the  circulation  and  property  budgets  of  the  wAP 
shelf,  are  heavily  influenced  by  the  ACC  and  the  complicated  bathymetry,  are  fully 
three-dimensional,  and  that  many  of  the  important  processes,  including  the  UCDW 
intrusions,  have  time  and  spacial  scales  small  enough  that  would  require  process  or 
modeling  studies  resolving  horizontal  distances  of  1  to  2-km  and  time  scales  of  a  day 
or  so. 
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Figure  4-20:  Effect  of  the  band-pass  filtering  of  the  velocity  records.  The  top  (bottom) 
panel  shows  the  along-stream  (cross-stream)  velocity  records  from  a  synthetic  time 
series  using  the  Rankine  eddy  model  described  in  Section  4.5.1.  The  solid  lines  are  the 
original  time-series  of  an  eddy  with  a  maximum  velocity  of  2  cm/s,  and  the  dashed 
lines  are  the  band- pass  filtered  time  series  (half-power  of  20-h  and  5  days). 
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Chapter  5 


A  First  Description  of  the 
Antarctic  Peninsula  Coastal 
Current 

5.1  Introduction 

The  west  Antarctic  Peninsula  (wAP)  continental  shelf  hosts  a  rich  biological  com¬ 
munity  centered  around  Antarctic  krill  ( Euphasia  superba).  The  region  has  been 
hypothesized  to  act  as  a  source  for  krill  populations  elsewhere  in  the  Southern  Ocean 
(Murphy  et  alM  1998;  Hofmann  and  Klinck,  1998;  Fach  and  Klinck,  2006),  as  well 
as  to  sustain  a  number  of  higher  trophic  levels  in  the  peninsula’s  food  web  (Costa 
and  Crocker,  1996;  Fraser  and  TYivelpiece,  1996).  Additionally,  multiple  studies  of 
the  regional  climate  have  found  that  the  surface  air  temperature  of  the  peninsula  is 
rising  at  one  of  the  highest  rates  in  the  world  and  causing  rapid  changes  in  the  region 
(Vaughan  and  Doake,  1996;  Cook  et  al.,  2005).  Ocean  dynamics  of  the  wAP  shelf 
are  expected  to  play  an  important  role  in  these  processes,  but  they  remain  poorly 
understood.  The  Southern  Ocean  Global  Ecosystem  Dynamics  (SO  GLOBEC)  field 
program,  conducted  between  2001  and  2003,  was  aimed  at  understanding  the  domi¬ 
nant  physical  and  biological  processes  in  Marguerite  Bay  and  the  adjacent  wAP  shelf 
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region  (Hofmann  et  al.,  2002,  2004). 


Prior  to  SO  GLOBEC,  the  shelf  circulation  had  been  addressed  by  only  a  handful 
of  studies,  mostly  involving  observations  of  the  hydrographic  fields  during  austral 
summer.  In  a  review  paper,  Hofmann  and  Klinck  (1998)  describe  the  circulation 
in  the  SO  GLOBEC  study  area  (Figure  5-1)  as  a  slow,  cyclonic  gyre  composed  of 
a  north-eastward  branch  flowing  along  the  shelf  break  and  a  southward  flow  along 
the  inner  shelf.  The  shelf-break  circulation  is  dominated  by  the  large-scale  Antarctic 
Circumpolar  Current  (ACC),  while  the  near-shore  southward  flow  was  thought  to  be 
part  of  a  shelf-scale  gyre.  However,  two  recent  studies  hint  at  the  presence  of  an 
independent  southward  inner  shelf  circulation.  Tracks  of  drifters  released  during  the 
summer  of  2001  and  2002  (as  part  of  SO  GLOBEC)  show  a  southward  current  flowing 
along  Adelaide  Island  and  into  Marguerite  Bay,  and  then  exiting  on  the  south  side  of 
the  Bay  along  Alexander  Island  (Beardsley  et  al.,  2004).  The  SO  GLOBEC  regional 
surface  dynamic  topography  for  the  fall  (April-June)  of  2001  is  consistent  with  those 
results,  and  the  low  salinity  water  found  on  the  inner  shelf  suggests  the  flow  is  a 
buoyancy-forced  coastal  current  (Klinck  et  al.,  2004). 


Here  we  use  the  SO  GLOBEC  data  set,  including  broad-scale  hydrographic  data 
and  mooring  observations,  to  provide  a  first  description  of  this  nearshore  coastal 
current.  We  will  refer  to  it  as  the  Antarctic  Peninsula  Coastal  Current  (APCC).  We 
discuss  the  vertical  structure,  transport,  along- shelf  variability,  and  seasonal  evolution 
of  the  APCC.  We  also  explore  the  possible  forcing  mechanisms  and  the  implications 
of  our  results  for  the  distribution  of  freshwater  on  the  shelf. 


The  paper  is  organized  as  follows:  we  summarize  the  available  data  and  process¬ 
ing  (Section  5.2),  followed  by  brief  descriptions  of  the  study  area  and  background 
oceanographic  and  meteorological  conditions  (Section  5.3),  results  (Section  5.4),  and 
discussion  (Section  5.5).  A  final  summary  is  presented  in  Section  5.6. 
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5.2  Data  Sources 


During  2001,  two  hydrographic  cruises  were  conducted  off  the  west  coast  of  the 
Antarctic  Peninsula  aboard  the  RV/IB  N.  B.  Palmer  as  part  of  the  SO  GLOBEC 
field  program  (Table  5.1).  These  broad-scale  cruises  attempted  to  cover  the  entire 
study  region,  and  provided  CTD  (conductivity  for  salinity,  temperature  and  pressure 
for  depth)  profile  observations  on  a  grid  with  a  nominal  cross-shelf  resolution  of  20 
km  (Figure  5-1).  The  fall  survey  (April  24th  to  June  5th)  provided  a  larger  spatial 
coverage  than  the  winter  survey  (July  24th  to  August  31st),  in  which  weather  condi¬ 
tions  and  sea- ice  coverage  impeded  navigation.  Sea  surface  temperature  (SST)  and 
sea  surface  salinity  (SSS)  as  well  as  wind  speed  and  direction  data  were  also  collected 
throughout  the  cruises  using  the  RV/IB  N.  B.  Palmer  continuous  sampling  systems. 
Details  of  the  processing  of  the  data  set  can  be  found  in  Wiebe  et  al.  (2001)  and 
Klinck  et  al.  (2004). 

As  part  of  these  broad-scale  cruises,  conductivity,  temperature  and  pressure  data 
from  depths  of  roughly  15  m  to  250  m  were  also  collected  using  the  Bio-Optical 
Multi- frequency  Acoustical  and  Physical  Environmental  Recorder  (BIOMAPER-II), 
a  towed  undulating  multi-sensor  system  (Wiebe  et  al.,  2002).  A  subset  of  these  data 
collected  on  June  1st,  2001  (during  the  fall  cruise)  were  used  to  construct  a  high- 
resolution,  cross-shelf  hydrographic  section  off  Adelaide  Island  (Figure  5-1,  Section  I) 
by  merging  a  near-shore  CTD  station  (Station  101),  the  BIOMAPER  data,  and  the 
underway  SST  and  SSS  data. 

Velocity  data  were  obtained  during  the  cruises  with  a  shipboard-mounted  RD  In¬ 
struments  (RDI)  150-kHz  narrow-band  Acoustic  Doppler  Current  Profiler  (ADCP) 
sampling  from  31  to  431  m.  The  ADCP  was  configured  with  fifty  eight-meter  bins, 
an  eight-meter  pulse  length,  and  five-minute  ensemble  averages.  During  most  of  the 
cruise,  the  ADCP  operated  in  bottom-tracking  mode.  The  data  processing  was  ini¬ 
tially  done  during  the  cruise  using  the  Common  Oceanographic  Data  Access  System 
(CODAS).  Final  processing  and  quality  control  were  done  by  by  J.  Hummon  (Univer¬ 
sity  of  Hawaii)  and  S.  Howard  (Earth  Systems  Research)  (Klinck  et  al.,  2004).  ADCP 
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data  return  during  the  winter  cruise  was  greatly  reduced  while  the  N.  B.  Palmer  was 
moving  through  sea-ice  and  data  from  those  times  are  not  considered  here.  For  the 
fall  cruises,  the  small  (<0.02  m  s-1)  barotropic  tidal  currents  were  removed  from  the 
data  set  using  the  Antarctic  Peninsula  tide  model,  version  04.01  (AntPen04.01),  a 
high-resolution  barotropic  model  based  on  the  linearized  shallow  water  equations  (L. 
Padman,  personal  communication  -  See  Chapter  3). 

Velocity,  salinity,  and  temperature  observations  were  also  collected  at  the  A1 
mooring,  located  16.3  km  off  the  west  coast  of  Adelaide  Island  at  a  water  depth  of 
509  m  (Figure  5-2).  The  Al  mooring  was  deployed  and  recovered  by  the  R/V  Lau¬ 
rence  M.  Gould  and  collected  data  from  March  2001  to  February  2002  (Table  5.1). 
The  mooring  was  equipped  with  SeaBird  Electronics  (SBE)  temperature  recorders 
at  134,  209  and  259  m  depth,  SBE  MicroCATs  measuring  both  temperature  and 
salinity  at  159  and  309  m,  and  a  conductivity  and  temperature  sensor  set  attached 
to  a  Vector-Averaging  Current  Meter  (VACM)  at  459  m.  A  SBE  SeaGauge  pres¬ 
sure  sensor  at  502  m  recorded  pressure  and  additional  temperature  and  conductivity 
data.  Velocity  measurements  included  an  upward-looking,  RDI  300-kHz  broadband 
WorkHorse  ADCP  located  at  169  m  depth  and  configured  with  sixty  two-meter  bins 
(ranging  from  47  to  165  m)  and  half-hour  ensemble  averages.  Four  bins  (107,  109, 
159  and  161  m)  were  replaced  by  linear  interpolation  due  to  contamination  by  the 
acoustic  reflection  from  other  instruments  farther  up  on  the  mooring  line.  Finally, 
two  VACMs  also  recorded  currents  averaged  over  fifteen  minutes  intervals  at  309  and 
459  m.  For  details  of  the  moored  array  design  and  data,  see  Moffat  et  al.  (2005). 

Ice  concentration  data  for  the  study  region  were  obtained  from  the  Special  Sensor 
Microwave/Imager  (SSM/I)  data  set  derived  from  the  bootstrap  passive  microwave 
algorithm  and  distributed  by  the  National  Snow  and  Ice  Data  Center  (NSIDC).  The 
sea-ice  data  grid  has  a  cell  size  of  25  x  25  km,  and  data  are  provided  as  daily  averages 
(Comiso,  1999,  updated  2005).  Ice  thickness  data  were  also  obtained  at  a  series  of 
ice  stations  during  the  winter  cruise  (Perovich  et  al.,  2004),  and  throughout  the  SO 
GLOBEC  mooring  deployment  (March  2001  to  February  2002)  at  the  A2  mooring 
(Figure  5-1  and  Table  5.1),  which  was  equipped  with  an  Ice  Profiling  Sonar  (IPS) 
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deployed  at  the  top  of  the  mooring  line.  The  IPS  measured  the  acoustic  range  to  the 
surface  (the  sea  surface  or  ice  bottom  when  present)  every  two  seconds,  which  was 
then  converted  to  distance  (and  thus,  ice  draft)  using  in  situ  sound  speed  and  density 
calculated  from  a  nearby  sensor  (Moffat  et  ah,  2005).  All  mooring  time  series  data 
shown  here  were  low-pass  filtered  using  a  33-hr  half-amplitude  filter  (Flagg  et  al., 
1976;  Limeburner  et  ah,  1983)  to  remove  tidal  and  inertial  band  variability.  All 
salinities  are  given  in  the  practical  salinity  scale. 


5.3  Background 

The  SO  GLOBEC  study  area  extends  nearly  500  km  from  the  northern  tip  of  Adelaide 
Island  to  Latady  Island  (Figure  5-1).  The  wAP  shelf  is  relatively  wide  (~130  km), 
deep  and  characterized  by  complex  bathymetry  that  includes  shallow  plateaus  and 
deep  depressions,  with  mean  and  maximum  depths  of  «450  m  and  1600  m,  respec¬ 
tively  (Bolmer  et  ah,  2004).  This  complexity  is  also  evident  inside  Marguerite  Bay, 
where  small  islands,  subsurface  reefs  and  troughs  are  ubiquitous. 

The  general  features  of  the  shelf  hydrographic  structure  have  been  described  in 
previous  studies  (Hofmann  and  Klinck,  1998;  Smith  et  ah,  1999;  Klinck  et  ah,  2004), 
and  are  illustrated  by  the  potential  temperature-salinity  (0-S)  diagram  and  vertical 
profiles  in  Figure  5-3.  The  data  were  collected  using  a  shipboard  CTD  at  the  A2 
mooring  site  during  the  2001  mooring  deployment,  the  fall  and  winter  cruises,  and  the 
2002  mooring  recovery  (Table  5.1).  The  surface  waters  of  the  wAP  region  are  occupied 
by  Antarctic  Surface  Water  (AASW),  a  water  mass  with  temperatures  ranging  from 
-1.8  to  1°C  and  relatively  low  salinities  of  33.0  to  33.7  (Smith  et  ah,  1999).  By  the  end 
of  the  austral  fall  (end  of  June),  strong  surface  heat  loss  leads  to  the  formation  of  a 
deep  winter  mixed  layer  and  a  sea-ice  cover  over  the  entire  shelf.  As  the  sea-ice  retreats 
and  later  during  the  ice-free  season,  the  surface  mixed  layer  is  partially  restratified 
by  surface  heating  and  freshening  by  sea-ice  melt  and  runoff.  The  remnant  of  that 
cold  winter  mixed  layer,  called  Winter  Water  (WW),  is  found  at  roughly  80-  to  100-m 
depth  during  the  ice-free  season,  with  temperatures  of  -1°C  or  less.  Subpycnocline 
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waters,  below  150  to  200  m,  are  dominated  by  a  modified  version  of  the  Circumpolar 
Deep  Water  (CDW),  a  warm  (1.5°C)  and  salty  (34.6  to  34.73)  water  mass,  a  source  of 
heat  and  nutrients  for  the  upper  layer  of  the  shelf  (Smith  et  al.,  1999;  Serebrennikova 
and  Fanning,  2004). 


5.4  Structure  and  Dynamics  of  the  APCC 

5.4.1  Spatial  Structure 

During  the  fall  cruise,  maps  of  SSS  and  SST  (Figure  5-4)  show  the  fresher  end  of 
the  AASW  near  the  coast,  where  salinities  are  typically  lower  than  33.25.  These 
fresher  waters  are  found  along  Adelaide  Island  and  particularly  in  the  region  off 
Alexander  Island  and  Charcot  Island,  where  the  surface  salinity  is  at  its  minimum. 
SST  (Figure  5-4)  is  above  zero  in  the  northern  half  of  the  study  region,  and  the  coolest 
temperatures  are  found  off  Alexander  and  Charcot  Islands,  and  along  the  shelf  break. 
The  north-south  temperature  gradient  is  likely  a  combination  of  the  shelf- wide  cooling 
during  the  four- week  cruise,  which  proceeded  from  north  to  south  (Klinck  et  al.,  2004), 
and  a  mean  gradient  reflecting  the  regional  climate  (Stammerjohn  et  al.,  2003). 

Potential  temperature-salinity  (0-S)  diagrams  (Figure  5-5)  illustrate  the  cross¬ 
shelf  hydrographic  structure  during  the  fall  cruise.  Far  offshore,  the  vertical  structure 
of  the  water  column  resembles  the  structure  characteristic  of  the  outer  shelf  and  open 
ocean  as  described  in  Section  5.3  (black  solid  lines  in  Figure  5-5).  Near  the  coast, 
however,  salinity  anomalies  of  the  order  of  0.3  to  0.7  extend  down  to  100-200  m  (red 
dashed  lines  in  Figure  5-5).  Also,  the  subsurface  temperature  minimum  that  is  the 
signature  of  the  Winter  Water  (WW)  is  much  weaker  near  the  coast  than  offshore.  At 
Section  C,  WW  is  completely  absent  and  temperature  increases  monotonically  with 
depth.  This  absence  of  the  WW  during  the  fall  is  unique  to  the  nearshore  areas  and 
Marguerite  Bay  (not  shown). 

A  detailed  hydrographic  and  velocity  section  (Figure  5-6)  taken  off  the  coast  of 
Adelaide  Island  (Section  I,  Figure  5-1)  shows  the  vertical  and  cross-shelf  structure 
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of  the  APCC.  As  suggested  by  the  fall  broad-scale  CTD  survey,  freshwater  is  con¬ 
centrated  near  the  coast.  The  front  separating  the  fresh  coastal  current  from  the 
shelf  waters  is  centered  20  km  off  the  coast.  The  salinity  anomaly  is  of  the  order 
of  0.6  at  the  surface,  and  extends  below  100  m  and  to  30  km  off  the  coast.  Unlike 
the  shelf  waters,  which  show  a  well-mixed  surface  layer  50  to  70  m  deep,  the  surface 
stratification  of  the  APCC  is  relatively  strong,  and  the  temperature  structure  differs 
markedly  from  the  temperature  structure  of  the  shelf  waters.  Winter  Water,  typically 
found  at  80  to  100  m  depth  during  the  ice-free  season  (Figure  5-3),  is  absent  from  the 
nearshore  waters,  and  is  replaced  by  a  tongue  of  warm  water  out  to  10  km  offshore  of 
the  front.  This  is  consistent  with  the  absence  of  WW  near  the  coast  described  above 
in  the  broad-scale  CTD  survey  data  (Figure  5-5). 

The  frontal  structure  of  the  salinity  field  is  typical  of  coastal  buoyant  plumes,  and 
is  reflected  in  the  velocity  field,  which  shows  a  well-defined  coastal  current  co-located 
with  the  salinity  front  (Figure  5-6).  In  order  to  obtain  the  velocity  field  above  the 
shallowest  ADCP  bin  (31  m,  Section  5.2),  we  used  the  thermal  wind  (geostrophic) 
shear  calculated  from  the  hydrographic  fields  to  extrapolate  the  ADCP  velocity  to 
the  surface.  A  comparison  of  the  mean  vertical  shear  observed  by  the  ADCP  and  the 
thermal  wind  shear  averaged  across  the  current  (Figure  5-7)  confirms  that  the  vertical 
structure  of  the  along-shelf  velocity  is  largely  determined  by  the  cross-shelf  density 
gradient  and  validates  the  use  of  geostrophic  shear  to  extrapolate  the  ADCP  data 
vertically.  Inertial  motions  and  other  wind-forced  currents  could  strongly  contaminate 
the  characterization  of  the  low-frequency  motions  we  are  addressing  here.  However, 
measurements  from  the  N.  B.  Palmer  indicate  wind  forcing  was  low  (<2  m  s_1)  while 
Section  I  was  occupied.  A  wavelet  analysis  (not  shown)  of  the  A1  mooring  velocity 
records  showed  negligible  inertial  frequency  energy  during  this  period. 

Section  I  data  show  surface  velocities  approaching  0.4  m  s-1  (Figure  5-6),  which 
decays  vertically  with  a  scale  of  approximately  80  m  and  have  a  horizontal  scale  (L) 
of  only  5  km.  The  Rossby  number  for  the  current  is  Ra  ~  V/fL,  where  /  is  the 
Coriolis  parameter  and  V  is  a  typical  velocity.  With  the  observed  scales  R0  is  «0.6, 
an  indication  that  nonlinear  effects  may  be  important  in  the  along-shelf  momentum 
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balance. 

Evidence  of  a  velocity  front  can  be  found  elsewhere  on  the  shelf.  Other  ADCP 
sections  off  Adelaide  and  Alexander  Islands  (Figure  5-8)  also  show  nearshore  flow 
moving  with  the  coast  to  the  left  (Figure  5-9).  Although  none  of  these  sections 
have  concurrent  high-resolution  hydrography  that  would  allow  examination  of  the 
dominant  frontal  features,  freshwater  fluxes  and/or  dynamical  balances,  the  structure 
and  scales  observed  in  Sections  II  and  III,  with  surface  velocities  of  0.15-0.40  m  s-1  and 
depth  scales  of  »T00  m  are  similar  to  those  from  Section  I.  Section  IV  is  different  in 
that  the  velocity  front  is  much  wider  and  weaker,  and  is  also  located  farther  offshore,  a 
displacement  perhaps  forced  by  the  shallower  bathymetry  of  the  region.  Although  all 
these  observations  suggest  a  continuous  current  along  the  shelf,  it  must  be  noted  that 
the  APCC  front  could  not  be  traced  within  Marguerite  Bay,  and  given  the  presence 
of  other  sources  of  freshwater  within  the  bay,  more  data  are  required  to  determine 
the  path  of  the  APCC  along  the  coast. 

5.4.2  Volume  and  Fresh  Water  Transports 

The  volume  transport  of  the  APCC  was  determined  by  integrating  the  ADCP  velocity 
from  Section  I  shown  in  Figure  5-6.  As  a  lower  vertical  limit  for  the  integration,  we 
used  the  34.4  isohaline.  This  choice  captures  the  frontal  structure  seen  in  the  velocity 
field,  and  is  also  a  reasonable  choice  for  the  freshwater  transport.  Offshore  of  the 
front,  the  34.4  isohaline  is  typically  located  in  the  upper  pycnocline,  above  150  m, 
and  is  deeper  than  200  m  near  the  coast  (Figure  5-6).  The  integration  from  the 
34.4  isohaline  to  the  surface  and  from  the  coast  offshore  to  24  km  results  in  a  total 
transport  of  0.32±0.13  Sv.  The  associated  freshwater  transport  of  the  APCC  in 
Section  I  is  given  by 

FWT=  [  vdA ,  (5.1) 

where  v  is  the  along-shelf  velocity,  S  is  the  observed  salinity,  Sr  is  a  reference  salinity 
(chosen  here  as  34.4),  and  A  is  the  total  area  enclosed  by  the  isohaline  and  an  offshore 
limit  of  24  km.  This  results  in  a  value  of  FWT  of  approximately  126±50  km3yr-1 
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for  Section  I  (For  convenience,  we  take  both  the  volume  and  freshwater  transports 
as  positive  quantities).  This  value  is  obtained  by  assuming  the  APCC  is  present  on 
the  shelf  for  seven  months  of  the  year.  The  work  presented  in  the  discussion  section 
provides  the  basis  for  this  assumption. 

The  uncertainty  in  the  two  transport  estimates  (about  40%)  can  be  explained  as 
follows.  Systematic  instrumental  errors  have  been  removed  during  post-processing  of 
the  data  (see  Section  5.2  and  references  therein),  and  other  instrumental  errors  are 
assumed  to  be  negligible  when  the  data  are  integrated  over  the  entire  section.  Tidal 
currents  are  small  (as  discussed  before),  and  have  been  removed  using  a  regional  tidal 
model.  The  Al  mooring  data  provide  a  good  estimate  of  the  synoptic  variability  of 
the  APCC  during  the  fall  season.  This  variability  might  be  due  to  temporal  variability 
in  the  transport,  as  well  as  cross-shelf  displacements  of  the  front  of  the  APCC.  An 
estimate  of  both  effects  can  be  calculated  from  the  ratio  of  the  mean  velocity  to  the 
standard  deviation  of  the  fall  along-shelf  velocity,  which  is  of  the  order  of  35%.  As  we 
are  unable  to  separate  temporal  variability  from  horizontal  displacements  of  the  front, 
we  assume  all  of  this  uncertainty  is  due  to  temporal  variability.  An  additional  source 
of  uncertainty  is  the  vertical  extrapolation  of  the  ADCP  data  using  the  thermal  wind 
shear.  Apart  from  this  calculation,  we  attempted  two  other  methods  of  extrapolation: 
linear  extrapolation  and  assuming  a  constant  velocity  from  the  shallowest  ADCP  bin 
to  the  surface.  The  resulting  spread  added  an  additional  5%  or  so  of  uncertainty  to 
the  transport  calculation,  giving  a  total  uncertainty  of  about  40%. 

5.4.3  Comparison  with  Theoretical  Scales 

Several  studies  have  explored  the  processes  and  dynamics  that  determine  the  relevant 
scales  of  a  geostrophically-balanced  coastal  buoyant  plume  (Yankovsky  and  Chapman, 
1997;  Fong  and  Geyer,  2001;  Lentz  and  Helfrich,  2002).  Using  the  data  from  Section 
I,  we  can  compare  some  simple  theoretical  scalings  to  those  observed  in  the  APCC. 
Among  those  scales  are  the  along-shelf  volume  transport  ( Q ),  the  total  width  of 
the  current  ( Wp ),  and  the  width  of  the  front  (VFS).  Assuming  a  buoyant  plume  in 
geostrophic  balance  over  a  motionless  deep  layer,  the  volume  transport  as  a  function 
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of  the  distance  offshore  can  be  estimated  from 


<?(*)  -5}(K-  m*)2)  •  <5-2) 

where  g'  is  the  reduced  gravity  gAp/p,  A p  is  the  density  difference  between  the  plume 
and  the  ambient  water,  g  is  the  gravitational  acceleration,  hp  is  the  depth  at  which 
the  front  touches  the  bottom,  and  h  is  the  depth  of  the  front  at  a  distance  x  offshore 
of  the  foot  of  the  front.  If  the  isopycnal  chosen  to  characterize  the  front  outcrops 
and  the  integration  is  carried  out  over  the  entire  region  enclosed  by  that  isopycnal, 
h  becomes  zero  and  the  volume  transport  Q  becomes  independent  of  x.  The  width 
of  the  front  scales  with  the  internal  radius  of  deformation  (Ws  =  \Jg'hp/ f),  and  the 
total  width  of  the  current  is  then  Wp  =  Wa  +  Wb,  where  Wb  is  the  horizontal  distance 
from  shore  to  the  isobath  hp  (Figure  5-10;  modified  from  Lentz  and  Helfrich,  2002). 

Using  the  same  34.4  isohaline  used  in  the  previous  section  provides  a  meaningful 
comparison  with  the  observed  along-shelf  transport,  as  the  velocity  below  the  chosen 
isohaline  is  reasonably  close  to  zero  (a  basic  assumption  in  the  derivation  of  Equation 
5.2).  This  isohaline  touches  the  bottom  at  hp  =220  m,  and  the  depth  h  at  the  offshore 
limit  of  the  front  is  «170  m.  The  reduced  gravity  is  computed  using  the  observed 
density  difference  across  the  density  front.  Using  either  the  density  difference  at  the 
surface  or  the  average  density  difference  above  the  34.4  isohaline  results  in  A p  ~0.45 
and  g'  ~  4.3  10-3  m  s~2.  Using  /  =  1.34 10-4s-1,  the  theoretical  volume  transport  Q 
is  ~  0.31  Sv  and  the  front  width  Ws  is  «7  km. 

The  width  of  the  front  is  estimated  from  the  observations  as  the  decay  scale  of 
the  along-shelf  velocity  maximum  on  the  surface  layer,  and  is  approximately  5±2  km. 
Both  the  width  and  the  volume  transport  calculated  from  the  data  in  Section  I  are  in 
good  agreement  with  the  theoretical  scales  (0.31  Sv  and  7  km  respectively),  suggesting 
that  the  observed  current  is  a  coastal  plume  in  geostrophic  balance  to  lowest  order. 
However,  the  observations  show  the  total  width  ( Wp )  of  the  APCC  is  larger  than 
expected  from  the  model  sketched  in  Figure  5-10.  In  the  model,  the  buoyant  plume 
adjusts  to  a  final  state  where  the  near-bottom  velocity  is  zero  and  the  front  can 
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be  found  attached  to  depth  hp.  Our  observations  suggest  that  the  front  is  farther 
offshore  than  expected,  so  that  the  total  width  is  greater  than  Wb  +  Ws.  This  can 
also  be  observed  in  the  hydrographic  data,  which  show  nearly  flat  isopleths  in  the 
upper  120  m  or  so  of  the  water  column  onshore  of  10  km.  This  is  a  possible  result 
of  the  APCC  adjusting  to  some  upstream  bathymetric  feature  or  local  forcing,  or  the 
violation  of  some  of  the  assumptions  of  the  simple  theoretical  model  used  here. 

5.4.4  Temporal  Evolution:  The  Winter  Shutdown 

As  the  winter  approaches,  heat  loss  from  the  ocean  to  the  atmosphere  increases  signif¬ 
icantly,  eventually  leading  to  the  formation  of  sea- ice  on  the  shelf.  As  a  result,  SST 
during  the  winter  months  (Figure  5-4)  is  nearly  at  the  freezing  point  everywhere. 
SSS  is  also  strongly  modified,  showing  higher  values  (partially  explained  by  brine 
rejection),  and  dramatically  reduced  horizontal  gradients  compared  to  those  for  the 
fall. 

Due  to  the  sea-ice  and  harsh  weather  conditions,  no  high-resolution  shipboard 
ADCP  or  BIOMAPER-II  survey  data  were  collected  during  the  2001  winter  cruise. 
However,  the  surface  fields  (Figure  5-4)  suggest  the  horizontal  density  gradients  are 
greatly  diminished  when  the  shelf  is  ice-covered.  This  is  confirmed  by  the  winter 
CTD  9/S  diagrams  (Figure  5-5,  right-hand  panels),  which  show  reduced  horizontal 
gradients  of  properties  during  the  winter  cruise  when  the  surface  layers  have  become 
colder  and  saltier.  Although  a  freshwater  anomaly  can  still  be  found  near  the  coast 
in  winter,  it  is  much  smaller  than  that  found  during  the  fall  cruise. 

Despite  the  lack  of  detailed  cross-shelf  sections  of  hydrography  or  velocity  during 
winter,  the  A1  mooring,  deployed  16.3  km  offshore  of  Adelaide  Island  (Figure  5-1), 
was  favorably  located  to  explore  how  the  changing  hydrographic  structure  and  sea-ice 
conditions  affect  the  evolution  of  the  coastal  current  throughout  the  year.  Low-pass 
filtered  ADCP  and  VACM  velocities  from  the  Al  mooring,  sampling  the  water  column 
from  65  to  459  m,  are  plotted  in  Figure  5-11.  During  the  fall,  the  mooring  data  reflect 
the  strongly  baroclinic  structure  evident  in  the  hydrographic  and  velocity  sections  off 
Adelaide  Island  (Section  5.4.1).  Velocities  in  the  shallowest  depth  bin  (65  m)  are 
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typically  0.15  to  0.20  m  s-1,  and  deep  velocities  remain  small  throughout  the  year. 
The  Al  site  mooring  was  located  in  a  deep  depression  with  a  downstream  (in  the 
sense  of  the  APCC  flow)  sill  depth  of  around  320  m  (the  mooring  depth  is  509  m.  - 
Figure  5-2),  which  might  explain  the  small  near-bottom  velocities  observed.  These 
near-bottom  velocities  are  not  unlike  others  found  elsewhere  on  the  shelf  in  the  SO 
GLOBEC  study  area  where  sills  and  deep  depressions  are  ubiquitous  (Moffat  et  al., 
2005). 

The  SSM/I  satellite  ice-concentration  time  series  clearly  show  the  rapid  onset  and 
slow  retreat  of  ice  cover  typical  of  the  wAP  region  (Figure  5-11;  Stammerjohn  et  al., 
2003).  As  the  shelf  transitions  from  ice- free  to  mostly  ice-covered  at  the  beginning 
of  July  2001,  the  along-shelf  current  is  abruptly  reduced  (Figure  5-11).  The  ice-draft 
time  series  from  the  A2  mooring  is  consistent  with  the  satellite  data  and  also  shows 
gradual  growth  in  thickness  of  the  sea-ice  cover,  reaching  a  draft  of  2  m  by  mid- 
September.  By  the  time  the  2001  winter  cruise  was  conducted,  there  is  little  sign  of  a 
coherent  along-shelf  flow  at  Al,  suggesting  that  the  APCC  has  completely  shut  down. 
Velocities  during  the  ice-covered  period  rarely  exceed  0.05  m  s-1,  and  northeastward 
flow  ( v  <0),  virtually  non-existent  during  the  fall,  is  frequent  during  the  winter. 

This  evolution  is  also  evident  in  the  seasonally-averaged  mooring  records  (Figure  5- 
12),  which  show  a  strong  along-shelf  velocity  profile  during  the  fall  (defined  here  until 
June,  before  the  onset  of  the  sea- ice  on  early  July),  with  velocities  in  the  shallower  bins 
of  0.20  m  s-1  or  more  and  a  depth  scale  exceeding  100  m.  The  winter  profile  shows 
significant  variability  in  the  top  90  m,  but  the  seasonal  mean  is  small  throughout  most 
of  the  water  column.  After  the  sea-ice  retreats  from  the  region  in  early  November, 
a  significant  along-shelf  flow  reappears  with  a  smaller  seasonal  mean,  reflecting  the 
reduction  of  the  velocity  and  brief  flow  reversal  in  mid-December.  Although  the 
seasonal  change  may  be  due  to  cross-shelf  displacements  of  the  front,  the  broad-scale 
CTD  surveys  during  the  winter  (discussed  above)  show  significant  reduction  of  the 
cross-shelf  density  gradients,  consistent  with  the  shutdown  of  the  coastal  current.  Ice 
breakup  and  retreat  starts  in  mid-  to  late-October  2001  and  is  evident  both  in  the 
satellite  ice  concentration  data  as  well  as  in  the  A2  ice-draft  time  series.  As  the  sea- 
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ice  retreats,  a  coherent  southward  along-shelf  flow  reappears,  and  is  only  interrupted 
in  mid-  to  late-December,  when  a  strong  reduction  of  the  along-shelf  velocity  was 
observed. 


5.5  Discussion 

The  freshwater  sources  which  force  the  APCC  are  difficult  to  quantify  due  to  the  lack 
of  reliable,  continuous  observations  of  freshwater  inputs  to  the  shelf.  We  will  attempt 
to  enumerate  and  quantify  these  sources  and  discuss  their  relative  importance  in  the 
context  of  a  freshwater  budget  for  a  coastal  region  upstream  of  Section  I.  We  finish 
with  a  brief  discussion  of  the  effects  of  wind  forcing  on  the  APCC. 

5.5.1  Fresh  water  budget 

Unlike  buoyant  plumes  formed  by  riverine  outflow  (e.g.,  Yankovsky  and  Chapman, 
1997;  Whitney  and  Garvine,  2005;  Lentz  and  Largier,  2006),  where  most  of  the  fresh¬ 
water  can  be  traced  back  to  a  single  upstream  source,  the  inputs  for  the  APCC  are 
multiple  and  are  distributed  along  the  coast.  A  freshwater  budget  for  the  APCC 
upstream  of  Section  I  can  be  written  as 

FWT  =  U+{R)  +  (P-  E) ocean  +  (Mice)  +  (O),  (5.3) 

where  FWT  is  the  freshwater  transport  measured  at  Section  I  and  U  is  the  fresh¬ 
water  transport  at  some  point  upstream  of  Section  I.  The  other  terms  refer  to  the 
influx  of  freshwater  between  that  upstream  point  and  Section  I:  R  is  the  runoff  from 
land,  P  —  Eocean  is  the  precipitation  minus  evaporation  that  occurs  over  the  ocean, 
and  Mice  is  the  sea-ice  melt.  Other  sources  or  sinks  of  freshwater,  such  as  cross-shelf 
Ekman  transport,  eddies,  etc.  are  represented  by  O  and  are  not  considered  further, 
as  we  have  no  information  about  them.  Integration  of  the  terms  over  a  given  region 
upstream  of  Section  I  is  represented  by  brackets  (()).  Because  there  are  no  simul¬ 
taneous  observations  of  the  terms  on  the  right-hand-side  of  the  freshwater  budget 
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above,  we  assume  instead  that  FWT  from  Section  I  is  a  reasonable  estimation  of  the 
climatological  freshwater  transport  of  the  APCC,  and  use  both  published  data  and 
results  from  an  atmospheric  regional  model  to  estimate  the  right-hand-side  terms  of 
Eq.  5.3. 

Evidence  for  the  APCC  upstream  of  Section  I 

A  brief  review  of  the  available  literature  on  the  circulation  north  of  Section  I  can 
shed  light  on  the  approximate  size  of  the  upstream  region  influencing  the  nearshore 
freshwater  budget  at  Section  I,  and  the  size  of  U .  Near  the  tip  of  the  Antarctic 
Peninsula,  the  Antarctic  Coastal  Current  (Fahrbach  et  al.,  1992)  has  been  observed  to 
flow  northwards  and  to  continue  counter-clockwise  around  the  tip  and  then  westwards 
into  the  Bransfield  Strait  (Heywood  et  al.,  2004).  Carrying  water  with  T/S  properties 
similar  to  those  of  the  APCC,  the  Antarctic  Coastal  Current  was  found  to  transport 
~1  Sv  into  Bransfield  Strait,  although  its  path  after  it  entered  the  strait  could  not 
be  determined  (Heywood  et  al.,  2004). 

Studies  of  the  surface  circulation  in  the  Bransfield  and  Gerlache  Straits,  using 
drifter  observations  and  hydrographic  data,  suggest  the  dominant  feature  in  this  re¬ 
gion  is  flow  crossing  from  the  northern  section  of  Gerlache  Strait  to  the  southern  coast 
of  the  South  Shetland  Islands,  although  meanders  and  recirculations  are  common  else¬ 
where  in  this  area  (Niiler  et  al.,  1990;  Huntley  and  Niiler,  1995;  Zhou  et  al.,  2002). 
This  northward  circulation  from  Gerlache  Strait  to  Bransfield  Strait  suggest  that  the 
path  of  the  Antarctic  Coastal  Current  might  retroflect  within  this  region  rather  than 
continue  along  the  shelf  and  towards  the  SO  GLOBEC  study  region.  This  is  consis¬ 
tent  with  summertime  observations  from  the  Palmer  Long  Term  Ecological  Research 
(LTER)  program,  which  show  minimal  amounts  of  mixing  between  Bransfield  Strait 
waters  and  the  waters  farther  south  (Capella  et  al.,  1992;  Smith  et  al.,  1999). 

For  the  region  south  of  Bransfield  Strait,  a  few  published  observations  show  a 
southward  circulation  between  Anvers  Island  and  Section  I.  A  drifter  released  in  Ger¬ 
lache  Strait  on  December  1991  moved  to  the  southwest  along  the  coast  (Huntley  and 
Niiler,  1995),  and  another  drifter  passing  by  Anvers  Island  also  moved  to  the  south- 
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west  in  the  last  few  months  of  1979  (Hofmann  et  ah,  1996).  This  is  consistent  with 
LTER  observations,  which  showed  the  presence  of  a  cross-shelf  density  gradient  and 
geostrophic  flow  to  the  southwest  near  the  coast  south  of  Anvers  Island  (Smith  et  ah, 
1999).  In  summary,  there  is  little  observational  evidence  of  a  continuous  southward 
circulation  from  Bransfield  Strait  to  the  shelf  south  of  Anvers  Island,  suggesting  that 
the  APCC  has  a  relatively  local  origin  somewhere  between  Anvers  Island  and  Section 
I.  Thus,  we  take  U  to  be  zero  and  consider  the  freshwater  budget 

FWT  =  (R)  +  (P-  E) ocean  +  (Mice)  (5.4) 

for  the  region  from  Anvers  Island  to  Section  I  (Figure  5-13). 

Runoff  from  land 

The  mass  budget  of  the  land  ice  on  the  Antarctic  Peninsula  is  given  by 

^-  =  (P-E)land-(R)-(C),  (5.5) 

where  M  is  the  mass  of  the  land  ice,  (P—E)iand  is  the  precipitation  minus  evaporation 
over  the  land  ice,  R  is  the  net  runoff,  and  C  represents  iceberg  calving  into  the  ocean. 
Although  there  is  no  quantitative  information  about  the  size  of  C,  icebergs  that 
melt  locally  are  a  net  input  of  freshwater  to  the  nearshore  ocean.  Lacking  better 
information,  we  assume  R»  C  (or  equivalently,  we  assume  icebergs  melt  primarily 
near  the  coast),  and  we  estimate  the  runoff  from  land  as 

(R)  =  (P-E)land-~.  (5.6) 

Thus,  net  runoff  from  land  is  understood  here  to  represent  the  net  precipitation  minus 
evaporation  over  land  ((P  —  E)iand)  minus  long-term  land  ice  change  dM/dt.  Note 
that  downstream  of  Section  I,  the  George  VI  ice  shelf  is  expected  to  make  a  significant 
contribution  to  the  budget,  as  studies  suggest  there  are  53  km  yr-1  of  basal  melt  for 
the  entire  ice  shelf  (Potter  and  Paren,  1985),  and  therefore  the  runoff  term  would 
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include  such  contribution. 

Both  the  mountainous  topography  of  the  peninsula  and  a  consistently  westerly  at¬ 
mospheric  flow  contribute  to  the  wAP  having  the  largest  precipitation  rate  in  Antarc¬ 
tica  (van  Lipzig  et  al.,  2004).  Since  in  situ  measurements  are  sparse,  van  Lipzig  et  al. 
(2004)  used  a  regional  atmospheric  mesoscale  model  to  quantify  the  spatial  variability 
in  surface  precipitation  minus  evaporation  (P  —  E).  This  Regional  Atmospheric  Cli¬ 
mate  MOdel  (RACM014),  using  a  horizontal  grid  spacing  of  14  km,  was  integrated 
for  the  period  1987-1993.  The  annually  averaged  model  output  (Figure  5-13)  shows 
P  —  E  values  in  excess  of  2  m  yr-1  along  most  of  the  west  drainage  basin  of  the 
peninsula  north  of  Marguerite  Bay.  In  the  SO  GLOBEC  study  region,  high  values 
are  found  along  the  west  coast  of  Adelaide  Island,  the  north-east  coast  of  Marguerite 
Bay,  Wordie  Bay  and  the  northern  coast  of  Alexander  Island. 

In  order  to  calculate  the  annually  averaged  input  of  freshwater  to  the  ocean  from 
the  land  (P  —  E)iand ,  we  need  to  define  the  eastern  boundary  of  the  region  between 
Anvers  Island  and  Section  I  where  the  integration  will  be  carried  out.  We  used  the 
RACM014  topography  to  find  the  spine  of  the  peninsula,  i.e.  the  limit  of  the  western 
drainage  basin  (Figure  5-13).  The  integration  of  the  7-year  average  P  —  E  over  the 
land,  integrated  from  Anvers  Island  to  Section  I  and  from  the  spine  of  the  peninsula 
to  the  coast,  results  in  a  influx  of  68±20  km3yr-1  (±1  standard  deviation).  Here, 
we  assume  the  freshwater  accumulated  over  land  during  the  winter  is  discharged  to 
the  ocean  preferentially  during  the  ice-free  season,  as  is  suggested  by  the  results  of 
Meredith  et  al..  The  uncertainty  above  includes  only  one  standard  deviation  of  the 
monthly  averaged  model  output,  and  not  the  differences  with  observations  or  other 
sources  of  error. 

The  Antarctic  Peninsula  is  one  of  the  fastest-warming  regions  in  the  world,  with 
mean  surface  atmospheric  warming  of  2.5  to  3°C  since  1951,  accompanied  by  a  1°C 
warming  of  the  ocean’s  surface  (King,  1994;  Vaughan  and  Doake,  1996;  Meredith  and 
King,  2005).  Of  the  244  glaciers  draining  the  Antarctic  Peninsula  ice  sheet,  87%  have 
retreated  in  the  period  between  1940  and  2004,  and  retreat  rates  are  accelerating 
(Cook  et  al.,  2005).  Although  few  published  estimates  exist  of  the  secular  freshwater 
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flux  resulting  from  this  process,  an  extrapolation  of  the  retreat  rates  of  glaciers  melting 
into  Wordie  Bay  (a  small  bay  within  Marguerite  Bay,  Figure  5-1)  to  a  region  including 
the  west  basin  of  the  peninsula  from  its  tip  to  Wordie  Bay  suggests  20±6  km3  of 
freshwater  may  be  lost  to  the  ocean  every  year  (Rignot  et  al.,  2005).  A  separate 
study  estimated  the  runoff  as  a  result  of  the  increasing  number  of  melting  days  on 
the  peninsula  to  be  between  2.9  to  20  km3  yr-1  (For  comparison  purposes,  we  convert 
this  range  to  an  average  of  11.5±8.5  km3yr_1).  This  estimate  included  the  entire 
peninsula  (including  west  and  east  drainage  basins)  for  the  year  2000  (Vaughan, 
2006).  Assuming  the  same  rates  for  the  ~  19200  km2  region  used  in  our  budget,  and 
averaging  the  estimates  of  the  two  studies  cited,  we  obtain  a  dM/dt  of  -6±4  km3  yT-1. 
Adding  the  estimate  of  (P  —  E)land  north  of  Section  I  and  the  size  of  dM/dt  discussed 
above,  our  estimate  for  ( R )  is  therefore  74±20  km3  yr"1. 

Precipitation  minus  Evaporation  over  the  Ocean 

Uplift  of  air  occurs  upstream  of  the  mountains  of  the  Antarctic  Peninsula,  resulting 
in  high  precipitation  rates  over  the  nearshore  ocean  of  the  wAP  (van  Lipzig  et  al., 
2004),  as  illustrated  in  Figure  5-13.  Cross-shelf  profiles  of  P  —  E  upstream  of  Section 
I  suggests  this  input  decays  from  the  coast  with  an  offshore  scale  of  «25  km.  In 
order  to  quantify  the  importance  of  this  component,  the  P  —  E  field  over  the  ocean 
is  integrated  from  Anvers  Island  to  Section  I  and  from  the  coast  offshore  to  the 
1  m  yr-1  P  —  E  isopleth  (some  20-50  km  off  the  coast).  This  results  in  27±7  km3  yr-1 
of  freshwater  input  to  the  nearshore  ocean  from  (P  —  E)ocean.  This  integration  only 
includes  the  influx  during  the  seven  ice-free  months  of  the  year  (from  November  to 
June),  and  not  during  the  winter  months  (July  to  October),  consistent  with  Section 
5.5.1  (see  below)  and  the  freshwater  transport  calculation  at  Section  I. 

Sea  Ice  Melt 

Sea-ice  is  formed  at  the  end  of  the  fall  and  throughout  the  winter,  creating  a  sea¬ 
sonal  cover  that  slowly  retreats  during  the  spring  (Stammerjohn  et  al.,  2003).  Sea- ice 
samples  collected  during  the  2001  SO  GLOBEC  winter  cruise  had  an  average  thick- 
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ness  of  62±44  cm  (Perovich  et  al.,  2004).  Relative  to  other  places  in  the  Southern 
Ocean,  the  wAP  shelf  region  has  high  precipitation  rates  and  relatively  thin  sea-ice 
(Perovich  et  al.,  2004;  van  Lipzig  et  al.,  2004).  This  leads  to  proportionally  deep  snow 
and  negative  sea-ice  freeboard,  setting  the  conditions  for  significant  surface  flooding, 
refreezing  and  snow-ice  formation. 

Quantifying  the  contribution  of  sea-ice  melt  Mlce  to  the  freshwater  budget  of  the 
APCC  is  particularly  difficult.  Ice  formed  from  sea-water  and  melted  in  the  same 
location  does  not  generate,  over  a  yearly  cycle,  a  net  input  of  freshwater.  Other 
dynamics,  however,  are  potentially  important.  Below,  we  discuss  two  mechanisms 
to  input  freshwater  near  the  coast  related  to  sea-ice:  shoreward  sea-ice  drift  and 
precipitation  near  the  shore  stored  throughout  the  winter  as  snow  and  snow-ice. 

Evidence  for  sea-ice  drift  has  been  reported  in  several  regional  studies  of  sea- 
ice  dynamics.  The  meridional  atmospheric  circulation  plays  an  important  role  in 
determining  the  extent  and  composition  of  the  sea-ice  cover,  both  by  advecting  the 
existing  sea-ice  and  opening  and  closing  leads  where  sea-ice  production  is  enhanced 
(Stammerjohn  et  al.,  2003;  Massom  et  al.,  2006;  Harangozo,  2006).  In  nearshore 
areas,  wind  forcing  can  result  in  a  significant  amount  of  compaction  and  rafting  of 
the  sea-ice.  This  process  was  especially  dramatic  during  the  2001-2002  season,  when 
anomalously  strong  north-westerly  winds  were  reported  to  cause  rafting  of  sea-ice 
near  the  shore,  reaching  a  thickness  of  10  to  15  m  during  an  October  cruise  in  mid- 
Marguerite  Bay  (Massom  et  al.,  2006). 

High  snowfall  rates,  another  source  of  freshwater  input,  are  also  evident  in  the 
SO  GLOBEC  sea-ice  measurements.  A  mass  balance  buoy  was  deployed  on  August 
6,  2001,  on  sea-ice  initially  0.65  m  thick  with  0.25  m  of  snow  on  top.  The  buoy 
functioned  through  November  10,  mostly  within  100  km  or  so  from  the  northern  shore 
of  Alexander  Island  and  measured  an  increase  of  snow  thickness  of  approximately  1  m 
(see  Figure  10  in  Perovich  et  al.,  2004).  Although  the  accumulation  for  the  rest  of 
the  year  is  unknown  and  the  density  of  the  snow  cover  was  not  recorded,  these  rates 
of  snowfall  are  an  indication  that  the  rates  of  1  to  2.5  m  yr'1  of  water  equivalent 
suggested  by  the  atmospheric  model’s  climatology  for  approximately  the  same  region 
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are  reasonable  (Figure  5-13).  Using  the  same  region  as  in  Section  5.5.1,  the  model 
output  results  in  37±12  km3yr~1  of  mean  freshwater  accumulation  over  the  winter 
months  (July  to  October). 

As  the  sea- ice  retreats  and  melts  throughout  the  spring  (and  into  summer),  the 
freshwater  input  from  both  shoreward  advection  of  ice  and  precipitation  will  con¬ 
tribute  to  the  observed  APCC  freshwater  transport.  The  resulting  meltwater  will 
presumably  leave  the  region  of  interest.  An  estimate  of  the  timescale  for  this  signal 
to  propagate  alongshore  away  from  the  region  can  be  obtained  by  estimating  the 
propagation  speed  of  a  buoyant  plume  of  finite  along-shelf  length  trapped  against  the 
coast.  In  the  absence  of  any  other  forcing,  the  leading  and  trailing  edges  (the  “nose”) 
of  the  plume  will  propagate  at  a  speed  given  by: 

^~TT WF.'  (5'7) 

where  Cp  is  the  nose  speed  of  the  plume  and  Cu,  =  (2 g'Qf)1^4  is  the  nose  speed  of  a 
plume  if  it  were  moving  along  a  vertical  wall  (Lentz  and  Helfrich,  2002).  The  ratio 
Wb/Ww  (Figure  5-10)  describes  the  importance  of  a  bottom  slope  in  the  dynamics, 
with  surface-trapped  plumes  having  ratios  closer  to  zero.  Using  Wb/Ww  from  0.1  to 
10  and  the  observed  transport  and  density  anomaly  from  Section  I,  Eq.  5.7  gives 
Cp  of  ~0.7  m  s-1  to  ~  0.07  m  s-1,  respectively.  With  this  scaling,  which  assumes 
the  density  anomaly  generated  by  sea-ice  melting  would  be  as  large  as  the  observed 
one,  the  freshwater  signal  would  propagate  through  the  ~270-km  long  stretch  from 
Anvers  Island  to  Adelaide  Island  in  5  days  to  6  weeks.  Even  at  the  slowest  estimate, 
it  would  seem  that  sea-ice,  which  decays  to  a  minimum  extent  during  the  summer 
(Stammerjohn  et  al.,  2003),  would  have  a  maximum  impact  on  the  freshwater  budget 
considered  here  during  the  spring  and  early  summer,  and  not  during  the  fall,  when 
the  data  discussed  here  were  taken. 

A  few  other  datasets  address  the  relevance  of  sea-ice  to  the  nearshore  freshwater 
budget.  Hydrographic  data  taken  off  Anvers  Island  from  1991  to  1999  suggest  the 
summer  nearshore  freshwater  is  largely  explained  by  runoff  detected  as  pulses  of  fresh- 
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water  in  nearshore  stations  throughout  the  summer,  and  not  melted  sea-ice  (Dierssen 
et  al.,  2002).  Also,  <5180  measurements  collected  close  to  the  British  Antarctic  Survey 
Rothera  Station  in  Marguerite  Bay  during  2002  and  2003  were  used  to  quantify  the 
relative  contributions  of  sea-ice  and  freshwater  runoff  from  land  to  the  local  upper- 
column  freshwater  content.  These  observations  suggest  runoff  from  land  is  largely 
responsible  for  the  freshwater  found  in  the  surface  layer,  and  that  this  runoff  has  a 
significant  seasonal  signal,  with  maximum  input  from  February  to  June  (Meredith 
et  al.).  However,  it  must  be  noted  that  the  oxygen  isotope  samples  were  collected  on 
the  leeward  side  of  Adelaide  Island,  so  snow  and  snow-ice  accumulation  over  the  sea 
ice  is  not  expected  to  be  as  large  as  on  the  windward  side,  where  Section  I  is  located 
(Figure  5-13). 

In  summary,  we  expect  sea-ice  melt  to  play  a  potentially  important  role  in  the 
freshwater  budget  of  the  APCC,  in  particular  as  the  sea-ice  melts  during  the  spring 
and  the  37±12  km3yr-1  that  the  model  output  suggest  have  accumulated  on  the  ice 
start  melting.  The  details  of  how  this  freshwater  is  discharged  to  the  ocean,  and 
quantification  of  the  sea- ice  drift  contribution  are  difficult  to  determine  at  this  point, 
but  sea- ice  is  likely  not  relevant  to  the  discussion  of  the  freshwater  budget  when 
using  observations  taken  in  June.  At  the  end  of  the  fall,  we  expect  most  of  the  sea- 
ice  to  have  melted  and  the  associated  freshwater  anomaly  to  have  propagated  away 
from  the  region  considered  here.  We  therefore  assume  this  component  to  be  zero, 
with  the  important  caveat  that  our  reasoning  here  must  be  tested  using  new  field 
measurements  and  high-resolution  models  of  this  region. 

Total  Freshwater  Input  Upstream  of  Section  I 

The  sum  of  runoff  from  land  ((P  -  E)iand  and  ( dM/dt ))  and  precipitation  over  the 
ocean  ((P  —  E)ocean )  could  contribute  101±22  km3yr-1  of  freshwater  to  the  coastal 
ocean  from  Anvers  Island  to  Section  I  (Eq.  5.3),  with  an  unknown  contribution  from 
sea-ice  and  other  sources  and  sinks.  Given  our  estimates  of  the  freshwater  upstream 
of  Section  I  (126±50  km3  yr-1  -  Section  5.4.2),  this  suggests  ( R )  and  (P  —  E)ocean  are 
the  primary  sources  of  freshwater  for  the  APCC.  In  the  previous  section,  we  argued 
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that  sea-ice  is  of  secondary  importance  to  the  forcing  of  the  APCC  during  summer 
and  fall.  The  freshwater  budget  is  certainly  affected  by  sea-ice  formation  and  melting, 
although  near  the  coast,  it  is  the  runoff  from  land  that  appears  to  play  a  dominant 
role  (Meredith  et  al.,  this  issue).  In  this  scenario,  the  driving  factor  is  the  reduction  of 
air  temperature  as  winter  arrives.  This  causes  both  sea-ice  formation  and  reduction 
of  freshwater  input  from  the  land,  but  it  is  the  latter  that  would  be  responsible  for 
the  shutdown  of  the  APCC  (Section  5.4.4).  However,  our  estimation  of  the  transport 
(based  on  a  single  shipboard  ADCP  section),  the  sources  for  our  estimation  of  the 
freshwater  influx,  and  our  knowledge  of  the  local  sea-ice  distribution  and  dynamics  are 
subject  to  significant  uncertainty,  so  additional  observations  are  necessary  to  better 
determine  the  relative  contributions  of  the  different  freshwater  sources  to  the  forcing 
of  the  APCC. 

5.5.2  Wind  Effects 

While  on  June  1st,  2001  when  the  data  in  Section  I  were  collected,  the  winds  were 
rather  weak  (<2  m  s_1),  strong  downwelling-favorable  winds  are  prevalent  in  this  re¬ 
gion  near  the  coast.  They  are  expected  to  significantly  affect  the  structure  and  dynam¬ 
ics  of  the  APCC.  Although  the  response  of  coastal  buoyant  currents  to  downwelling- 
favorable  winds  is  not  well  understood,  it  is  known  that  it  will  force  a  cross-shelf 
Ekman  transport  near  the  surface  and  an  additional  along-shelf,  geostrophically  bal¬ 
anced  flow  that  can  increase  the  observed  transport  of  the  current  (Fong  and  Geyer, 
2001;  Whitney  and  Garvine,  2005;  Lentz  and  Largier,  2006). 

Assuming  a  simple  model  with  linear  bottom  stress  balancing  the  wind  stress,  the 
wind-driven  along-shelf  transport  is  Qw  =  r/(p0r)A,  where  r  is  the  along-shelf  wind 
stress,  p  is  a  reference  density,  A  is  the  cross-sectional  area  of  the  plume  and  r  is  a 
bottom  drag  coefficient,  taken  here  as  5  10~4  m  s-1  (Lentz  and  Largier,  2006).  Wind 
measurements  collected  in  Marguerite  Bay  (Moffat  et  al.,  2005)  from  June  of  2001  to 
April  of  2003  have  a  mean  along-shelf  wind  stress  of  -0.04±0.14  N  m-2  (±1  standard 
deviation),  and  a  maximum  value  of  -1.71  N  m-2  (Here,  the  coordinate  system  is  the 
same  used  for  Section  I,  so  negative  values  indicate  downwelling-favorable  winds). 
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Over  the  same  area  enclosed  by  the  34.4  isohaline  and  a  24-km  offshore  limit,  which 
resulted  in  a  volume  transport  of  0.32  Sv  (Section  5.4.1),  Qw  is  0.28  Sv  for  the  mean 
wind  stress.  Although  this  model  is  very  simple,  it  illustrates  the  sensitivity  of  our 
transport  calculations  to  wind  forcing,  and  highlights  the  importance  of  concurrent 
atmospheric  measurements  on  the  wAP  shelf  when  considering  the  dynamics  of  the 
APCC. 

5.6  Summary  and  Open  Questions 

Observations  collected  during  2001-2002  as  part  of  the  Southern  Ocean  GLOBEC 
program  reveal  the  presence  of  a  buoyancy-forced,  geostrophically-balanced  coastal 
current  along  the  west  coast  of  the  Antarctic  Peninsula,  from  Adelaide  Island  to 
Alexander  Island.  The  Antarctic  Peninsula  Coastal  Current  (APCC)  is  strongly 
seasonal  and  forms  during  the  spring/summer  when  the  shelf  is  ice- free.  The  APCC 
is  characterized  by  a  surface  density  anomaly  of  ~0.6  kg  m-3  and  surface  velocities  of 
the  order  of  0.30-0.40  m  s-1.  One  high-resolution  hydrographic  and  velocity  section 
off  the  coast  of  Adelaide  Island  suggests  the  volume  transport  of  the  APCC  is  of 
order  0.32db0.13  Sv  with  a  freshwater  transport  of  order  126±50  km3yr_1.  Three 
additional  ADCP  sections  along  the  coast  show  a  velocity  front  consistent  with  a 
coastal  buoyant  current  off  Adelaide  and  Alexander  Islands  (Figure  5-14). 

Moored  observations  show  the  virtual  disappearance  of  the  APCC  during  the 
winter  months,  after  sea-ice  forms  on  the  shelf  and  the  melt- water  fluxes  from  the  coast 
are  drastically  reduced.  As  the  sea-ice  breaks  up  and  melts  with  increasing  insolation 
and  air  temperature,  the  APCC  reappears.  Output  from  a  regional  atmospheric 
model  and  a  review  of  the  literature  suggest  runoff  from  land  and  precipitation  over 
the  ocean,  when  integrated  from  Anvers  Island  south,  can  provide  roughly  enough 
freshwater  to  the  nearshore  ocean  to  explain  the  observed  freshwater  transport  of 
the  APCC  off  Adelaide  Island.  Other  sources  and  sinks,  like  ice-melt,  wind-driven 
Ekman  transport  and  eddies  were  not  considered,  but  may  play  an  important  role. 
Previous  studies  reviewed  here  suggest  the  wAP  region  presents  particularly  favorable 
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conditions  (high  P  —  E,  retreating  glaciers,  downwelling-favorable  winds)  for  the 
formation  of  a  coherent  seasonal  coastal  buoyant  current. 

Many  questions  remain.  Although  it  is  likely  that  the  continuous  freshwater  influx 
along  the  peninsula  makes  the  APC’C  a  ubiquitous  feature  along  the  coast  of  the  wAP 
south  of  Anvers  Island,  more  detailed  observations  are  required  to  determine  its  path 
in  the  SO  GLOBEC  study  region.  Of  particular  interest  is  the  path  (or  paths)  of  the 
APCC  inside  Marguerite  Bay.  which  could  not  be  determined  with  our  data  set. 

An  important  point  that  needs  to  be  investigated  is  the  interannual  variability  of 
the  APCC.  Researchers  have  found  significant  interannual  variability  in  upper-ocean 
characteristics  (Dierssen  et,  al.,  2002;  Meredith  et  al.,  2004)  and  sea-ice  (Stainmerjohn 
et  al..  2003;  Perovich  et  al.,  2004;  Harangozo,  2006)  among  others,  and  much  of  that 
variability  is  likely  to  directly  impact  the  APCC.  Observations  designed  to  monitor 
its  strength  and  along-shelf  variability  seem  necessary  to  gauge  the  overall  importance 
of  the  APCC  on  the  wAP  system. 

Also,  we  have  not  addressed  the  connections  of  the  APCC  to  the  rich  biological 
communities  present  on  the  wAP.  Several  recent,  studies  of  the  biological  communit  ies 
in  this  region  have  observed  strong  cross-shelf  gradients  in  biogeochemical  variables. 
Chlorophyll  (Dierssen  et  al.,  2002),  nutrient  deficits  (Serebrennikova  and  Fanning, 
2004)  and  plankton  concentrations  (Ashjian  et  al.,  2004;  Lawson  et  al.,  2004;  Law- 
son,  2006)  are  high  along  the  coast.  This  suggests  that  the  APCC  might  be  key  to 
providing  a  favorable  environment  for  biological  production.  This  might  be  due  to  a 
variety  of  reasons:  a  coastal  buoyant  current  like  the  APCC  retains  freshwater  and 
nutrients  (like  trace  metals)  contained  therein  near  the  shore,  it  has  relatively  high 
stratification,  and  could  be  host  to  a  biological  community  well  adapted  to  the  cold 
and  low-salinity  water  maintained  by  the  current.  Also,  the  APCC  provides  a  path 
for  the  along-shelf  transport  of  biogeochemical  tracers  during  the  ice-free  season.  The 
ongoing  synthesis  effort,  of  the  SO  GLOBEC  data  sets  and  further  research  into  these 
and  other  questions  should  help  address  some  of  the  many  issues  raised  here. 
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Figure  5-1:  The  west  Antarctic  Peninsula  (vvAP)  eontintental  shelf.  The  gray  lines 
indicate  the  500-  and  1000-m  isobath,  and  the  permanent  ice-shelf  cover  is  indicated 
by  light  blue.  Also  shown  are  the  A,  B  and  C  CTD  lines  visited  during  the  fall  and 
winter  survey  cruises  of  2001,  the  high-resolution  Section  I,  the  locations  of  the  2001 
Al  and  A2  Moorings  and  of  the  British  Antarctic  Survey  Rothera  Station. 
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Figure  5-2:  Schematic  of  the  Al  mooring,  deployed  during  2001  and  early  2002  oil 
Adelaide  Island.  One  of  the  VAC  Ms  had  a  temperature  recorder  (VACM/T)  and  the 
other  a  pressure  recorder  (VACM/P).  The  dashed  line  shows  a  profile  of  the  cross-shelf 
bathymetry  ~3  kin  southwest  of  the  mooring  site. 
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Figure  5-3:  Evolution  of  the  hydrographic  structure  on  the  mid-shelf  of  the  wAP  (A2 
mooring  site).  The  left  panel  shows  a  6- S  diagram  of  profiles  taken  throughout  2001 
and  2002.  The  gray  boxes  illustrate  the  typical  range  of  the  Antarctic  Surface  Water 
(AASW),  Winter  Water  (WW)  and  Circumpolar  Deep  Water  (CDW).  The  dashed 
line  indicates  the  freezing  point  of  seawater.  The  right  panels  show  the  corresponding 
potential  temperature  (0),  salinity,  and  potential  density  (ae)  profiles. 
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Figure  5-4:  Sea  surface  salinity  (SSS)  and  sea  surface  temperature  (SST)  for  the  fall 
(CTD  and  underway  data,  top  panels)  and  winter  (CTD  data,  bottom  panels)  of 
2001,  collected  during  the  SO  GLOBEC  broad-scale  cruises  on  board  the  RV/IR  N. 
B.  Palmer. 
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Figure  5-5:  Potential  temperature-Salinity  (0- S)  diagrams  from  line  A  (top  panels), 
B  (middle  panels)  and  C  (bottom  panels)  for  the  fall  of  2001  (left  column)  and  winter 
of  2001  (right  column).  The  solid  black  lines  are  stations  from  the  mid-  and  outer 
shelf,  and  the  dashed  red  lines  are  nearshore  stations.  The  dashed  line  in  each  panel 
indicates  the  freezing  point  of  sea- water.  Station  locations  are  plotted  in  Figure  5-1. 
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Figure  5-6:  Gridded  temperature,  salinity  and  along-shelf  velocity  (ADCP)  data  from 
cross-shelf  Section  I,  occupied  on  June  1st,  2001.  The  red  dots  indicate  the  approxi¬ 
mate  location  of  BIOMAPER-II  dives  (top  and  middle  panels)  and  the  ADCP  5-min 
averages  (bottom  panel).  The  thick  black  line  in  the  middle  panel  indicate  the  34.4 
isohaline.  The  dashed  black  line  in  the  bottom  panel  indicates  the  approximate  cross¬ 
shelf  position  of  the  Al  mooring.  The  along-shelf  velocity  is  defined  as  the  component 
perpendicular  to  the  section  shown  in  Figure  5-1,  with  negatives  values  denoting  flow 
to  the  southwest. 
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Figure  5-7:  Mean  ADCP  velocity  profile  for  the  velocity  front  depicted  in  Figure 
5-6  (solid)  and  mean  velocity  calculated  from  thermal  wind  shear  (with  an  arbitrary 
offset)  over  the  same  cross-shelf  region  (4  to  32  km  from  the  coast). 
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Figure  5-8:  Cross-shore  ADCP  sections  collected  during  the  fall  2001  SO  GLOBEC 
Cruise.  The  lines  indicate  the  location  of  Section  I.  and  sections  conducted  off  the 
southern  coast  of  Adelaide  Island  (II),  and  off  the  northern  (III)  and  western  (IV) 
coast  of  Alexander  Island. 
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Figure  5-9:  Along-shelf  velocity  ADCP  sections,  with  negative  values  indicating  flow 
with  the  coast  to  the  left.  The  locations  of  Section  II  (top  Panel),  III  (middle  panel) 
and  IV  (bottom  panel)  are  shown  in  Figure  5-8.  The  thin  solid  lines  near  the  topog¬ 
raphy  delineate  the  boundary  for  the  velocity  estimates. 
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Figure  5-10:  Schematic  of  the  relevant  scales  of  a  coastal  buoyant  plume.  The  plume 
front,  of  width  Ws,  separates  the  shelf  water  of  density  p0  from  the  lighter  (by  A p) 
nearshore  water.  The  front  depth  is  h  and  the  front  touches  the  bottom  at  hp  and  a 
distance  from  the  coast  114  (modified  from  Lentz  and  Largier,  2006). 
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Figure  5-11:  Time  series  of  along-shelf  (positive  north-eastwards)  velocity  at  the  Al 
mooring  (68°  01.217'  W,  67°  01.134'  S).  The  colored  lines  show  velocity  at  different 
depths  from  the  upward-looking  ADCP  (65  to  165  m)  and  VACMs  (309  and  459 
in).  The  shaded  areas  indicate  periods  when  the  fall  and  winter  broad-scale  cruises 
were  conducted.  The  black  line  shows  the  sea-ice  concentration  (in  percent)  for  the 
nearest  grid  point  to  the  Al  mooring  site  from  the  Special  Sensor  Microwave /Imager 
(SSM/I)  Satellite  data  (Comiso,  1999.  updated  2005).  The  magenta  line  shows  the 
sea-ice  draft  (in  meters)  measured  with  the  IPS  at  the  A2  mooring. 
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Figure  5-12:  Average  of  along-shelf  velocities  from  the  Al  mooring  records  for  the 
fall  (February  to  June),  winter  (here  defined  as  the  ice-covered  season,  from  July 
to  October)  and  spring/summer  (From  November  to  the  end  of  the  deployment,  in 
February)  shown  in  Figure  5-11.  The  lines  indicate  the  mean  values  for  the  season, 
and  the  shading  the  standard  error. 
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Figure  5-13:  Modeled  precipitation  minus  evaporation  ( P  -  E)  map  for  the  wAP 
region  (data  from  van  Lipzig  et  al.,  2004).  The  thick  gray  line  indicate  the  location 
of  the  spine  of  the  peninsula  (from  the  RACM014  topography),  and  the  height  of 
the  topography  along  the  spine  is  plotted  on  the  left  panel.  The  dashed  black  line 
indicates  the  region  upstream  of  Section  I  used  to  compute  the  freshwater  budget- 
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Figure  5-14:  Schematic  of  the  path  of  the  APCC.  The  solid  lines  indicate  sections 
along  the  coast  where  the  velocity  front  and  a  cross-shelf  density  gradient  was  ob¬ 
served.  The  dashed  lines  are  suggested  paths,  not  derived  from  observations. 
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Table  5.1:  Summary  of  the  cruise  dates  and  number  of  CTD  stations  conducted,  and 
of  the  dates,  location  and  local  depth  of  the  moorings  deployment. 


Cruise 

Objective 

Start 

End 

#  CTD  Stations 

LMG01-03 

Mooring  Deployment 

03/18/01 

04/13/01 

6 

NBP01-03 

Fall  2001  Survey 

04/24/01 

06/05/01 

84 

NBP01-04 

Winter  2001  Survey 

07/24/01 

08/31/01 

70 

LMG02-1A 

Mooring  Recovery 

02/06/02 

03/03/02 

9 

Mooring 

Deployment 

Recovery 

Longitude 

Latitude 

Depth  [m 

Al 

03/26/01 

02/13/02 

69°01.217' 

67°01.134' 

509 

A2 

03/31/01 

02/13/02 

70°00.683" 

66°5 1.883’ 

561 
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Chapter  6 


Discussion  and  Summary 


Research  presented  in  this  thesis  has  substantially  increased  our  knowledge  of  the 
circulation  and  dynamics  on  the  western  Antarctic  Peninsula  continental  shelf,  the 
study  site  of  the  Southern  Ocean  GLOBEC  program. 

In  Chapter  3,  the  basic  properties  of  the  circulation  are  described.  The  results  sug¬ 
gest  that  the  correlation  scales  are  relatively  large  for  tidal  and  near-inertial  motions 
on  the  mid-  and  outer-shelf  off  Adelaide  Island,  while  across  the  mouth  of  Marguerite 
Bay,  where  the  bathymetry  and  the  coastline  are  significantly  more  complicated,  the 
scale  of  the  motion  appears  to  be  smaller  than  the  mooring  separation  (25  km  or  so). 
Also,  it  is  shown  that  Marguerite  Trough  has  a  strong  influence  on  the  circulation, 
determining  the  orientation  of  the  mean  flow,  steering  the  tides,  and  sustaining  vor- 
ticity  waves.  Circulation  in  the  trough  also  shows  a  significant  wind-driven  response, 
with  downwelling-  (upwelling-)  favorable  wind-stress  generating  a  anticyclonic  (cy¬ 
clonic)  response  in  the  northern  side  of  the  canyon  and  at  A3.  In  the  exploration  of 
the  warm- water  intrusion  hypothesis  (Chapter  4),  the  trough  also  plays  a  role  in  the 
cross-shelf  advection  of  heat. 

Chapter  4  revealed  some  surprising  results.  Before  this  research,  it  was  thought 
that  the  UCDW  intrusions  on  the  western  Antarctic  Peninsula  continental  shelf  were 
of  relatively  large  scale  and  that  they  happened  5  to  6  times  a  year  (Klinck  et  al., 
2004).  The  mooring  observations,  however,  show  otherwise.  The  warm  UCDW  in¬ 
trusions  are  short,  frequent  (four  per  month  at  a  fixed  site,  on  average)  and  have  an 
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apparent  horizontal  scale  of  3-5  km.  They  are  advected  in  Marguerite  Trough  by  the 
low-frequency  cyclonic  flow,  although  they  are  not  observed  in  the  mooring  near  the 
mouth  of  Marguerite  Trough.  It  is  unclear  whether  they  mix  away  before  reaching  the 
deeper  part  of  the  trough  in  Marguerite  Bay,  or  whether  the  Lagrangian  circulation 
is  such  that  the  intrusions  are  advected  elsewhere  on  the  shelf. 

Further  investigation  is  necessary  to  understand  the  underlying  dynamics  of  the 
intrusions.  No  attempt  was  made  here  to  explain  the  generation  process,  although 
their  properties  indicate  they  originate  in  the  ACC.  Also,  their  horizontal  scales  and 
velocity  structure  were  inferred  using  simple  dynamical  models  and  scaling  arguments, 
not  by  direct  observations.  A  moored  array  similar  to  the  one  used  by  SO  GLOBEC, 
but  consisting  of  more  closely  spaced  moorings  across  the  northern  edge  of  Marguerite 
Trough  (at  the  A2  site)  could  help  resolve  many  of  the  questions  left  open  by  the 
research  presented  here. 

Chapter  5  provided  the  first  description  of  the  Antarctic  Peninsula  Coastal  Current 
(APCC).  Although  the  uncertainties  are  large,  observations  and  modelling  results 
from  the  region  suggest  that  during  the  ice-free  season,  the  intense  precipitation  and 
melt  water  from  the  coast  results  in  the  formation  of  a  coastal  buoyant  current  in  the 
study  region.  The  current  is  roughly  in  geostrophic  balance,  and  is  the  most  coherent 
and  energetic  structure  observed  on  subtidal  time-scales. 

The  SO  GLOBEC  moored  array  has  provided  tremendous  amounts  of  new  infor¬ 
mation  about  the  circulation  on  the  wAP  shelf.  Some  basic  ideas  about  the  general 
circulation  and  scales  of  variability  are  now  in  place,  there  are  also  implications  from 
these  results  that  should  be  tested  more  extensively  with  new  observations. 

One  can  speculate  on  what  those  new  efforts  might  include.  On  the  structure  and 
origin  of  the  APCC,  an  important  question  is  what  is  the  origin  of  the  fresh  water, 
which  was  addressed  here  with  model  output  which  itself  needs  further  testing  with 
observations.  A  relatively  straightforward  way  to  make  progress  on  this  question 
would  be  to  make  extensive  measurements  of  oxygen  isotopes  during  the  ice-free 
season  on  the  shelf,  which  would  determine  the  relative  importance  of  sea-ice  vs 
glacier  melt- water  since  they  have  distinctive  isotopic  signatures  (e.g.  Meredith  et  al., 
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2006). 


The  APCC  is  also  interesting  in  that  our  scarce  observations  indicate  the  structure 
of  the  front  is  not  exactly  what  one  would  expect  from  the  theories  dealing  with  the 
evolution  of  coastal  buoyant  currents  on  the  shelf:  the  prediction  is  that  the  foot  of 
the  front  will  move  across  isobaths  (across  the  shelf)  until  the  near-bottom  velocity 
is  zero  (Lentz  and  Helfrich,  2002).  The  data  presented  here  shows  the  front  of  the 
APCC  to  be  located  several  kilometers  farther  offshore  than  expected  (i.e.,  there  is 
a  region  of  relatively  flat  isopycnals  between  the  coast  and  the  front).  This  suggest 
than  other  processes  (along-shelf  topographic  variations,  wind-forcing,  etc)  are  acting 
on  the  APCC  front,  apart  from  interactions  with  a  relatively  constant  bottom  slope. 

The  APCC  is  also  similar  to  other  coastal  buoyant  currents  in  the  northern  hemi¬ 
sphere,  like  the  East  Greenland  Coastal  Current  and  the  Alaskan  Coastal  current, 
that  are  relatively  large  (compared  to  their  small  counterparts  forced  by  river  out¬ 
flows  in  mid-latitudes)  and  subjected  to  strong  downwelling-favorable  wind  forcing 
and  relatively  continuous  buoyancy  inputs  from  the  coast.  All  of  these  currents  raise 
questions  about  the  influence  of  wind  forcing  on  deep  fronts  near  the  shore  which 
remain  relatively  unexplored. 

Both  the  APCC  and  deep  intrusions  form  fundamental  components  of  the  wAP 
ecosystem,  and  it  would  seem  that  they  explain  important  aspects  of  the  horizontal 
distribution  of  plankton  on  the  shelf.  Looking  at  the  relative  abundances  of  differ¬ 
ent  types  of  phytoplankton,  Prezelin  et  al.  (2004)  noticed  that  off  the  west  coast  of 
Adelaide  Island,  diatoms  dominated  over  phytoflagellates,  while  near  the  coast  there 
was  co-dominance  of  these  two  groups.  Diatoms  tend  to  dominate  in  environments 
where  there  is  sufficient  supply  of  silica-rich  water,  as  is  the  case  with  the  oceanic 
water  intrusions  studied  in  Chapter  4  (Prezelin  et  al.,  2000;  Serebrennikova  and  Fan¬ 
ning,  2004).  This  suggests  that  the  APCC  front,  which  separates  the  relatively  fresh, 
presumably  silica-poor  water  near  the  shore  from  the  regions  on  the  mid-shelf  more 
likely  to  receive  the  high-silica  intrusions,  might  create  a  favorable  environment  for 
phytoflagellates  near  the  coast,  by  limiting  the  supply  of  silica  to  the  diatom  popula¬ 
tion. 
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Overall,  what  makes  the  study  of  regions  like  the  wAP  challenging  is  the  combina¬ 
tion  of  a  deep  shelf  with  complicated  bathymetry,  and  a  relatively  weak  stratification. 
This,  combined  with  the  high  latitude  location  of  these  sites,  results  in  flows  with  rel¬ 
atively  small  dynamical  scales  following  complicated  paths  strongly  influenced  by  the 
bathymetry.  Future  observations  and  numerical  models  must  be  designed  to  resolve 
these  small  scales  in  order  to  understand  some  of  the  processes  outlined  in  this  thesis. 
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